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ABSTRACT 


This thesis consists of two parts: 


Part I: Densicy Matrix Simulation of Two-Dimensional and 
Multiple-Pulse Experiments 


A simulation program to analyze multiple-pulse and 
two-dimensional FT NMR experiments using density matrix 
theory has been implemented. The program computes 
analytical expressions for NMR signals arising from an 


LS) 


m°n #Pin system (where I and S are spin 1/2 nuclei and 


1 <n, m< 3), after the application of a pulse sequence 
specified by the input data. The simulation program has 
been designed to be quite general and incorporates a 
number of features including irradiation from a broadband 
decoupler, effects of phase cycling of the pulses and 
receiver, and pulses of any given flip angle. 

The state of a spin system during a multiple-pulse 
experiment is followed by monitoring the evolution of the 
density matrix. During each step in the sequence, the 
effective Hamiltonian operator will be time independent in 
a suitable rotating frame, and the equation of motion for 
the density matrix in this frame can be solved in a 
straightforward manner. The evolution of the density 
matrix through a pulse sequence is thereby reduced to a 


series of successive transformations on the density 
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matrix. The simulation program has been implemented in 
such a manner that the transformations of the density 
matrix, under radiofrequency pulses, free-precession 
periods, periods of broadband decoupling, etc., are 
effected by various independent program modules which have 
access to the density matrix data structure. 

A detailed investigation of the effects of ordinary 
noise decoupling and spherical randomization decoupling on 
the elements of the density matrix for 1S, spin systems 
has been carried out. It has been shown that under strong 
decoupler irradiation, the density matrix elements reach 
steady state values in the rotating frame of the decoupled 
nuclei. The steady state values are found to be linear 
combinations of the density matrix at the beginning of the 
decoupling period, and often involve mixing of populations 
with multiple quantum coherences, and mixing of the 
perpendicular components of the magnetization with higher 
coherences. The transformation of the density matrix 
during the decoupling period is carried out by the linear 
combination procedure using coefficients stored in the 
program. 

The application of the program in the analysis of 
multiple-pulse and two-dimensional experiments has been 
demonstrated by simulating the distortionless enhanced 


polarization transfer (DEPT) experiment, a two-dimensional 
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version of the DEPT experiment for multiplet-selective 
heteronuclear shift correlation, the INADEQUATE experiment 
and three variants of the gated decoupler 2D J-& resolved 
experiment which illustrate certain unexpected effects of 
spin decoupling. These examples demonstrate the 
usefulness of the program in simulating experiments in 
which the spin dynamics cannot always be followed using 


classical magnetization vector diagrams. 


Part 2: Nuclear Relaxation Study of HexaMethyl- 
Phosphoramide 


The spin-lattice relaxation times of 1 


31p, 144 and 


H 
in hexamethylphosphoramide (HMPA) and of 315 ana 7H in the 
deuterated analog (HMPD) have been measured over the 
temperature range 290 K -—- 450 K. The relaxation of the 
31p nucleus in HMPD has been assumed to be entirely due to 
spin-rotation interactions and the angular momentum 
correlation time (tz) has been determined. It has been 
found that the quadrupolar interaction in 144 in HMPA is 
modulated by a time dependence characterized by two 
correlation times: the overall molecular tumbling 
Correlation time (1, )@and that for the internal rotation 
about the P-N bond ee tg was determined using the 
Hubbard relation with ty obtained from the 31p gata in 


HMPD. EN has been determined from the 14, data. The 
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activation energy for the internal motion was found to be 
374.2 0-2) kcal/mole. 


The relaxation rate of 2H in HMPD was found to be 


BN 


= ano. that, Lor 


governed by three correlation times: tg, 


the rotation of CD3 groups Gr )s 


4 The rotation of CD3 


groups was found to possess a relatively high activation 
energy 3.0 + 0.1 kcal/mole. The delineation of the 
contributions to the relaxation rate of 3!P in HMPA was 
carried out using the above results and data from the 


temperature studies of the diffusion coefficient of HMPA. 


ee 


7 —F . . : 
a at, 
° eer ee Gonae’t one etree bison 


2 ¥; 
Gd  9aw nian ae te me Zz Sani yt “a 
mm. vents @ dialer sends thiviiore 
¢ i 
P21 eck agee, racageers omy. 
y 
a¢ Tae) Gole (Levis eter | ty. “ies | Gav oy 


en Ki wt :ciseht ine GiOeS Bet veer 


ah qa) 
6 
ee ne 


nic) ces pal 


ay 20m I 
gti P edt) nibs let og? 2 Ueienytereee 
ve 3 sibiesc 4) ie Ae heatara 


e123 an ity 14 @2)9eteeRee 


hte ry 


ACKNOWLEDGEMENTS 


I take this opportunity to express my deep gratitude 
tomProfessom RvE.D.. McClung sion his guidance, and 
enthusiasm during the course of this work. He has 
patiently spent countless hours leading me through the 
rough terrain of density matrix and spin relaxation 
theories. 

I owe special gratitude to Dr. T.T. Nakashima for 
many stimulating discussions and collaboration in many 
experiments. 

Loam also, gratetul to Proressor R.B. MJordan for his 
help with instrumentation and experimental aspects. 

I thank Dr. J. Ripmeester for measuring the 
phosphorus spectrum of solid hexamethylphosphoramide. 

Special thanks are due to Annabelle Wiseman for her 
expert help in typing this manuscript. 

To all my friends, I express my sincere appreciation, 
especially to Dipak Ghosh and Swapan Chowdhury for their 
continuous support. 

Financial assistance from the Department of Chemistry 
and Natural Science and Engineering Research Council of 


Canada is gratefully acknowledged. 


ox 


—- 


nhitiaece eee (al Gta © yt) oleae ee | 


ee @cnsthie ob: ©: (pee 1.20% aveaeetoel BF 


at a ee &r = L_ Cviei eged @rAene 


a Liwmae @ Boihes' eu Twine Ci enecreg: 


aasseeatees,.. caym Sie 4 elec f:e!t Gas 


«6 @e paaaes 

en oe: . - 
sam af 4 Pega (ts. Sei j é *5 shh. Qee 
vegas 

aie Wh cape .- 5 mS 
aeqbhe, 1)-+fret 7 r 5) Sanece Ho dep qa 

ad 64 re. 
ase ore FR? ot 4 2 : 41) @i2%peqeoeet 

nent wo) wameore? 37! amt ©) Ae 2 oh ' 7. La 
lef i $1ak fete A/M. SLAVE ee Cie 

sat eo reegyr esate » i} at [ a *) Gm cic gt 


+92 on \veuleorny.” a0 oe tose “ae } wait ae Vt hie : 
jp SQ2Rit. 0 Ces : 
eepelaet> tO saeco” tao Oe acne (ase Ip poten - 


tn ddan) MAMwERE un cogesh {UT bre goemEsy sah 
> er ise 
ai 


i 
: wl tenth 


- 


= 


: 


: 
a 


TABLE OF CONTENTS 


CHAPTER PAGE 


PART 1: DENSITY MATRIX SIMULATION OF TWO-DIMENSIONAL 


AND MUDTLPLE=PRPULSHaANMR EXPERIMENTS « 0006 o 6 61s 6 1 
dite APPLICATION OF (DENS DUTY MATRIX THEORY IN 
PULSED FT NMR smal ce eNisMeh aed eh edie: Sulers whe stile) sie 'd s.¢) rele sel eta ser'e Z 
ay oak Bohai erohee Wey gy conv; aaa Disa Gera Fee Vem ta 2 
Tz Definition of Density Matrix and 
Eguations for Physical Observables...... 3 
ee Equation of Motion of the Density 
Mat il Kirra Siewanct ret evere < srs eras Seas oh ceame: 3) SMen ay cL eh bee hse 5 
T.4 ROtating Erame @iransformat lone. «22 sae. D 
25 Inclusion of Relaxation Effects 
ange ChemLCa Wee xchange es sees ot ase oo ee 12 
6 Density Matrix Method and Bloch 
EQUA CIONS PAPP LOAN a acm so oh hey ele! wo ane oa has 14 
Dea, Applications StorPulsed rT NMR. & 2423016 sere 7 
58 Analysis of a Simple FT Experiment...... 22 
rao SPUTUM LV, stetacsteremets t8e. opelhs tees ores oiler eteristelehe ie. ois e084) 6 28 
eis DESCRIPTION OF THE SIMULATION PROGRAM.......2ee.. 30 
Tie E Introduction... 6. 2A, OD ober eects Soahon ers 5 o.o 30 
Dis 2 Objectives of the Simulation Program.... 31 


vd 


aren + Tezte 


oe te "ise 


(rs. ay 


Ol Wile, tie Teel ee Sour Ke) GPS TTS A 


> wes 


ae ee re 
: oo seu) mew shirae os idew ( ? 
as. tase”. | = : 7 
pes speeg/ pet) iat). OT - 
vous ; 
a - +e ya 3-36" 
: eg isine§ Of t Ofer? yt eow Th 
JVOe 6h» TA U 
east Ll a fete 
f= 2% Ge yf ot 9 oc (c® “Oawes § 
-ecae lS Lala ited 6@i-5 
ede . GLP eet =) Molina Ug? 3 
ee eee fa!) Bry) ow a 3 
er . ba hy: TT asi en 
sw oir 6 seas wma nym ere “5 
7 i, Nea: deer cassie 


asi aneaslg 


, we on 


7 a oe) 


a 


CHAPTER 


Tiss 


LL.6 


lass) 


16 ied HO) 


eek 


ee 2 


iS 


II.14 


Implementation of the Simulation 


Program. e e e e e e e ° . e e @ 10>)6)"'0' (6. 6 °O) @ (8. 0) 8.6) Oe Hot 2 Se Da PO Lee) 
dupe Desceriptson Of (Input Commands... ..4).. 
Choicey obs basic. Clk icc. «ess ceies See We 


Construction of Spin Operator 
Representation wMatind COS scumie.ie.s se os pss ove 
Computation of the Thermal Equilibrium 
Denstty Maeai xn ue eres Gels Weibslion oie c) eile kel el cneutertels 
Transformation of the Density Matrix by 
git Web MY-y UOlel wer yer obey ceadig SL Shae OI eee eee 
Transformation of the Density Matrix 
Puringerree. Precession Perlods...9..s.... 
Transformation of Density Matrix During 
Broad =BandeVeCCOUP LLG tec «steels one oa cette seis 
Computation and Simplification of Signal 
[ie qeb ereY crawl @ lest oy oC eee ae Cs SCTE OR Oe 
Identification of vy; and vg as Frequency 
OLrisetem frome Iransmattrerds sere ss ccs: 5 sicl sis tone 
Other Features of the Simulation 

PLOG EAM coc tls aces 6 Ageistems «kee oh Spode (oi aiielalee stains 
Limitations of and Suggested Extensions 
Comtiewolml Lat vOnPrOO rams si sslsls ssc s ieee 


SUMMALEY sie 21s. «s\6 6 silelleeiele ess 616 6 vis 6) 0's 160 4 6 e915 


x 


PAGE 


83 


BS) 


49 


50 


=) 


55 


61 


62 


63 


66 


Us 


76 


78 


VP 


“9 


ae 


* 


7 reset tha vile 


> : 7 : 
saitghinze +t! 40 dt a iniiel a : 
»>e¢ Se r= ee eli s . meyer" : 


ebyieenes we ' a6 ,<=000 ° 


ote 2 6 Ops 


; = so 74 fap ies = 
_ 
ie woe OR ever pd) 26 estan’ or ( 7 
; ; 
9? pimege eri ghee (i) oy gl! <2 DAL? af 
. 
? or te 7 ! 4 La 
grits s je ¢* set 4 ar 
poo tay! >. ie 
jaca De oi wie ih Sates . Sh. es 
: ee ee a 
pomiuiied? we Wi on WY G7 WE Pee Pa eat 
eT err eee ee, ee 7 
wlielkes= « ere geran aon tease ci.pe 
Py eer a rt eee Ey uieygllt 


gene st salt begeengve O18 we emp 
sank tel + ppeberty neg ¢) 1, wits 
et 


= 


CHAPTER PAGE 


III. THEORY AND IMPLEMENTATION OF THE HETERONUCLEAR 


NO TS Ee DE COUPTING «SCH M birs e sn sivas oars 1s. eee ectocdme: Chace 80 
Hf ot ARS IVC TOU Cte POM iereetersia te a cele eee ere ele clea: bao ne 80 
Pi. 2 THEOL Pi Ca teLE Ga GMON Es <.ecd.c eles ele Pies ie 81 


a. IS9 Spin VIS CT sist ere vierss ots: wlohe Wie) e i's: 6 ets 83 


Ope ES ana tS wepiiesyStems. ...- sn see ot 
¢. Spherical Randomization Decoupling.. hs 
DoS General Effects of Noise Decoupling on 


The MMenSsi eye Matis crs ad eels wwe cre, LOO 
III.4 Implementation of the Decoupling 

SCIVCMO Crees ots tetra seer sions cnet ieesi el sisters ee ele sere so OZ 
Sea A gs) SUMMA ey aster. A Get tote OD Oot Gre Oi 6 Gd oo om eee 


IDV ABP ESVGCATAON Stes chet coments! chele lower cteketen ster eters) siicleten creche erste c 106 


DVL pipelesateye sbWeh en Wey ete 5 Asin. 6 tert) 0 0.0 Cote OL Ota uray nD 106 
EV 22 Simulation of Experiments... ...- Peek ecetmeerne, «OT 
as “DEPT -Experiment 4a. 6.85.5. CpooomoounG, thew 


b. Selective 2D Heteronuclear Shift 
Corre lLaGrom sh Xperymentes 608 6 ore sec) she oes 130 
Cn “INADEQUATE EXperiment.4.50ce0.ne. oes «=§©6AO 
dad. Unexpected Effects of Spin 
DEGOUDDLNG is ce ctelse ee cleiee se sila este wees) | 49 


LV .3 Summa eye ot stebodere veers ele cole ee etalon rete ste ot haiieeehe 174 


X1li 


a 


> | 
— a 


nate meENT SF aie 7 worry Cebgiene AA ba dees “sat, 
era Ganesan ora 


ae pa? £397 


: 
a? ean ri == 


e& 


“e oe ‘se 8729 


aya) 86 Sie oti 


10 Acw )...emepe? 
wnt 6 Gh Qt a : 
<? ° ; 
2 
°@ — - rT A lA {sé 
ra. ° 
*@ Jah il eer arr w'apentiaaes i =k Fant Ze : 
u Ss 
see ow 1. Oe 59 Ad , 
a> * . an) ts if “) 40 
i ,, ory 
® a i) 
. © 


iangeti®) a JURY 


ov)? ale we ; 7 ee = y 
nesta )) Ot Se 
Dieters aT ad | 7 - 


sc al 
| oa ee tewndeeat 0 ahs eee ige c.¥1 
- 7 7 7 1 
oFPt sas wee + itagel yard Ya a 
Pia a 
sticge, ssi Guise sn Or os oe fed a 
kobe. ee 
ove ee Oo ees + dip gee cot beat 
- a> 2 a - 
we) °« aman Peg tases pe AS ; 
ae J 
nit ts «ET 
: me 4 wae e 
Or . 2 2eS'8 6 es 5 *= . es = . $7 L 
age 
a 
ee, 


ocr or —— ids veteG A ~ 7 


- * as 
a a a - — 7 
oe _ 


7 a 


CHAPTER 


PART 2: 


Vie 


NUCLEAR RELAXATION STUDY OF HEXAMETHYL- 


PHOS PHORAM IDE ete tote tote rore te tote fo tate 4m 6 te tei te Meme nie tetette, Yel 


DESCRIPTION OF, SPIN RELAXATION. MECHANISMS... .%.. 


Ven i ahengreyel ter 2 ah olele: 3 4 Sirol Ome nec Rone ee len 0 
Viz Spin Relaxation and Applications of 
Relaxation ts CUGiHES <r-is reer tas ehh ROSA OLD ene 
Veo Spin Relaxatwons Mechanisms. .4.444547444% 
an) DipOlarstintenract@ons. «64.0 5..00% ate resee Ns 
Te oRadCOmeROta tl Oliemes crore cs tte sistt se sisi 6 
20m ea Llansd ae Lona ee MO bd Oriya ist). is ts Saewes 
Di Mopin-hoOtatuGhe i nteract 1 Ons <-46 11s ers eis 
G. MNeetear Ouadrupolesinveraction. 15... . 
V.4 Brief Review of Rotational Diffusion 
MOG CUS. sence) ol cietietetelis to ome tee etcl oscars Rife Motay ei ccetettes wiles 
V.5 Ob JeCUrVeS One ther er OgeCts sce vives + eres 
V.6 SUMMamy Met tae ces es oie eee Petey ote te tel’ teenie tens 
BAPERIMENTA Di 6 stele sbec Siete ve "slo Sie) on 0 ts) tei i tele wile ie 6.10: SS 6 ae 
Wate ll iRejescalshblebe Wapato Toy cc tec eo Ler OnG Sy ciel anton 
Via SanpLeyPreparatious....s.4ss6vs sews a ss es 
Vile Instrumentation........ SN OG Sere Oe ateue 
Vine Measurement of T) Relaxation LimMecce sis) ors 


a. Inversion Recovery Experiment..... ay 


Dee tTeiplet SequenGe rss ccc dsiesvsect evans 


YL 


PAGE 


176 


Ly a/. 


RIT 


are), 


182 


184 


186 


188 


She) 


oz 


193 


L968 


9 


200 


200 


200 


201 


204 


204 


205 


- 


7 paaen ol 


wr, ‘ees 
i ,% 


_ 


- a 
eer d “n ale sutdenla gots | ¢ 


® er orn ; : _ ea gppacar Rea 4 


> 
ee ee ee ee 7 
eae ee j-% 


- , ’ 
ie egniieel lifyh Bag xe-” NO ee = 
are ou s ae 
- DO °¢ 
_ aa Pa er oT vor (Lee fine de 
~ D 
enter (oge(se. MaKe) : 
i 
_- 
@ 7 '« @ ry 
Tia “ae ae 1S >| 4 2 } a 
. : s 
Tomi 
e ® 
= re be (8S oH i Ca¥ 
o < Lom? (= * ad 
s © = @ e s 
wa O.,80TT @¢ Ft fates vee -t® 
» 8 eS ¢ Ran 4 Bie © cs 
4: Wes ) a ghaleadil “t rom I isi¥ 


o> Se As y ae lel? eee de ie 
Aer sista pene eae = Aeted aintegeant 


#2 yeoe tr pene 

ara. VIR 
eielbaatho 
— 


ae 


a 


CHAPTER 


ViI.6 


VII. RESULTS 


Vie 


Vallee 


Vib 


VII.4 


WERE GS, 


Vil. 6 


or Lies? 


REFERENCES.... 


Measurement of Translational Diffusion 


COCTRLCLEN tS wer. s.- 
SUMMATLY sce severe es ore » 
AND ®DISCUSSION 333. 


Pitroduceron es... 


31p T, (HMPD) Results and Analysis 


14, T, Results and Analysis....... 


2 


H Ty) Results and Analysis 


oeoe ee ee ee we we ew 


31p T, (HMPA) Results and Analysis...... 
a. ‘Calculation of -Spin-Rotation 
Reva XaC TONE RACES «seis ere sere os 


jor" Calculation of Intermolecular 


Dipolar Relaxation Rates..<.- 


eee ee 


eC. Calculation of the Intramolecular 


Dipole Relaxation Rates..... 


di Comparison of Calculated and 


Experimental Relaxation Rates. 


COnCHUSBOMNS sca sues ste 


e 


SUMMA Vises) s suslct ets! seususieleienels clielcersre 


kkkkk 


O60" (610) OF 6: 0 6 6) 6: 6 6 6116 04 0); 0) 26:0" 1617.0) 6 +0) (6) 10:6) (0) 6) 0: 6) 6 Oe. 


Sov. 


PAGE 


207 


PED 


22 


212 


ZZ 


221 


230 


236 


239 


241 


243 


244 


246 


249 


250 


eis 


as 


oon 


esiaudt ie’ teemitelaees |) s9tietwens- 


> Gs. OSS 8 OS): 


® suet dent) 


> : ; m ed a sr 


-~—= §eOGes4e?*.' 


ateves* 7 ,TT > 4 Oe CEs 4 
rane ye sae | 459 TTI 9 
sasetient ‘Se Site" ae f 
e., of ras rehen lilpoe AF OG 4h L yes 
0d 09 wa, 2 
elaei segs Faa © . as-7 BRrY «2S pa’ 
i¥* ec 
ee * © «@ 2 we Be a 
=p Dnt erve sas 19 £7000! er ls n 
‘ sence «Oe 766 j 
90 [cis (afas on: ’ 1. =i eo : 
so SOs wy? sada 6 PAG 


——=— Ge © & 


See CogeifateT no ite tape 


© 
v 


on (lenee ce Dretetes Cf a at! seen 


a ee as Je re @viatls Os 8 190" = 


yy eer ee ee ee TS eal is 


<.ce 


». Vv - 
op pied, - STV. 


i a 
£40" a 
-*@ ay 


i itv 


a 


Lisl OF TABLES 


TABLE DESCRIPTION PAGE 

1 Program input format for density 

Mate Kee Vee ie cee un eesiee cu abewele snus) ois Las oha pe cos ee rerers 49 
2 Relaxation data and results for 3!p in 

PIMP Decisis pemsuete tals cMetaler syeksuclspeusycusks us eguiindehsstel ei stecs) «ers 214 
3 Relaxation data and results for 14n in 

FIM ED, is sve ce wkeke us, sodeele fos couse le sks te teile sci ls,s (seus 664) 6 ole < 222 
4 Relaxation Data and Results for 7H 

INUICHLOLIS:. one ueie wou epee pehvaspsacde is boks oivitelobcds fells (6 is els) osirela!s 231 
5 Relaxation Data for 7!p in HMPA......... Ses | ASKS} 
6 Comparison of calculated and observed 31p 

relaxation rates Sli MPA i grers « syekele eA tr Ohta 241 
il Translational diffusion coefficients of 

HMPA gen EUy ah Oi Ones Oe Fer ee. ie Mg MOR On eae 


XV 


weed 


Pe ee ee jugs “eI HoOs 4 


anata ‘VLIN4 


s @ © &e 


FIGURE 


AL 


2 


10 


iE 


DIS OF) FIGURES 


DESCRIPTION 


Block diagram of the simulation program.... 


Pulse sequence for calibration of pulse 


WPL Hiteras acc etokce teens (o’ 6) ies 'e' 6 6 6) 6 6 oeeeeee ee ee © © © © eo 


Pulse sequence and input data for 


simulation of the CH DEPT experiment..... 


Output of the DEPT ‘simulation -for -CH 
SVS SI che epee cten Melons remensber otetenewest: s0e.6 oy rites of on 
Magnetization diagrams for CH DEPT 
EXPSELIMSTU ease csr ter steeds wn ates fire sie ls wi ene 


Input data and output from the 


simulation of CH coupled DEPT experiment... 


Simulation of §-pulse angle dependence 


PIMCH SER PI SxpeiwePMen Crags. ts as aborts 6 es s¢.6 6 6 


Density matrix and output from the 
Simulation of CH5 DEPT experiment...... 


Density matrix and output from the 


simulation ¥ory CH. DEPT experiment ..-..%... 


Magnetization diagrams for CHy DEPT 


EXPCFiMent.....eescves Ghee sn etal Nel ates 3 sks 0 le tes sis 


Magnetization diagrams for CH3 DEPT 


experiment...... oeoee ee oeoeee ee ee we ee oeoee ew ee @ 


XVi1 


PAGE 


108 


ie) 


iee2 


118 


BOS 


2 


ee 


123 


124 


eantars 2 Tex2 


mt Sala Der 


» Go (ekaeeIS Gale W. t 


i «PIrIOIG ve Ses omis oes 


- j yen seis ay : _ 


pater So Ani eeFit = 


Ga eo » @ » 4% ¢ 7 ee 
i f iHjeut- a ei ° 
- oe - g en ry 4.@ ¢ To 
ci sB--s ¢ mieeere a Ve mow ft . 7 
vp 7 PI SSu ice » aoe oe de a 4 
aE ' = , aee oF Ga q 1 ° ig 
" 522 2% née @ 
s/ a é ; o DTS tree 
, 
ewe ail perc: a. Roe Frm a 
Oy’ ie ye pet sane | cha twAr® 
taht r= & fs 6.1 OER 7 
- - 
Ore = TT éekisfaa’." Fe ue 
ia wet eyufah Vie eoatgr yee wed e 
sGi > ito. 2 J, MOR LAH *eS mL - Bf Gey OSSS fin? 2 = 


er peers peg lve Pes viaser > @opea 
oer yeagée srry. ivr joanne gt ta ans 2eWunce 
SOD) gi a : Soames pti raak tena | 


= 
ek this pSepesy- ese meee reas om Yan 
: : > - 


aa uae z 
Sp 145 Ry . i or _ 1 Geass eS 7 


_ a 


FIGURE 


TZ 


vs 


14 


15 


16 


1) 


18 


Tg 


20 


DESCRIPTION 

Output from the simulation of CHo 

coupled DEPT © experimentl «..% Niet oe een ee wo 
Simulation of §-pulse angle dependence 

PnP CH, DEPP  SxXperUinent ss iets ws + sos cee ele 
Output from ther simulation of CH, 

coupled = DEPT SExXperUMeMN cis «1516 6 e's 6 wie les es sens) os 
Simulation of 6-pulse angle dependence 

in CHs DEPT Vexperiment......... Se a eee , 
The selective 2D heteronuclear chemical 
shite "corretation pulse sequences. .% 1.1... 
Tnput data sandsoutput fromthe 

Simulation of CH selective shift 
COLCESTALCION "CXPeHLMeN tyes cis sxe 0 s1¢1s losis se sels os 
Input data and output from the 

simulation of CHy selective shift 
GOrrelat@won Experiment 25 26 eres vies eo ie Sch a 
tiput data®and = output “from the 

simulation of CH, selective shift 
GOLrrelatiion EXperUMent ss. Vie 2 Nie cs 6 cus le oi eis. 
Mputedatayand output £rom the 
4-phasecycle CH selective shift 


correlation experiment....-.eeseeee oes Goa 


KV 


PAGE 


126 


E30 


132 


134 


s5 


136 


oe 


oJ 


“= 7 _—_ . 7 ‘ 
> 7 


. 
: a 
; “puneget ie 
et? ST salgatiente ois er +2: wena 
if - out Rpenteatid opt, tins Capos eS : 
ovate sibs oe ti ow ‘Sedtel teh e 


( s=4R rae on Wa : 


» © 


: "9 naltwewn err inst euqiuy 


feed) za foxes Sie 
7 re = - a 
Liege mt Oe arlahh 1 aca Gans 422 en, 
ok Ue oeeep! Pegae? TSee Gian OL : 
rey ® 164 foounhzasbY “A Wr erseeees ; ag 
he) 23:00 TUAl. 
a be eee sea ‘sys | 69 Ge PO 7 
L pee!) al, 15 1e= 1S Th Sal rah eines : 


és 


. be ys 7 
cS) « OPPs oF syuwirl oi 


R ' 
7 (ete : =2\4 yas 28 wil te lewse 
nT die mea) AT ae aut 
is 


we i= *y Agee thas un eo Le 
Pe ee 
rea ‘see imyto lesa 

ead nv? or ae aay cola 


een a beaut ~) We 


eee Subetears Gist $4, 1g 
see-e i. - ites 4 * 


a 


FIGURE 


2 


22 


23 


24 


25 


26 


Pah 


28 


Yas: 


DESCRIPTION 

Pulse sequence and input data for 
Simulation of the INADEQUATE experiment 
on 13¢- 


Output data from simulation of 


INADEQUATE experiment on 13c-13¢ 


GE ACME HE « cylowemereine «eb meu ene ccsleetaMete tevaders (orig “a eho tone rie ly Lo 


Output data E£rom simulation of 


INADEQUATE experiment on 13c-l2c fragment.. 


Output data from, simulation, of 


INADEQUATE experiment on 13c-l3¢-1l3¢ 


ERAGMOIE io. s tacos ye too yore ees AEB Gal SEE Ce ee enn Ged ees 


2D gated decoupler 130 j-6 resolved pulse 


sequences........ ALAC OC OE DOS CRMC IO ree 


Inputs datayandwoutpur, Eroms simulation: of 


gated decoupler pulse sequences A and B 


olg) Otel oy seus eco Poe O Om De) he O ONCE RRO nO 


Output data from simulation of gated 


decoupler pulse sequence A on CHy system... 


Output data from simulation of gated 


decoupler pulse sequence B on CH» system... 


Outputs datas Eromesimulativon) of gated 


decoupler pulse sequences A and B on CH3 


SYS COM sisi cicve os eoevwewe#e LT ok wD ot Dat Tea gel Sok } OL18s O70): 8 Oo Oe) 621920 18! 48: 18 0: 


SOVAlALL 


136 CCM LUN ONS, oe ee eee RO et ea 


PAGE 


141 


143 


146 


148 


iS 


ie a)S) 


E56 


We 


peepee Steml fat siemens apie 


pean: — atin wt: te nowmetente 


Fispt om 
oi; w = @ ox 
omy ; : 7 papa © ae aOee 


Sa) neLse! ini: «=~ eh@® see? we 


lalate GrastoedAdt 


> ay i) 4 - 
‘ oe, ee 
oy © Ont ake waie® B88 ‘nae © ts: — 
x6 ay tthe? Ss af = | onde Een ieee 
7 hy te tue & i a8 >be ® s&s 
t ij o¢m FSA 
: ; scanQews 
ry . 
- 
be tavcousd tet f°! sy i gota Beles +3 
- 
, sO CU es - 
0 _ 
rygsiihs ate al et ws 
e 
wuths, bal | yee Pero 7 
: - 
é ‘ oo vs i fae es i ; 
ney, a uikertesd ea » Jugend ae 
am -¢..8e0dne® BO | f pes jerue atquargek. 


Weocl te? adhangs bo det = 
264 jt» Bahay ary ae wlygee ia anton welyioust: _ 
chet ot 
Se ee Y PN poe 
: : sn and: cot we 
me CEB regnedoe ate wie ; = eee 
rey 


FIGURE 


30 


od 


32 


SiG) 


34 


35 


36 


37 


DESCRIPTION 
2D J-& resolved spectrum of 2-butanol 


obtained with pulse sequence C (with 


decoupled aC eds DOM) mer encts + Sherer ese ghee enchoerers 


impute data andPoutput from simulation of 


gated decoupler pulse sequence C (with 


decoupled acquisition) on CH system....... 


Output data from simulation of gated 


decoupler pulse sequence C (with 


decoupled acquisition) on CHy system...... 


Output data from simulation of gated 


decoupler pulse sequence C (with 


decoupled acquisition) on CH3 system...... 


2D J-6 resolved spectrum of 2-butanol 


obtained with pulse sequence C (with 


GOUDMeEC ACG sr OM) sa) sister cusiersiersr enero eter sneer ote 


Output data from simulation of gated 


decoupler pulse sequence C (with coupled 


aeqwist t LoOmy sony CH SYS teM err acuel otersh sl siers cre! che 


Output data from samulation of gated 


decoupler pulse sequence C (with coupled 


acquisition) on CHy SYSTEM... .eeeeeeeeeees 


Output data from simulation of gated 


decoupler pulse sequence C (with coupled 


acquisition) on CH3 system.............06. 


xix 


PAGE 


163 


164 


TGS 


166 


170 


taf. 


isd 


174 


sanaioe ss e. Gs 7a ovthe og ane ae 
AY, Lae a denee sate 1G fede AAC 
(wo) oles jae bi i.qoe 78": 


an? a eae 
et <1 01) oe Ce ide “ig eg peace /€ 
ilgtie’ ’ 5 Te ae ( aseret pore 
ed j : or Gy fe 8 - nehbee . : 
— v 7 
ard Fed ey (P ' 7 : ab al 7 oe 
7 
nfist 1a fat} ee ee ee 
| | 7 
» so) > 2 nc Lie if at 4.ia aj ed - 
Le 5 soe! ye 
_ — 
S in} 1,203 ">? Pia 3 >? ear vy 
ip Fae a yo 6b SOHTESH Syés © ope 7 
7 iv _nie ee 
- Jette AFis sAi POKES 
ee ak eines : 
aim 3 297 ne — aie ranges 
' 
Sipe, gta ke) Se soni maine oast 
sal A 49-10 ont re oe Ne: 4 whaleespity 


Sesau, 1 co 84 Nim et? Ain sone 


nies Sie 


Dehigaey S208) * Os 


PIGURE 


38 


39 


40 


41 


42 


43 


DESCRIPTION 

Observed relaxation rates of 3!p in HMPD... 
Observed 3lp Spectrum OL SOlLa HMPD? 1... ; 
Observed relaxation rates of !4N in HMPA... 
Variation of 4,4 with +t, for HMPD........... 
Observed relaxation rates of *H in HMPD.... 
Comparison of calculated and observed 3lp 


HelaxvatLony Date sate MMP Arter. tel ciel olere lev sles siele ele e 


XxX 


PRESS 


27 


Ee 


Hab 7 


- 


- 


7 


aie «(awe uth ay wari pet ‘as ; « Peres 
OL) Sa PATE nse WE Samat 


eaten eee 
7 


° si yi Ta an 


alee 1 Poo 6 Pay 7 by enriatrey 


oe 


"The outcome of any serious research can only be to 


make two questions grow where only one grew before." 


Thorstein Veblen, 1919 


XxX1 


- 
> * 


Jka 2 ere pay om Or Ais AAT 
7 
: 1. wesbseagh bees eae ek 


oe ‘| ee 


ee Een > ae 


PAR ae. 


DENSITY MATRIX SIMULATION OF TWO-DIMENSIONAL AND 


MULTIPLE=PULSE NMR EXPERIMENTS 


CHAPTER I 


APPLICATION OF DENSITY MATRIX THEORY IN PULSED FT NMR 


ar | Introduction 


This chapter discusses the application of the density 
matrix method in analyzing pulsed FT NMR experiments. A 
brief review of the basic aspects of density matrix theory 
is given, starting with the definition of the density 
matrix, prescriptions for the derivation of physical 
observables from the density matrix, the equation of 
motion of the density matrix and its solution, followed by 
a discussion of the incorporation of relaxation and 
chemical exchange phenomena into the equation of motion. 
The differences between the density matrix method and the 
classical approach, Bloch equations, are pointed out. 
BPinally, the application of density matrix theory to 
pulsed FT NMR experiments is illustrated by following the 
evolution of the density matrix through a typical pulsed 


FT NMR experiment. 
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I.2 Definition of Density Matrix and Equations for 


Physical Observab les 


The time-dependent Schroedinger equation for a system 


described by a wavefunction Ww is given by 
fv = inde , [1-1] 


where H is the total Hamiltonian of the SV CECM mil) amelie 
wavefunction ¥Y can be represented in a complete 


orthonormal basis set {|m>} as 
W(t) =r C(t) im, [1-2] 
m 


where Cri) are time-dependent coefficients. 
The value of any observable property 0, associated 
with the corresponding quantum mechanical operator Ge the 


fe pieley lene teil ae) 
O(t) = <w(t)]O|w(t)> (1-3) 


= * ~ 
~ z c#(t) C(t) <mi|O) n> ee 
m,n 
Since. Olt) isegiven in terms of products of coefficients, 
it is convenient to define a matrix P whose elements are 


the products of coefficients arranged so as to give 
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Pi = Crea) c4(t) : [1-4] 


Equation [1-3] can then be written as 


II 
M 
ig, 


Oe) (t) <m]Q|n> [1-5] 


nm 


= 2 Pomlt) Qmn = 2(P(t)Q) ny = Trace(P(t)Q). 

m,n n 
Equation [1-5] describes the property Q(t) for a 
particular molecule, but experimental measurements give 
the average <Q(t)> per molecule in the ensemble. 


Averaging over the ensemble, we obtain 
SO tee racer (or) Oa. Dod 


where the density matrix p(t) is the ensemble average of 


P(t) and has elements 


nie = * 
p44 ft) Cc, (t) c5(t) ; 
where the bar denotes an ensemble average (1,2). Equation 
[1-6] gives the prescription for deriving average values 


of properties from the density matrix p(t). 
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The density matrix p is the quantum mechanical analog 
of the classical ensemble density. pp contains all of the 
time dependence of the system and defines the system 
completely. The density matrix po is the representation of 
the density operator , in a particular basis set, with 
elements 
aS <i o)j>". 

Since the representation of any operator depends on the 
particular basis set used, the density matrix has 
different forms when represented in different basis 
sets. s Although the form off the matrices, 6 and, @ depend 
on the particular basis set chosen, the trace of their 
product, i.e. the average value of the observable, is 


independent of the basis set. Hence the density operator 


can be represented in any convenient basis set. 


I.3 Equation of Motion of the Density Matrix 


The time-dependence of the density matrix can be 
ascertained from the time-dependent Schroedinger 
equation. The time derivative of the density matrix 


element pn, 25 given by 
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= ace) dc*(t) + dc (t) c* ft) ; [1-7] 
at dt 


From Equations [1-1] and [1-2], one obtains 


‘ diem 
H | m> = 1h pata | m> ’ [1-8] 


from which the time derivatives of the cocfficients Abe) 


can be obtained as 


tome : 
ab Oo vere ons air <n|H|m'> 
m 
and [1-9] 
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Substitution of Equation [1-9] into [1-7] gives 
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which can be rearranged to give 
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Ory un, Operator form, as” (1) 
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In the absence of any explicit time dependence in al 


the=solution to Equation [1-10] is given by 


o(t) = exp{-ift/nh} 6(0) exp{ift/n} , [1-11] 


where 6(0) is the density operator at time t = 0. 

Provided the Hamiltonian H is time-independent, the 
density matrix at any further instant in time can be 
computed from the matrix representation of Equation 


PTs) 


T.45) Rotating Frame Transformation 


When the Hamiltonian contains a periodically time- 
dependent term, the equation of motion for 0 is solved by 
transforming the density operator into a suitable rotating 
frame so that the effective Hamiltonian in the rotating 
frame is time-independent. When the Hamiltonian Het was 


expressed as 
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the equation of motion for the density operator is given 


by 


Ales 3 [orHy + Ay(t)] . iii] 


OQ. 


In a frame rotating at frequency w about the 
laboratory z-axis, the density operator o*(t) is related 


to the laboratory frame “density operator by 
Peers exp G—Wutet ep Gtimexn iat Ulli [1-13] 
where it is the z-component of the total spin angular 


momentum of the system. The equation of motion for o* is 


obtained from Equation [1-12] by recognizing that 


Gat ele oat ot +255 era ie 
cane ee Oe exp( 1515) ae exp (iol t).. [1-14] 
One obtains 

in Wav ares Tol H,*(t) [1-15] 
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where 

Ay (t) = exp(-iwI,t) A,(t) exp(iwI,t) PLS Liol 
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is the time-dependent interaction in the rotating frame. 
In all cases of interest in magnetic resonance, Ay (t) has 
a time-independent part which causes important changes in 
o*, and a time-dependent part which is far off-resonance 
and has no influence on the magnetizations of the 


system. One therefore takes the equation of motion for p* 


eo, De 

do* A Aa 

ae = = lo*, Hegel - [1-17] 
where Hote = Ho + Hol, + time-independent part of Bee 


Since Hore is time-independent, the solution to Equation 


[1-17] is given by 
p*(t) = exp(-iHgget/h) p*(0) exp(iteget/n) [1-18] 
and the density operator in the laboratory frame will be 


i tom=nexD ole inn p*(i)mexp(-iel stn. fil=19 | 


tt is* instructive to illustrate some of >the” features 
of the rotating frame transformation by considering some 


simple nuclear magnetic resonance experiments on a single 
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Spin 1/2 particle. The following examples also serve to 

illustrate the similarities between this quantum 

mechanical rotating frame and the rotating frames 

imtroduced in the solution of the Bloch equations. First, 

the case of a time-independent Hamiltonian is 

considered. The spin 1/2 particle is subjected to a 

constant magnetic field of strength Hp, along the z- 

axis.) ine Hamiltonian for this system can be written as 

Hy = -yhyt, [1-20] 
= -hugl, ’ 

where wa is the Larmor frequency (in radians sec1), ys 

the magnetogyric ratio of the nucleus, and es is the z- 

component of the nuclear spin angular momentum operator. 


The Hamiltonian operator of the system in the frame 


rotating at angular frequency w, is given by 


Hee eed ul tel 28) 
= Al(wo-w)I, . 
When the system is exactly "on resonance", wo = w and Here 


=p silisc leags to the result fEquation [i—17]) 
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which implies that o* is time-independent (in the rotating 
frame). The Bloch equations give exactly the same result 
in terms of magnetization vectors in the classical 
rotating frame (3,4). 

The case where the Hamiltonian has a periodic time- 
dependence will now be considered. The system is 
subjected to an additional oscillating radio frequency 
field in a direction perpendicular to the constant field 


Hos the Hamaltonian is given by 
H = —huol, = Vas COS suet i A [1-23] 


where H, is the amplitude of the time-dependent magnetic 
field oscillating at a frequency w. Transformation into 


the rotating frame yields (Equation [1-16]) 


Ay = - —y— - —z (1, cos Qut - 1, sin 2ut). [1-24] 
The first term in Equation [1-24] is time-independent and 

is the "on resonance” Component of the rf field. The last 
term represents a field oscillating at frequency 2w in the 


opposite direction. This counter-rotating component is 
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far from resonance and has virtually no effect on the 
mMagnetizations. The last term can therefore be 
neglected. The effective: Hamiltonian in the rotating 


frame is then written as 


Hoff = -h(wo-w) I, = 5) iL ’ [1-25] 


where Hore is time-independent. For the exact "on- 
resonance" condition, it reduces to 

i -yhH, re 

Meppee So ty [1-26] 


Equation [1-26] shows that, at resonance, the effective 
field in the rotating frame is just the H, field along the 


X-axis. Bloch equations also give the same result. 


Peo ine lusion Of (Relaxation Bitects “and “Chemical Exchange 


The effect of spin relaxation caused by a stationary 
random perturbation rip eae has been incorporated into the 
equation Of motion of the density matrix by Redfield (5). 
Using perturbation theory and making use of the ensemb le 
average properties of H,(t), he has shown that the 
equation of motion of the density matrix due to relaxation 


is given by 
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relaxation 


where Roa’: 8B! is the relaxation matrix element given by 


1 


Raa’: BB! FS Re Ta co ce Jupaigt (2-8) 7 
oe ; ay a owe : Jey igy toh a ’ [1-28] 


where a-B is the frequency corresponding to the transition 
states a and 8, and Toatppt io) are the spectral densities 


CL); 


Jyatgg(o) = f° <alHy(t)[a'><p* [Hy (t+7)| p> e°2°*dc[1_-29) 
where the bar indicates an ensemble average. 

The equation of motion of the density operator has 
been further modified to account for the chemical exchange 
processes (6). A chemical exchange process has the effect 
of permuting the nuclear spin states, and the average rate 


of change of the density operator can be given by 
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where an exchange event occurs on the average every 1 
seconds and P is the permutation operator. This approach 
to chemical exchange is elegant for simple product basis 
sets, but tends to be complicated when product bases are 


unsuitable (strongly coupled systems). 


26 Density Matrix Method and Bloch Equations Approach 


The Bloch equations provide the classical description 
of the motion of non-interacting spins. These equations 
(1,3,4), are the equations of motion for the magnetization 


components M,, M,, and M,, in the presence of relaxation 


M 

effects: 
ee, ee eee 
dt Mi Zz an i 
aM M 

23 vee at WR ae AS 1-31 
dt Cs H), i ’ [ i 
yy "Y 
at a y (M x He aa we ’ 


where Mp is the thermal equilibrium magnetization and Tj 


and Ty are the longitudinal and transverse relaxation 


times respectively. 
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The Bloch equations can be derived from density 
matrix theory in some situations by evaluating the 
expressions for the expectation values of the 
magnetization components. Both methods, the density 
matrix theory and Bloch equations approach, yield 
identical results, in the description of many experiments, 
for the observable components of magnetization in weakly 
coupled spin systems. However, the applicability of the 
Bloch equations is limited. Some of the major differences 
between these methods are amplified in the following: 

(i) Exact expressions for the relaxation terms can be 
derived using the density matrix method, as was done in 
the Redfield relaxation equations. In contrast, the 
relaxation terms included in the Bloch equations are 
purely empirical. The inclusion of these terms is 
justifiable only because they seem to explain experimental 
observations. Density matrix theory provides, in 
addition, the conditions under which the Bloch equations 
are applicable. 

(ii) The classical description of NMR experiments, given 
by the Bloch equations, are not applicable for strongly 
coupled spins. Here it is necessary to introduce quantum 
mechanical energy levels and operators rather than 
classical magnetization vectors. For weakly coupled spin 


systems, the simple products of the eigenfunctions of the 
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separated spins are eigenfunctions of the time-independent 
Hamiltonian Ho for the complete spin system, and hence the 
transitions between various energy levels are identifiable 
as corresponding to the spin flips of individual nuclei. 
In such cases it is possible to write down Bloch equations 
for the magnetization components of each individual 
nucleus. However, in strongly coupled spin systems, the 
Simple product functions are not eigenfunctions of the 
Hamiltonian Hp and the transitions cannot be identified as 
apleing serom the spin’ flip of “any particular nuclei, 
Density matrix theory must be used to describe the time 
dependence of such strongly coupled spin systems. 

(iii) The Bloch equations are concerned only with 
magnetization components (one quantum coherences), not 
with zero quantum and multiple quantum coherences, and 
hence fail to describe the creation and evolution of these 
invisible (directly unobservable) coherences in multiple- 
pulse experiments. These coherences can be produced ina 
spin system initially at thermal equilibrium by 
application of two radiofrequency pulses separated by a 
delay (7). Multiple quantum coherence phenomena have been 
utilized in multiple pulse and two dimensional FT 
experiments to provide additional information about the 
spin system. In such situations, the density matrix 
theory must be used in order to describe the signals and 


the quantum mechanical state of the system. 
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They Applications to Pulsed FT NMR 


The density matrix theory can be used to describe 
pulsed FT and CW NMR experiments. Since the theory is 
capable of explaining the origin and behaviour of multiple 
quantum coherences, it is especially useful in the 
analysis of many multiple pulse FT experiments. However, 
the explicit computation of the density matrix can become 
tedious depending on the size of the spin system and the 
pulse sequence employed. 

A typical pulsed FT NMR experiment consists of the 
application of a set of pulses to the spin system, 
separated by periods of free precession followed by an 
acquisition period where the transverse components of the 
magnetizations are sampled. The pulses may act on all 
nuclei (homonuclear systems) or may be selective and act 
only on one group of nuclei (heteronuclear systems). The 
flip angles and phases of the pulses, duration of the free 
precession periods, spin decoupling during precession 
and/or acquisition periods, etc., can be used to select 
the information about the spin system which is obtained in 
the measured response. The transverse components of the 
magnetization are measured as a time varying signal 


referred to as the free induction decay (FID). In order 
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to improve the signal-to-noise ratio, the signals from a 
number of scans are co-added to obtain the final FID. 
During the individual scans, the phases of the pulses may 
be varied to suppress undesirable or spurious signals or 
to correct for pulse-length imperfections, pulse-phase 
abnormalities, and for differentiating magnetization 
arising from various orders of coherence, i.e. zero 
quantum, single quantum, and multiple quantum 
coherences. The final FID is Fourier transformed to 
obtain the conventional frequency domain NMR spectrum. 
Recently, two dimensional FT experiments have been 
developed (8). Here the signal is collected as a function 
of two time variables. It is then Fourier transformed in 
both time domains to yield a frequency spectrum which is a 
function of two independent frequency variables (9-12). 
The advantage of two dimensional experiments over the 
conventional FT experiments is the simplification of 
comp lex spectra which can be achieved by selectively 
presenting a particular type of spectral information (e.g. 
chemical shift) along one frequency dimension while 
Simultaneously displaying another spectral parameter (e.g. 
spin-spin coupling) along the second dimension (13-15). 
Thus the total information content of the conventional 
spectrum is partitioned along two dimensions, which 


greatly simplifies the interpretation of spectra, 
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particularly those of large molecules. Another 
application of two dimensional experiments is in 
correlating the chemical shifts of coupled nuclei in 
homonuclear and heteronuclear systems (16). The 
correlation of chemical shifts is effected through some 
mechanism for magnetization transfer between the nuclei 
involved (17). The different coupling mechanisms utilized 
for magnetization transfer are scalar J coupling 
interaction (18-22), dipole-dipole interactions (23-25), 
and chemical exchange processes (26,27). Depending on the 
shift correlation information desired, a two dimensional 
experiment can be performed utilizing any of the above 
mechanisms for magnetization transfer. Pulse sequences 
can also be designed to utilize the multiple quantum 
coherences in spectrum editing (7,28-30). It is possible, 
for example, to edit the spectrum so as to differentiate 
signal components arising from various spin subsystems. 
This technique has been used to distinguish between 
signals arising from CH, CH», and CH3 spin systems (eam) s 
Pulse sequences for all types of NMR experiments can 
be represented as a series of radiofrequency pulses and 
free precession periods applied to a spin system which is 
initially at thermal equilibrium. In applying the density 
Mato xetheony to the aescription Of athiese experiments, one 


begins with the thermal equilibrium density matrix. It is 
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convenient to choose a basis whose members are the 
eigenfunctions of Hos the Hamiltonian operator in the 
absence of radiofrequency pulses. The Hamiltonian 
Operator corresponding to a multiple pulse experiment 
would be different at each stage in the pulse sequence. 
Hence wml telse preferable tomthink iof the entire pulse 
sequence sas’ a sequence of «steps, “such «that during any 
particular step, the Hamiltonian operator has a definite 
time dependence. The state of the spin system during the 
experiment can then be followed by monitoring the 
evolution of the density matrix in each successive step. 
In a multiple-pulse sequence, there are basically two 
types of Hamiltonian operators acting on the spin system 
in a sequential manner: radiofrequency pulses and free 
precession periods. Solving the equation of motion of the 
density matrix during the free precession is 
straightforward since the Hamiltonian is time-independent 
(Equation [1-11]). The time dependence of the Hamiltonian 
dumungethe tappliicataon: Of saneck pu lsencan bemremoved by 
transforming the system into a rotating frame with an 
angular frequency equal to the frequency of the rf pulse 
irradiation and the solution of the equation of motion for 
the density matrix (Equation [1-17]) is again 


straightforward (Equation [1-19]). 
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Thea Cimes evolution, of therdensityamatrix in a 
multiple pulse sequence can be represented as follows: 
Let p(0) be the thermal equilibrium density matrix 
operator and let ae Sipy. Seances Hn be the sequence of 
Hamiltonian operators representing the pulse sequence. 
The operators are time-ordered with the earliest one 
leftmost and successive ones to the right. Here it is 
assumed that the time-dependence of the operators has been 
removed by appropriate rotating frame transformation. The 
final density operator at the end of the pulse sequence is 


given by 


ee eer Ls (Vo (0) Dea) ere) hies24 
where ny is the propagator during the ith time interval. 
Using this) prescription at is then, possible sto follow the 
evolution of the density operator under any given pulse 
sequence. However, at the computation stage, it is 
preferable to work with the matrix representations of the 
operators rather than the operators themselves since it is 
easier to apply the matrix transformations rather than 


evaluate operator transformation equations. 
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I.8 Analysis of a Simple FT Experiment 


The application of density matrix theory in the 
description of pulsed FT experiments is nicely illustrated 
by considering the pulse sequence used for routine 
measurement of Spectra. The sequence consists of a 90° 
pulse followed by signal acquisition. The detailed 
Svouutaon o£ , during this sequence will be considered for 
an Gnsembie Or (systems containing a single spin 1/2 
nucleus. 

The system is subjected to a static magnetic field of 
strength Hj, along the z axis. The Hamiltonian for the 
system is given by Equation [1-20] with eigenfunctions |a> 
and |p>, corresponding to eigenvalues -fiw,/2 and +fu,/2 
respectively. From statistical thermodynamics, the 
thermal equilibrium density operator 0 (0) can be written 


as 

p(0) = (1/z) exp(-Hp/kT) , [1-33] 
where zis the partition function of the system. The 
density matrix will be diagonal at thermal equilibrium. 


In the high temperature limit, 


exp(-Hj/kT) = 1 - Ho/kT , [1-34] 
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and the density matrix elements are 


GOR a = CA 2:5 p(0) gp = (1-A)/2 and p(O) gg = p(0) py = 0, 


where A = yhHo/2kT. 
The magnetization per molecule for the ensemble at 
thermal equilibrium can now be determined from Equation 


[1-6] and has components 


See tsace (0M; ) ele eK Vp [1-35] 
where 

em 

and Tie ie-eX oy, 2 ace wivenm Dy. the Pauli Gspin 


matrices. Evaluation of the magnetization components 
gives 


<M,> = 0, <M,> = 0, and <MZ> = (1/2) yha . (ls Gn 


As expected, the x and y components of the magnetization 
vanish since the off-diagonal elements of the density 


Mathiesaresatilazero. Note that <M,> 1S proportional to 
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(PagtPga) and <M, > 1S PDLOpOEeronal to (Dag - PBq) so that 
transverse magnetizations are present only if p has non- 
zero off-diagonal elements. The value of SMO is 


proportional to Pag) and, in this case, is simp.ly 


(Pag - 
the thermal equilibrium magnetization. 

Inv tne = presence of the rf pulse, the Hamiltonian, for 
the system is given by Equation [1-23]. As discussed in 
Section 1.4, the time dependence of the Hamiltonian H can 
be removed by transforming the system into a frame 
rotating at frequency w about the field Ho, such that 
YH) >> (w-w5)- The equation of motion of the density 
Operator in the rotating frame, then becomes 
-iyHy 


die fori esa 


Wlell SOMULLON, 


o*(t) = exp/( 5 VCO 2 Vie [1-38] 


ed Nora tiny cy is the angle through which the magnetization is 
tipped by the pulse and will be denoted by 6 (in 

radians). The representation of the exponential operator 
can be obtained from its series expansion, from which we 


have 
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exp(16é Ix) = 5 [1-39] 


POL cater “piu lsenoL lip angileso0os.6 = 1/2y.eand the 


representation matrix in Equation [1-39] reduces to 


= ares . [t=40)j 


For Simplicity the density matrix is viewed in the 
rotating frame in the following discussion. After the 
transformation indicated in Equation [1-38], the density 


Masi ithe trotating: frame t1segiven by 


o*(t) = 3 Te . [1-41] 


Evaluating the components of the magnetization per 


molecule in the rotating frame, we obtain 
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aM) Mal OG <M (a) >o = Se, land <M, (2)> = 0% [1-42] 


The 90° pulse along the x-axis of the rotating frame has 
rotated the magnetization from the z-axis to the y-axis of 
the rotating frame in keeping with simple magnetization 
vector diagrams. 

The spin system is now allowed to precess freely. 
The effective Hamiltonian is given by Equation [1-21]. 
The representation of exp(iHg¢¢ ea) is diagonal in the 


Simple product basis, given by 


exp{i(w-w,)t,/2} 0) 
exp(4Hore ee) = ’ 
0 expt si (G>a5) &y/ 2) 


[i=43] 


where ty ise Ghesacquiisiti10n, time ivamrapite.,) Ine density 


matrix during the free precession, therefore becomes 
0 —iAexp{i(w-we)t,} 


LOCXDP(= TCs a) yy 0 


[1-44] 
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During the free precession period, the non-zero off- 
diagonal elements develop time dependence. 

The NMR Signals in the rotating frame, can now be 
computed using the expressions for expectation values of 


magnetization components to give 


OM Girrtaess = CA/ 2) San toma. yee oF 
<My(ty+t)> = -(A/2)cos{(w-wg)ty} » [45] 


0. 


These calculations suggest that the magnetization lies in 
the transverse plane (xy plane) throughout the detection 
period t,. If the relaxation effects were included, we 


would have obtained 


<My(ty+t)> = -(A/2)sin{(w-w9) ty} exp{-ty/T2} , [1-46] 


<My(ty+t)> =(A/2)co0s{(w—w,) t4} exp tay To. . 

Spin-spin relaxation damps the oscillatory functions, with 
the characteristic rate constant (1/T>)- Equation [1-46] 
gives the time domain magnetization functions for this 
simple experiment. The conventional mode of presentation 
of the NMR spectrum is in the frequency domain. The 
frequency domain spectrum can be obtained by Fourier 


cosine transformation of SMe ie or <My(t+t) )> time 
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domain functions. Fourier transformation of the time 
domain functions in Equation [1-46] gives a Lorentzian 
line “centered at frequency VV Cin Na) WLC ert ewrat it 
at half-height of 17 Cpt) (1n8Hhz). ~ Te showld™ be noted 
that w = 2x v, in Equations [1-45] and [1-46]. 

The example of the spin 1/2 system illustrates the 
way in which the density matrix theory is used to describe 
pulsed FT experiments. “An analysis using *therc lassica! 
Bloch equations would have yielded identical results in 
this case. For larger systems with strong coupling or for 
sequences which involve the generation of multiple quantum 
coherences, the Bloch equation approach is inadequate. 
Even with larger spin systems, the basic guidelines in the 
application of density matrix theory remain the same, 


however the computational complexity can increase rapidly. 


I.9 Summary 


In this chapter a brief description of density matrix 
EHeOry and sitseapplications in) the: description of (pulsed 
FT NMR experiments was presented. The equation of motion 
of the density matrix was set up and the application of 
the rotating frame transformation in the simplification of 
the equation of motion described. Modifications of the 


density matrix theory to take into account the relaxation 
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and chemical exchange processes, were briefly discussed. 
Fina Vly, the application of the density matrix theory in 
the description of pulsed FT experiments was illustrated 
by considering a simple example of an ensemble of spin 1/2 
systems subjected to the most basic FT experiment, the 

90° pulse-acquire sequence. The next chapter describes 
the implementation of a computer program to carry out 
density matrix computations on more general spin systems 
under a variety of multiple pulse and two dimensional FT 


experiments. 
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CHAPTER II 
DESCRIPTION OF THE SIMULATION PROGRAM 


leet Tncroduct lon 


The implementation of a simulation program to analyze 
multiple pulse and two dimensional FT NMR experiments, 
based on density matrix theory, is discussed in this 
CGuapter, Therobj ectives "of the program are’ Jaid out, the 
SELUCTULe SOL thesprogran sis 1ilustrated with) a bp lock 
diagram and brief descriptions of the major subprograms 
are given. The data structure for the representation of 
the density matrix is described. The computation of the 
thermal equilibrium density matrix is illustrated with 
examples of the IS and IS, spin systems. The 
implementation of the density matrix transformation by a 
rf pulse and during free precession periods is discussed, 
with the transformations illustrated for the IS and IS, 
systems. 

The effects of broadband decoupling are considered 
briefly. The equation of motion for the density matrix in 


the presence of a broadband decoupler irradiation, and its 
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solution are given. The simulation of the signal 
acquisition process is considered next. The algorithms 
used to compute and simplify the signal expressions are 
briefly described. The interpretation of the frequency 
terms in the computed signal is clarified by carrying out 
the analysis of a simple pulse sequence. 

The implementation of various other features in the 
program, such as phase cycling of the rf pulses and the 
receiver, rf pulses of any given tip angle, are 
described. Finally, some of the limitations of the 
current version of the program are discussed along with 


suggestions for future extensions. 


iio peOoReCti Ves: Ob tnegSimulation Program 


The program has been intended to simulate the effects 
of multiple pulse sequences on a variety of spin 
systems. The development of the program has been 
motivated by the conspicuous absence of a device for 
analyzing the rapidly increasing number of two-dimensional 
NMR (2D NMR) and other multiple pulse experiments. The 
simulation method is based on the density matrix 
description of pulsed NMR discussed in the last chapter. 
Some work on the computer simulation of multiple pulse 


experiments has been reported using Bloch equations (32) 
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and density matrix theory (33-35). Recently another 
approach to simulation of multiple pulse and 2D NMR 
experiments has been devised, making use of the super-spin 
formalism (36-38). Simple rules for following the 
evolution of the density operator under rf pulses and 
during free precession periods have also been given (39- 
41), but the operator expressions tend to be prohibitively 
complex in a typical multiple pulse sequence. The 
computational complexity associated with the density 
matrix method grows rapidly with the size of the spin 
system. The need for automated density matrix calculation 
is therefore apparent, and the aim of this thesis has been 
to develop a rather general system to accommodate a large 
variety of pulse sequences and to simulate many of the 
features used in modern high resolution NMR spectroscopy. 
The spin systems that can be analyzed by the program 
are homonuclear and heteronuclear I,S, systems of weakly 
coupled, spin 1/2 nuclei with 1.< m,n < 3.. The input 
format is similar to the pulse microprograms used with 
high resolution NMR spectrometers. The specification of 
the pulse sequence constitutes: the definition of the 
spin system; specification of the pulse angles; pulse 
phases and nuclei on which the pulses act; delay periods 
and signal acquisition periods with or without broadband 
decoupling; options to carry out phase cycling of pulses 


and detectors; and so on. 


seen! pl ieeeee 2 Es 


| 7 ¢ Sle ytd 


age : ova nwt eal ne 
3 Oo - 
ysis Sele isriere 20teluiees lLeoqoterug = 
_ an 
- a : “wee aioe 
’ c ve ey t~— = 8 : oo i - » _ = ©) are 7s ‘ 
as ® 304 137 i c3 ay) later oe eager ee * ere? - 
ui SAT Gainwie: oy 
rc oy poss @4 47 i “zee 
(4 it, ied mn ® pee eo te i. 
' ong i wae (4? ae leg@ Clan ant 


a 
ad. oun’ 4 a? a 1 sau dewronnt ot 


1 A, aie aL | Gat: SL aaae 
=f ’ a a +9 je meals oe _ sat he al nI2 
bacatilies aad td 


- + af: Tiet oe te 8A _ ach es ap pe 


; Sad 4.) wre Taal =e ; 
salen dares 410g att Ceeereetls 
. ‘ 7 

ofard) 4). i Ye ~ 359 - L.¢ 


a ~ > ~ 
figsed Pn: 3s de as hd [ 


Pulses can be, applied tito either dor nuclei “Was in 
ayhemeronuc lear system)s, orm to: both nuclei (as sin a 
homonuc lear system), and along the +x, -x, +y, -y axes of 
the rotating frame. The program is also capable of 
Simulating the effects of pulses of any given flip angle, 
besides the most common 90° and 180° pulses. Since our 
early experience showed that examination of the elements 
of density matrix at different stages in a pulse sequence 
enabled us to investigate the behaviour of magnetizations 
and multiple-quantum coherences, a facility to output the 
density matrix was incorporated into the simulation 
routine. 

The computed signal can be simplified to some extent 
so aS to make the output more comprehensible. Another 
important simplification feature has been incorporated to 
allow explicit substitution for the durations of 
precession periods in terms of 1/(nJ), where n is an 
imveger)andaj. is, thejscalar’ coupling constant. between ithe 


tT .and S nuclei. 


II.3 Implementation of the Simulation Program 


In order to have a program with reasonable speed, 
efficiency, and portability, the FORTRAN language was 


chosen rather than an algebraic language like REDUCE2 (42) 
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for the implementation of the simulation scheme. The 
proyram is organized in a modular fashion, with a short 
main program and a number of subprograms, each with a 
distinct function. This feature simplifies the task of 
expanding the program, since modules to perform new 
functions can be added on with only relatively minor 
changes to’ the short main program. The main program calls 
on different subprograms to handle various steps specified 
in the input pulse sequence. Subprograms pass control 
back as soon as their function is performed and the main 
program chooses another subprogram depending on the next 
input specification. Communications between the main 
program and subprograms are made mostly through common 
blocks of memory. 

The density matrix elements are the most important 
and complicated items to be represented in the simulation 
program. The elements of the density matrix are, in 
general, complex and depend explicitly on the durations of 
precession and decoupling periods. In order to strike a 
balance between the complexities of pulse sequences which 
could be handled and the storage requirements, the number 
Or timervablab.wJes. 26 limited to, fourm: . Ti, T2, DL, and 
Does Bach ih these .ariab les: can, be jwsed. to.define the 
duration of a period of free precession or a decoupling 


period, and each variable may be used no more than twice 
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within a given sequence. Each density matrix element, at 
any particular point in the pulse sequence, will consist 


Of a-Sum Of terms of the” form: 


Oe eS) exp 2m vie 12 ey Dl ye D2) 1 orale 


where A + iB is a complex number and Vir Vor Var V4 denote 
frequencies which describe the evolution during time 
Pericdseil, 12, Dl; D2 respective ly." *Ihe frequencies Vue 
vor V3r Vg are equal to frequencies for transitions 
between particular spin states so that the specification 
of each frequency variable requires the specification of a 
pair of indices (encoded into a single integer) which 
identify the states involved in the corresponding 
transitions. Each term can therefore be represented as a 
complex number and four encoded integers which specify the 
frequency variables in each of the time domains. A 
collection of these terms, each with distinct exponential 
factors, completely defines a particular element of the 
density matrix. In order to economize on computer memory 
requirements, the density matrix is represented as a 
collection of elements, each consisting of a number of 
terms, in single dimension arrays with appropriate arrays 
of index variables to specify the density matrix element 


to which the particular terms belong. The numbers of 
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terms associated with each of the density matrix elements 
are also stored in a single dimension array. This data 
structure for the representation of the density matrix 
Wadi iepee rec erred tOtase datasaruay... FOR convenience min 
accessing and manipulating the density matrix, the index 
of the starting row and the total number of the rows of 
tne data array" “occupiediby the density matrix, ‘at a 
particular stage in the simulation are stored in another 
array Tevrerred toeas the "Status array". Since sthe 
density matrix is Hermitian, only the elements in the 
upper triangle are stored. The data structure 
representing the density matrix resides in a common block 
and 1s thereby accessible to the subprograms which 
manipulate and update its elements. 

A block diagram of the program is given in Figure l. 
The functions of the main program are to read and echo the 
input pulse sequence, and to control the flow of execution 
by selecting the appropriate subprogram to handle the 


particular step in the input. The main program also 


Carcies, out sthne task, of Inltlainzing «the values Of scertain 


variables, constants and arrays. The subprogram SPNSYS 
computes the initial thermal equilibrium density matrix 
appropriate for the spin system specified in the input 
data, PULSES transforms the density matrix by application 


of the pulse specified in the input, and the 
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Figure l. 


PULSES MEASUR 


SPNSYS 


SUBSTR DECOPL 


OUTPTO 


TLEDYUP 


Block diagram of the simulation program. 
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transformation of the density matrix during the free 
precession periods is carried out by the subprogram 
PRECES. The magnetization components are computed by 
MEASUR and simplified to a certain extent by the 
subprogram SMPFLY, which recognizes the relationships 
between exp(io) + exp(-io) and cos » and sin 6. The 
parities of the sine and cosine functions are also invoked 
my EnIS Simplification procedure. After simplification, 
the final output may have considerably many fewer terms 
than are present in the internal data structure. The 
printing of the signal expressions is carried out by the 
subprogram OUTPT2. The subprogram SUBSTR allows the 
Simplification of the density matrix at intermediate 
stages by explicit substitution for the duration of one of 
the vcime. perrodas ol, I2, Dl, Of D2 by 1/(nJ). where 1s 
the I-S spin-spin coupling constant andn is an integer. 
Such substitution often causes certain terms of the 
density matrix elements to vanish and produces significant 
SimpJreication. “The transtormation of the Gensity matrix 
Guring broadband decoupling is carrred oul by DECOPL, 
which computes the density matrix during the decoupling 
period using an averaging process discussed in detail in 
Chapter III. The phase cycling of pulses and the receiver 
is handled by the subprogram PHYCYCL, which to a certain 


extent duplicates some of the functions of the main 
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program by controlling the flow of execution inside the 
piase cycling loop. As indicated *by the plock dlagram, “1t 
interacts with most of the other subprograms. It also 
keeps track of the signals generated during” each of the 
individual phase cycles and passes on the accumulated 
signal for simplification. Besides the major subprograms 
mentioned above, there are a number of other subprograms 
witch cancy Out various other functions. “One of the most 
important among these is TIDYUP, which organizes the rows 
of the “data array", and carries out some simplification 
by collapsing mutually cancellable terms, simplifying 
terms which can be summed together, and deleting blank 


rows formed during the simplification procedure. 


a. Description of Input Commands 

A description of the input commands to the simulation 
program is given below. Most of the commands consist of 
twoestrings. (Gach containing Foun characters), typed 
consecutively (without blanks or any other characters 
between them) on the same line. There must be only one 
such command in a single line of the input file. Command 
Nanecewhichballo short of four characters are tosbe typed 
in left justified fashion. The input is read from FORTRAN 
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Command 


SPINImSn 


HOMO 


PULSabcc 


PULSabccabe' Cc... 


Description 
Sets up the thermal equilibrium 
density matrix for the Tem Pe 


system with 1 < m,n < 3. 


Defines the spin system to be 
homonuclear. HOMO command 
should be specified before the 


SPIN command. 


Applies a pulse transformation 

to the density matrix. a 

Specifies the pulsed nucleus and 

can be Tor S. “b specifies the 

pulse length and can be any of 

the jfo Plowing: 

A : denotes a 90° pulse 

Bs denotesva 180° pulse 

C-F: can be used to assign any 
desired pulse length (see 
SET command below). 

cc specifies the direction in 

the rotating frame along which 

the pulse is applied and can be 

ANVSOLsE tN Xk, GV7o onUey. 


The command (3b) shows how the 
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PRECCG 


MEASax Yb 
MEASac 
MEAS ax YbaxXyYb* 


MEASacac"’... 


4] 


PULS command is to be used in a 
phase cycle sequence. c'c' 
Specifies the phase of the pulse 
in the second pass through the 
phase cycle (sec PHASECYCLE 


command below). 


Applies the free precession 
transformation to the density 
matrix during the time period 
tt, where tt can be any of Tl, 
12; Dl On D2.  qlt is preterable 
COmuse LUO LL? LO si Giart 
acquisition periods since the 
program carries out more 
extensive simplification on 
terms, Containing tl Om eTo. sseach 
of the time periods can be used 
no more than twice in a pulse 


sequence. 


Simulates the signal acquisition 
by measuring the magnetization 
components. Command (5a) meas- 
ures the complex quadrature 


magnetization. a specifies the 
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nucleus (I or S). b specifies 
the particular quadrature 
combination as given below: 


De=0 2 xX + 2y 


b= 1 : -Y + ix 
b= 2 : -X - iY 
b = 3 >: Y - 1X 


The output is onto FORTRAN unit 
6. Command (5b) measures a 
Single magnetization 

component. a denotes the 
nucleus (I or S) and c denotes 
the magnetization component (xX, 
Yor 2) SSc-and Sbeiliustratce 
how the commands 5a and 5d are 
to be used inside a phase 
cycling Joop. Here b: and Cc! 
specifies the component to be 
measured in the second pass 
through the phase cycle (see 


PHASECYCLE command below). 


Applies the transformation 
appropriate for spin decoupling 


during a precession period. a 
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SUBStt=1/(nnJ) 


PHASECYCLEnn 


specifies the nucleus to be 
decoupled (I or S), bbb denotes 
3 blanks, and tt specifies the 
time period during which 

decoup ling=1s "carrred out (LL, 


D2 aD, 2 Or 2))s. 


Allows the substitution of the 
Specified time variable by a 
term proportional to the 
reciprocal of ‘the coupling 
constant J. tt specifies the 
Cine Varad able s(T le ..l2> 1DiljroOr 
D2) and enn) 1s va twoedigit 


Integer Maine i272. Lorma t,). 


Allows phase cycling to be 
carried out. This command is 
used to set up a loop, 
terminated by another such phase 
cycle command. The loop 
encloses that part of the pulse 
sequence in which phase cycling 
is carried out. nn specifies 
the number of cycles and is such 
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format). The pulse phases and 
receiver phases to be used for 
the different cycles are to be 
Speciiied sequentially. inwal! 
PULS and MEAS commands. If the 
number of cycles in the loop is 
larger than 16, then 
Speciiveations fon 16 cyeles 
should appear on each line, the 
next 16 on the next line leaving 
the first four characters blank, 


and so on. The last command 
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inside the phase cycle loop must 


be a MEAS and there can only be 
one MEAS command in the loop. 
There can only be one phase 
cycle, loop) in ‘the puise 


sequence. 


Defines any desired pulse angle 
for thespulsessC—F 3 ep speciiies 
one: of pthe spulses  (C, Db, EB, OG 
F) and angle represents the flip 
angle (in degrees) to be 


assigned to P. The angle is to 
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Desspeciivedsindrse.!] Ttormac. 
bbb denotes 3 blank characters. 
b. SETspecaa Allows the user to set states 

(ON/OFF) of certain 

specification switches. aa is 

the switch setting and is either 

ON or OFF. spec is one of the 

following: 

(a) MOVE<s > 1 1€) 1S desired) ‘to 
store the density matrices 
abtereeachrstep, for print 
out or debugging purposes, 
the OFF setting can be used 
to prevent overwriting ot 
previous density matrices 
within the data table. The 
default setting is ON. 

(b) DBUG: Intermediate output 
for debugging purposes can 
be produced with the ON 
setting. ~The default 


setting is OFF. 


Oe WRDM Can be used to print out the 


density matrix after the 
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previous command in the pulse 
sequence. = Thehoutput sis onto 


FORTRAN@UNLt 68 


bes DUMP Gansbe used to print out. the 
internal representation (data 
Structure) Of the density matrix 
after the previous command in 
the pulse sequence. The output 


is onto the FORTRAN unit 6. 


2. COMMSstring Echo prints string as a comment 


Onto FORTRAN unit 6. 


53% STOP Stop calculation and terminate 


execution. 


14. NEXT Terminates the calculation on 
the current pulse sequence and 
Starts calculation on the pulse 
sequence specified, starting 


from the next line of input. 


15s STDM Stores the internal 
representation of the density 
mMatruix Onto, FORTRAN Ante 7s eit 
can be read in later in another 


rin to continue ‘the calcu lation. 
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16. 


RDFL 


NEDM 


(subcommand ) 


Reads sin. an) initial density 
matrix from FORTRAN unit 4. 

This feature can be used to 
carry,oOut, calcu lations. starting 
with non equilibrium density 
matrices. The FORTRAN unit 4 
should contain the SPIN command 
(HOMO and SET commands if 
applicable), the» NEDM and CLFL 
commands as well as the density 
matrix elements. 

This subcommand follows the SPIN 
command mentioned above, and 
initiates reading in the density 
Matrix. The density matrix 
element is represented as 
described im Section 1.3. Only 
the non zero and upper 
triangular and diagonal elements 
should be specified. There 
should be one line for each 
"TERM" (Equation [2-1]) of every 
density matrix element. Density 
matrix elements are to be 


defined one after the other. 
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Each line should contain the 
data described in Table 1. 
After all matrix elements have 
been specified, the input mode 
is terminated by inserting a 
zero in the field corresponding 
COMNE GMP [CTab lem l))r. 

CLFL (subcommand) Instructs the program to 
terminate input from FORTRAN 
unit 4 and to continue to read 
input from the normal input 


device (FORTRAN unit 5). 


Ll.45 Choice Of Basis Sets 


The basis set for representation of p is the 
eigenbasis of Ho which, in the weak coupling limit, is 
taken to be the direct product of the irreducible basis of 
the 1,, subgroup and the irreducible basis of the Ss, 

Sib group.) the construction Of ithe 1nreducibie basis in 
the simulation program uses a building up process which 
Starts with the e€igenfunctions of the single spin 1/2 
system, forms uncoupled simple product basis, then 
transforms the simple product functions into the 


irreducible basis functions using appropriate 


48 


> — 7 a 
on a a ra 
- a > 
, AP Alerne ot eo piveibal 
~¥ meite®, (eh nea} Te ei 
view «icnate haben Sabot | 


aes git) SRS perivqunge wee 


es 2nryrewnt if imetAmirree@: of 


- 
a 
' 
f 4 = ! 7 -_ re 
i? ; a. ) e¢ » ; 


(ene (Mt errnenr) bene vckes 5 MIs 6 


(Urge tT te leheid wee oe Save 


4 J : asl wee. hare 
e ror" 

~ w 
row ™ 5) eel Ser? 


ta ,408° to ene fel 


we Al Tre |, capt Suen’ 7a aye -agyes, wee. 
) rom? funwig Spy apc AEA ie ft: ataniel olay 
sae pl foe rid 
oe ats Ae ioe tent avai Spine! 
nd eA Ob pabars)) wet. “a bast cece - 
ie, (ere, WH ims jue deed 
(ol slug shpate gas + 


= —— 


49 


Table 1 


Program Input Format for Density Matrix Elements 


Variable 


NCOMP 


LDM 


COREE 


TAG] 


TAG2 


TAG3 


TAG4 


JCOEFF 


Format Description 
16 Specifies the number of "TERMS" in 


the particular density matrix element 


I6 The row and column numbers of the 
matrix element in the form: 
(i x 100) + j, where i is the row 


number and j is the column number 


2H14e 7 A complex number specifying the 
magnitude of the density matrix 


element 


16 The Tl time dependence encoded as 
(i « 100) +24, representing! the 
SxXOUO Tig Vel ee ta Chores) mc). pe b tol 


ij 
16 Corresponding variable for T2 
16 Corresponding variable for Dl 
16 Corresponding variable for D2 
16 Encodes information regarding the 


spin coupling constant J during 
various precession periods. Should 
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transformation matrices. The eigenbasis of Ho is the set 


a 


of eigenfunctions |1,M,;>|S,;Mg> of the operators 17, I,, 
S2, and Su with eigenvalues I(I+1), Mz, S(S+1) and Mc 
respectively. Here T is the total Spin angular momentum 
Operator for the Tn Sp rns sand S is the total Spin angular 
momentum operator for the S, spins. Since the basis 
functions are never explicitly generated by the program, 
Dut the representations of spin angular momentum operators 


are, the building up process will be illustrated in 


connection with the representation of the spin operators. 


II.5 Construction of Spin Operator Representation Matrices 


The generation of representation matrices for spin 
angular momentum operators is accomplished by starting 
with the matrix representation of the required operator 
for a single spin 1/2 nucleus, forming the representation 
of this operator in an uncoupled product basis and 
transforming this product representation into the 
representation in the irreducible basis by appropriate 
transformation matrices. This process is illustrated by 


aA 


considering the representation of ie andy Soin I,S}, and 


11S) spin systems. 


The matrix representation z of the z-component of the 


~ 


spin angular momentum operator for a single spin 1/2 
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z= (1/2) [2-2] 


inert representation of iL for the 1,8) spin system is 


constructed by computing the direct product of z with a 


~ 


Pee aU matrix U5 corresponding to the S spin, giving 


- [2-3] 


0 

= = = 

Tz = 2@) Uy = (1/2) 0 
0 


Sey @Q) l= 
a) || as) =) 


The basis functions [I,M;>|S,Mg> are in the order (M,,Mg) 
SG 2 2 172) 1/2) en ty 27 1/2) ly / 2) ee he 
Feopeccemearion OLS > 1s constructed by erevercing) the order 


~Nv 


of direct products in Equation [2-3] and 1s given by 


So = Ur ®) 2 ° [2-4] 


A more interesting case, involves the representation 


Of Ss 


~ 


a for the 1,89 system. The representation Zoe of the 
z-component of the spin angular momentum for the S5 Spin 


subgroup in the uncoupled basis is generated by 


Zo = Z@)Uy + Un @z ; [2-5] 
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where the first term is the representation of aoe and the 
second term is the representation of Se The matrix 
representation z> must be transformed to obtain the 


Gepresentation in the irreducible coupled basis of the S-, 


spin subgroup. This transformation requires the matrix 
il fe) fe) 0 

aaa ie 1/v2. Oo ye 

cad bie! 1/v2 @) -1//2 ; 
@) 0) al 6) 


which maps the uncoupled basis onto the coupled basis. 
The representation z4 of the z-component of the spin 
angular momentum operator for the S5 subgroup, in the 


coupled basis, is given by 
Zo Dz 5D ‘ [2-6] 


The direct product of the 2 x 2 unit matrix U5 with z5 


Gives) the representation S, for the full 1)S85 basis: 
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Here the basis functions |I,My;>|S,Mo> are arranged in the 
Order ey Mole= (1/21) ad / 2) 12,0), (= 1/2,0)), 
Ci2,—1)y <-1/27-1), (1/2,0), .(-1/2,0) with the first (6 
funGlvons associated with the triplet states of the S5 
subsystem, and the last 2 with the singlet state. 
Similarly the representation of De for “the 1155 "spin 
system can be constructed by taking the direct product of 


BZ witheU,, a 4 4% 4 unit matrix “corresponding to S55 


~ 


subgroup: 
Tz =2@)U% - [3-8] 
The matrix representations of other operators oe eye oe 


and ce are likewise constructed, building up from their 
corresponding representations for the single spin 1/2 
system. The matrices D for the transformation of spin 


Operators from uncoupled to coupled basis for 2 and 3° spin 


1/2 subsystems are stored in the program. 


© Computation of the Thermal Equilibrium Density Matrix 


Tievcirse, step in the Simulation, Of any pulse Sequence 
is the specification of the spin system I,S, and the 


construction of the thermal equilibrium density matrix. 
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In the high temperature limit, the thermal equilibrium 


density operator given in Equation [1-33] reduces to 
e(0) = (1 - Ho/kT)/n , [2-9] 


where 1 is the identity operator and N is the number of 
spin states (22*™ ¢ h 

P or the I,S, spin system). In the case 
of weak coupling between the I, and S, spin groups, the 
eigenbasis of Hg, ‘constructed in Section 112.3,) is the set 


A 


of eigenfunctions {|I,My>]S,Mg>} of the operators [2, 3 


S? and Be with eigenvalues I(I+1), M;, S(S+1), and Mc, 
respectively. Hence p(0) is diagonal in this basis and 


has elements 


p44(0) = [1 - (yzhvyMz3 + yghvgMg3)/kTI/N , [2-10] 
where the contributions from the much smaller spin-spin 
interaction terms have been neglected, and Mr 5 and Ms 5 are 
the eigenvalues of is and Se for’ the j-th basis 

function. Since the term 1/N Ls. common! to all vorpethe 
elements of 9(0); it is unaffected by rf pulses and 
decoupling fields, and it gives no contributions to 
measurable properties of the system, it can therefore be 


ignored. We take 
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p55 (0) = ~AIMr45 — AgMgs , Pole ey 


where A; = (yyhv;/NkT) and Ag = (yghvg/NkT) are 
constants. Using the above considerations, it can be 
shown that the thermal equilibrium density matrix for the 


I,S} system can be written as 


AqztAs fe) O 0 
oO) = (1/2) pee 212 
O 6) 0 pATeic 
In the simulation program, Ay = Ag = 100 for homonuclear 
spin systems and, for heteronuclear spin systems, Na Se 


and Ag = LOO 


Bi ip hranstormatton Of “the Density Matrix Dy eiiebulse 
Operators 


During the presence of an rf pulse, the Hamiltonian 
operator is time dependent. This explicit time dependence 
of the Hamiltonian operator in the laboratory frame, is 
removed by a transformation into a coordinate system 
rotating at an appropriate frequency as shown in Section 


iespeand tne solution to the equation of motion is) given 
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byPEquatven (1-181. “Anvanalysis of the effects on an 1,8] 
spin system, on the application of a rf pulse of definite 
duration to the S nucleus, is illustrated below. 

The Hamiltonian for the Tren sein system, in a frame 


rotating at frequency v about the z-axis, is 


De a(t very) Sed ic 


(yay / 2m) “eye anjson. [2-13] 


The last two terms arise from the rf field which has 
frequency v and strength 2H,. The spin-spin coupling term 
can be ignored since its magnitude is very small compared 
to the Zeeman interaction terms. When | vc-v| is small 
compared to Y¥gH,/2n, the rf field is on-resonance with the 
Suspins ana as) far from resonance for the 2 spins so thet 
the frequency offset term for the S spins and the 

ayy zt) le term can be ignored. The effective 


Hamiltonian then reduces to 


Since the first term in Equation [2-14] describes the free 
precession of any I spin transverse magnetization, which 


is not really of interest here and it can be treated 
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separately, we will ignore this term during the 
application of the rf pulse. The density matrix after the 
application of the pulse is therefore related to the 
initial density matrix by Equation [1-18]. Hence the 
transformation operator, Ry (6), which describes the rf 


pulse is 
R, (6) = exp(iygH,tS;) = exp(idS,) : [2-15] 


where 6 is the angle through which S magnetization is 
flipped by the pulse. In order to carry out the 
transformation in Equation [1-18] the matrix 
representation of the transformation operator exp (ies, ) piel 
the 1,S, basis set must be computed. 

The transformation matrix is constructed by a 
building up process which is exactly analogous to the 
process described in Section I1.5 for the construction of 
the spin operator representations. The transformation 
operator Ry (6) for a pulse along the k-axis of the 


Botatanguceame,, for.single spin 1/2) nucleus; is) (4743) 
R, (6) etexp (ies yl cos (0/2) at 2a Sesin(e/2)e [2-16] 


The direction index k can take on values +x and ty 


corresponding to the four phases commonly available on 
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spectrometers. In general, the transformation matrix for 
theeparticular I, or Ss, “subgroup “is obtained by taking 
direct products of the single spin representations 
(Equation [2-16]) to produce a representation in the 
uncoupled simple product basis, then transforming this 
Simple product representation to the irreducible basis 
using an appropriate transformation matrix. The matrix 
representation of the pulse for the full Inon system is 
obtained by taking the direct product of the 
transformation matrix for the I, subgroup with a Pee 
UnYC Matrix -Gpulse gon. lenucilei*on ly); Vthe \airect. *product 
Qf gae2 i 2 init tmatrix with (the transformation amat Giux 
for thews, subgroup (pulse on S nucle1 only), or the 
direct product of ithe transtormation matrix for the i 
subgroup with the transformation matrix for the S, 
Subgroup. (pulses. ‘to both, P and S nuclder in a nomonucJear 
system). 

The case of an 1,S, system shall be cons@dered..)  iie 
pilsesmarnix hCG) elon aart pulse applied along the x-axis 
of the rotating frame for a Singlei=spin 3/2 nucleus, in 


nem? 1/72 ply 27-2) Basis nis 


Byat cy) = ’ [2-17] 
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where C = cos(06/2) and S = sin(@/2). By analogy with 
Equation [2-4], the direct, product of Uy with R,(6) yields 


the pulse transformation matrix Roy (6) (meyer galas: | See hiesh 199 


Systems: 

c O is O 

O C O 1S 
Rax(0) = U>@)Rx.(8) =] : [2-18] 
a * a LS O C O 

O es O e 


HOmetene Ge pulse applied to the -1 yspins, «the ,corresponding 
transformation matrix Ryx(6) is constructed by reversing 


the ouder of (direct products: 
Rrx (0) = Ry(8) @) U2 - [2-19] 


If the rf pulse affects both nuclei (a homonuclear spin 


system), the corresponding transform-tion matrix is 


c2 ics ics -s2 

ics c2 =S- ics 
Ree era a) Re Oe fo= 20) 
ee OES: 1GS= ae? Ce ics 


Po2> S4cs) ics ic 


since both spins experience the effects of the pulse. 
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Pulse representation matrices for larger spin systems 
are also built up in a similar stepwise manner in exact 
analogy to the construction of spin operators described in 
Section 11.6. 

The subprogram PULSES, which handles the rf pulse 
transformations, analyzes the input command to identify 
the nuclei to which the rf pulse is to be applied, the 
dupavrvonrol ther rf pulse, “and the direction im the 
rotating frame along which the rf pulse is to be applied 
(referred to as the phase of the rf pulse). The 
appropriate pulse transformation matrix is computed and 
the transformation indicated in Equation [1-18] is carried 
out, by retreiving the elements of the previous density 
Matrix, one by one and carrying out the matrix 
Multa piicarvons OF Equation [1=18)|.0 “Arter each jelement ot 
the final density matrix has been computed, a preliminary 
Simpieercation and organization ol the terms of =the 
density matrix element are carried out. At the end of the 
transformation procedure, bookkeeping information 
concerning the final density matrix, such as the start 
index and Size in the "data array”, is updated in the 
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LLeSs Transformation of the Density Matrix During Free 


Precession Periods 


During a free precession period, the spin Hamiltonian 
is time independent, and the density matrix in the 
moOtating frame evolves according to Equation Pl-—l))},. The 
spin Hamiltonian during a free precession period is given 


as 


= -vzI, = v8 ap (ufdl S ° [2—214) 


a itaefes 


Computation of the matrix representation of the 
exponential operator exp{-iH(t5-t,)/nh} is relatively 
straightforward, since the {|I,My;>]S,Mg>} basis set is an 
eigenbasis of H. The representation matrix is diagonal, 


and has elements 
<j|exp[-i(t,-t,)/m]] i> = exp[-2ni Ej(tp-t))] , [2-22] 


where E5 we, the energy (in Wz) ofethe jen state. 9le can 


be shown that the effect of the transformation in Equation 


[1-11], is to change the density matrix elements as 


Pog (t2) = Ppg'ty) exp[2niv,g(t2-t))] , [2-23] 
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where a a Eg = Ep: is the frequency of the transition 
betweennthe states p andrq;: 

As pointed out in Section II.3, the time-dependence 
of the elements of the density matrix is represented by an 
encoded pair of indices which specify the states whose 
energy separation corresponds to the frequency of 
OSGi pationsduring 4 particular time interval. (During a 
fueehprecession period sof duratione?l§ (wheresTl = t5-t) ), 
for example, the terms in Ppq will all be multiplied by 


exp(2ni v Tl) so that the encoded integer array 


Pq 
associated with the Tl time domain must be updated by 


incorporating the encoded integer pair p,q into each term 


in Pg: 


ILe9eatransformation of Density Matrix During Broad=Band 
Decoupling 


It will be shown in Chapter III that the behaviour of 
the density matrix, when a noise decoupling field is 
applied at the resonance frequency for the S nuclei, is 


described by 


Non re ee pistel 


[2-24] 
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where the Po'g'spaq are numerical coefficients given in 
Seerion Pit J3¥and es is the “decoupled frequency" 
corresponding to Vog (1.6. the frequency of the transition 
from state p to state q in the absence of scalar spin-spin 
coupling between the nuclei I and S). 

In the implementation of the decoupling subprogram, 


DECOPL, the coefficients P are grouped into a 


p'q';pq 
number of small matrices and stored in the subprogram. 

The) manipulataons of tthe density matrix datarstructure for 
a decoupling period are straightforward - one computes the 
appropriate linear combinations, then updates the encoded 
integer array associated with the decoupling time variable 
in the same way as for a free precession period. In 
aqdition, a tiilaqd is Set “to indicate to the outpul stages 
of the program that the spin-spin coupling constant J is 


topberset to zero in this time period, so that decoupled 


frequencies will be computed at that stage. 


Mad0 | Computation wand Simp lieication of (Signal Pxpressions 


The signal measured by the receiver in NMR 
experiments is proportional to the transverse components 
of ithe magnetization for a particular) type of nuclei.» The 
expression for the physical observable 0, is given by the 


Equation [1-6]. In NMR experiments, the physical 
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observables Q are either the individual transverse 
components, for example in the case of I spins “ou 
<I>, or, in quadrature detection mode, the complex 
components formed from these transverse components. The 
various complex magnetizations that can be formed from 


these transverse components are <I,+il,>, <I,+iI_,?, 
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The simulation of the receiver in NMR experiments is 


accomplished by constructing a matrix representation of 


the appropriate spin angular momentum operator La le 

Le hil ey ieee (spin I magnetizations 
f , 4-3 a +4 ‘ tel 

y y Sy pica tes iS_y S_y 15.9 spins 


magnetizations only), or TS, I. +S., (tees ENR Ee 


-x! 
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ey 


Only!) G.0G S 
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x 
Ute >) (both I and S magnetizations for homonuclear 
systems). The construction of the matrix representation 
of a spin angular momentum operator is discussed in 
Section II.5. In homonuclear systems, the sum of the 
representations for the spin operators for the I and § 
Sones) constructed. “The calculation (Of ane 
magnetization is effected using Equation [1-6], with 0 
replaced by the appropriate spin angular momentum as 
listed above. For completeness, the simulation program 
allows one to "measure" the z-component of the 


magnetizations although these cannot be measured directly 


with a spectrometer. 
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The subprogram MEASUR, which carries out the signal 
computation, first analyses the input command to determine 
the magnetization component which is to be extracted for 
Lhesspeci fied Spins. Ttethen iconstructs, thesmatrix 
representations of the appropriate spin angular momentum 
operators. In the computation of the magnetization 
components, since only the trace of the product of p and 
the appropriate spin operator matrix is needed, the 
diagonal elements of the product matrix are computed and 
summed together. Simplification of these terms is carried 
Out by recognizing the relationships between exp(io) + 
exp(-id) and cos » and sin @, where » refers to the 
encoded integer arrays associated with the time variables 
Pigeslay Di ang. D2.) The subprogram  SMPFrLY scans “througa 
the rows of terms in the "data array" associated with each 
element of the density matrix and picks out pairs of rows 
which can be simplified and can be expressed in terms of 
sine or cosine functions. The information regarding the 
ey pe sOLmtne stuncCtdoOn w(1.e, 8 sine. ,— COSING OG exp) sis 
encoded in appropriate integer arrays for communication to 
Eheaprinting subprograms. | Furthermore, the four “integer 
arrays associated with the four time variables are decoded 
in the output stage and the frequencies Vg associated 
with each time domain are given as explicit linear 


combinations of Vie Vs and J so that the time dependence 


65 


jectece bd Jun. Pele ear aolite 
piiaroiadoiM sapere Sige | ong 
«it Epieadaew hee SY Mere see 
oi nian @) abumerae nae Fi amnlge S01) cage ad 
wire eligns Agr sap Ciee ees te enol deinsesagex’ 7 
ine) tee éeie atts $300) @iegees, pf . etoseange 
ive . te Pa, Bilt Packlee eit) Yee Ovair ,eseeogney 
aes ee Ab ALS tine uke sae! ohms COATES ans 
ap -etuseos a MM wane ioe ot a bay ers o, Lanegai® 
Tia =e mit Sot eigid al? Tyee . peetepgor Geom 
joe ees. are riet Sey Tne sR Qadri qgerss yates 
~j ewan o \peanwiys (ote Pie eo emg tne iel~sqxe 
-i¢ am Bete eepelieece Syeven segetn) Setoses 


pore enou: vost gia) Cyt Seg ‘Tei Ght 
sy thop Detarmen Mae a th at PEjainet dn aver ond 
na ts welt con oun’ hab Mena adhimel Seite Teeaale 
viet nh paSrrGs “4 nn ane DY ani ed @ac ache ; os 
qurte Secmy Veneers 1 ar eae one easaeo-sp 00 
, | Ate "er Serbabe” J reales att ‘pel hom au hs se 
J) IPI oS Lees “1019 Syeear 7 a \ a - 
cayeyins qu). wall \ otras itp slate poe 
beedice ond GRR MAINED meee 


Ne Sas S0FO ay 
names} 


can be ascribed clearly to the chemical shift offsets 
and/or spin-spin coupling. This explicit dependence on 
Vin Vs and J is also given when the density matrix is 
Ise eOUCDUL. It Shou ldspe, pointed out. chat vz and vo in 
the output represent offsets from the corresponding 
Eransmitter, frequencies. —Theajustitacation for this 
interpretation of frequency terms is discussed in the next 


section. 


Piel identieication, Oly, and) ycpas meequencys,Orrsets 


from Transmitter 


It was indicated above that the frequencies vz and vo 
in the output represent offsets from the corresponding 
transmitter frequencies. In this section, a detailed 
analysis of a simple example will be presented to 
illustrate that despite the numerous transformations 
between various frames of reference required in the 
density matrix treatment, the signal output by the 
Simudation program is justified if “the? frequency terms in 
the output are considered as offsets from corresponding 
transmitter frequencies. 

Tt) is important to’ recognize that the density matrix 
and the Hamiltonian matrix which appear in Equation [1-17] 


must be in the same frame of reference. The 
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straightforward procedure for relating the density matrix 
in the laboratory frame before and after a pulse requires 


the transformations 


Pulse 
~(lab) n et) Transformation  «*(rot) " ~( lab) 
Sa 
ce) (t,) P (=e) Cee) 


[2-25] 


As an illustration, the evolution of density matrix will 
be followed in detail, by carrying out the sequence of 
transformations in Equation [2-25] for each pulse in the 
sequence shown in Figure 2 applied to the 1,8, spin 
system. This pulse sequence is employed for the 
determination of pulse lengths on nuclei, when they are 
not under direct observation (44,45). 

The elements of the density matrix in laboratory and 


rotating frames at time t are related by 


Gror) 


(lab) 
p = (8 
Pg 


(cy) (t) exp[-2ni(M -M )vt] , [2-26] 
Pp q 


where Mp and Mg are the eigenvalues of the z-component of 
the total angular momentum operator cs in the p and q 
States. At t = 0, the laboratory and rotating frames 


coincide and the initial density matrix is 
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fe) fe) 0 hee 


where A and 6 are the parameters defined above. A 90° 
pulse is applied, at time zero at the I resonance 


frequency vy and the density matrix in the frame rotating 


ab y= Vie after the pulse is 
A 16 ) 0 
(rot) -1i6 A @) O 
Pi ) <= ’ [2-28] 
= *p 0 @) -A ié 
6) O -1i6 IN 


where T, is the time required to effect the 90° pulse. 


Transforming to the laboratory frame, we obtain 


A a O O 

avs A 0) 0) 
as = ' [2-29] 
ee 0 © eh & 

0) @) a™ SAN 
where a = i6 exp[2ni v} Td and the asterisk denotes 


complex conjugation. According to Figure 2, the system is 
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allowed tojtprecess freely for a period t,; during which the 
density matrix evolves under the time independent 


Hamiltonian and becomes 


A b 6) 6) 
cet) b¥ A O 0 
5 (+ +t) = 0 0 : ' [2-30] 
O @) ore ih 


where b = i6 exp[2ni(vy Th + V12 eg Ny 

e— 16 expl2ni(y;, Th + v34 cond Vi5 is the frequency 
of the transition between states i and j. At the end of 
the precession period, a 90° pulse is applied at the §S 
resonance frequency. To handle this, the density matrix 
must be transformed into the frame rotating at the 
frequency Ve Of -thesS transmitter to facilitate the 
transformation by the S pulse. The density matrix in the 


S rotating frame before the pulse is given by 


A d 6) @) 
ax A O 0) 
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where d = i6 exp(2nilvic. + vy9t1 - vg(tptty) I}, 
e = ié exp{2nilvtt, TP Vaacy es vg (tptt,) I]. 

Por sthe sake Of ‘simplicity, only thes density mater 
elements of interest p,5 and p34 which are associated with 
transverse I magnetization will be considered further. 


After the application of the 90°, S pulse, the elements 


x 


P12 and p34 in the S rotating frame are given by 


ORT (tty tas) =p ae eee a = 


i 34 a: 


0 0 
i6 exp{2ni(v_t_tyv (ttt, )} {exp (2niv 


ONS t,) + exp(2niv 


34tz)}- 
po=30] 


de 


Transformation of these density matrix elements into the 


laboratory frame yields 


pa ee Fe (lab) 


= S. + = 
ma ie, 1 su P34 ses 5 Gee? 


5 +eexp( 2 tiny talon, 


1 34 1 


0 0 
i6 exp{2ni(v, 1 +vot,.)} {exp2niv 


[2-33] 


where Tp! denotes the duration of the 90; S pulse. During 


the signal acquisition period tj, the density matrix 


elements develop further time dependence as 
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(lab) : : e 
+t. + = 1 
Pio (+ t +t itt.) ié exp{2ni(v,t tye Wha 


{exp(2niv,,t, ) - exp(2niv,,t, )}exp(2niv,5t,) F 


[2-34] 


(lab) i ee : 
P34 (ttt) +t. +t.) =a6 exp {2ni(vyt+ve t+) } . 


{exp(2niv,,t,) + exp(2niv )}exp(2niv ) 


1nd Baie 34°2 


Since the heterodyne receiver in an NMR spectrometer 
measures the Jijmagnetization in (the rotating frame, 575 


and p34 in the I rotating frame are required: 
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{exp(2nidv,,t,) + exp(2nidv,,t,)} x exp(2nidv,,t.) ’ 


where Avi is the frequency offset of the transition 


frequency Vig from the transmitter frequency Vite 


Rhe=ivspiny signal 1, + il jemeasured is given by 


y 


Signal @wiéfexp(2niAvj5t) + exp(2Qnifva,ty)) * 
exp{2nilvet,: - VET) } x [2-36] 
{exp(2nidv)9t9) a exp(2niAv3zy4t) } 


The §factor involving Th and Tp! is just a phase shift and 
this is compensated in real experiments by "phasing" the 
spectrum to get pure absorption mode presentation. Since 


this phase factor is the same for all signals, we can 


ignore it and write 
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Signal a id{exp(2nidv,5t,) + exp(2niAv3zygt,)} x 
[2-37] 


{exp(2niAvjot5) + exp(2niAvgyts)} . 


The signal expression in Equation [2-37] obtained 
from the above analysis, shows that the transition 
frequencies in the measured magnetizations appear as 
offsets from the respective transmitter frequencies. The 
Signal output by the simulation program is exactly similar 
to Equation [2-37], but the program represents the 
difference vy;-vy as vz (similarly vo-vg as vo). Hence vo 
and vy; are to be interpreted as offsets from the 
corresponding transmitter frequencies. The 
transformations of the density matrix between various 
frames of reference, described above, are not really 
necessary since any phase shifts built up in the duration 
of the pulses are unimportant. These phase shifts do not 
arise if all pulses are taken to be infinitely short so 
that laboratory and rotating Erames, coincide for “the 
duration Of a pulse. In this view, all calculations are 
performed in the laboratory frame, but the frequencies 
which appear in the signals must be interpreted as offsets 


from the respective transmitter frequencies. 
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Il.12 Other Features of the Simulation Program 


A facility to carry out the phase cycling of pulses 
and receiver is implemented in the simulation program 
since phase cycling is a common feature in most modern 
pulse NMR experiments. Phase cycling is usually employed 
to Suppress unwanted or spurious Signals in the 
spectrum. For efficiency in the simulation program, the 
density matrix data structure at the point in the pulse 
sequence where phase cycling begins is retained so that it 
can be reused in each pass through the cycle. At the end 
of each pass through the cycle, the measured 
Magnetizations are accumulated with some preliminary 
Ssimpada Eication. 

Another feature of the Simulation program allows 
explicit substitution for the duration of precession 
periods in terms of 1/(nJ), where n is an integer and J is 
the coupling constant between the I and S nuclei. The 
substitution feature is especially useful in pulse 
sequences having multiple refocussing periods, since it 
usually simplifies the density matrix considerably and 
keeps the number of terms in the density matrix elements 
within manageable limits. 

miesettect. of —Gt pulses of any) Given sintermediate 


durations can also be simulated by explicitly defining the 
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tip angles of certain pulses to desired values in the 
input.terhere,are..6 pulse ghengths: A, By, .C, Dp En .and-2, 
wach ght c90 75 B=, LOOn ~awhadenC .-D, «By and Fi can wbe sset to 
any desired value. This feature is useful in studying the 
effect of imperfect pulses, devising means to reduce the 
undesirable effects of imperfect pulses, and in analyzing 
experiments with composite pulses (46) or those 
experiments such as the DEPT pulse sequence (31), which 


involve pulses with intermediate flip angles. 


II.13 Limitations of and Suggested Extensions to the 


Simulation Program 


The model used for the simulation program excludes 
spin relaxation and chemical exchange processes. This 
limitation is reflected in the inability of the program to 
simulate experiments involving Nuclear Overhauser 
Enhancement effects or chemical exchange experiments such 
as 2D chemical exchange and saturation transfer 
experiments. 

Another limitation of the program is in the type of 
spin systems which can be analyzed. Only I,S, systems of 
spin 1/2 nuclei can be handled. This domain excludes many 
interesting spin systems such as A,M,xX; systems of spin 
1/2 nuclei and spin systems containing nuclei with spin 


greater than 1/2. 
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An extension of the program, which has been 
considered and for which preliminary versions have been 
generated, is the expansion of the set of analyzable spin 
systems to include a third spin: 1/2 nucleus, so as to make 
Am<nM. systems within the scope of analysis using the 
simulation program. Extension to include nuclei with spin 
greater than 1/2 has not been investigated, but should not 
present sgreateditficulty. 

Another modification would be the implementation of 
the feature of variable (user defined) phases of pulses 
and receiver, so as to provide phase values of 2n/k where 
k is an integer. Variable phases could then be described 
in the input pulse sequence by defining a reference phase 
and specifying phase increments. This feature would be 
useful in multiple-quantum experiments which identify 
particular spin systems by detecting the highest order 
coherence present in the system, after filtering out the 
lower order coherences. The shortage of memory space does 
become a limiting factor when simulating some long pulse 
sequences especially on large spin systems. It would be 
useful to introduce a facility to zero out certain density 
matrix elements which may not be of interest, when the 
arrays containing the density matrix get nearly full. 

With this feature it may become possible to simulate 
certain pulse sequences which are exceptionally long and 


Cannot be handled with the present program. 
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bi. at Summary 


In this chapter, the implementation of a simulation 
system for multiple pulse and 2DFT NMR experiments, on a 
digital computer has been described. The objectives of 
the program and details of the implementation were 
discussed. Brief descriptions of the basis sets used, the 
procedures for building up the representations of spin 
Operators and pulse transformation operators were given. 
The important details of the procedures for computation of 
the thermal equilibrium density matrix, the transformation 
of the density matrix during free precession periods and 
during periods of broad-band spin decoupling, and the 
simulation of the signal acquisition process were 
described. The identification of the chemical shift 
frequency terms in the output as offsets from the 
corresponding transmitter frequencies were demonstrated in 
the context of a simple example. Some of the other 
features of the program which allow phase cycling, 
substitution for the durations of free precession periods 
ana the facility to handle rf pulses Of any user-defined 
tip angle, were described. Finally, some of the 
limitations of the current version of the program were 


exposed along with suggestions for further extensions. 


enagnlesleelve Oia geri ath eit 
» is ATnon eae ow THuT Nee 8 
+0 bay Vevailae ay ‘Pit itnaed | wut ear ees it 


oie shen Dia ald 50 dnlgoah tow earv0%e . 


oY | a heehee aoe, aa= 7) «6G ning) sem apis0 9 : ae 
ins ow he peSpbauaet eabt ape. pi bint wt a THONG: - 
waTay e Pe ianbore paiayee aquest embtey low sos me 

» noire teuremp al! fee pomie, 4 bo os bedal airtel BOP 


Jp yee ttew yo Lantes ned Itt Dupe Sener wAy ia 
awit, hatin 23a Alouae ihe 


B24) 9 wt ¢ ~~ 3° eae 


rt 42 eis Haney ne i Can be abe bie", pean 7 


let [asinem eh e painghsinaagake ot? i json 


ony aoe: tee 1 UhigR pase af? nd swrtee 

cd hradeteseriaet. “wane ee ns inne 
{ore Gt? Tm ahr 

“inissyn eter le ape om 

aioas | my, AOD enwasr ty, e232 = nent d schon 

bean Tee a abe - te 


a 


ee 


The next chapter gives a detailed analysis of the effect 
of broad-band decoupling on the density matrix and the 
method of solution of the equation of motion for the 
density matrix in the presence of decoupler irradiation. 
Applications of the simulation program in the analysis “o£ 
multiple pulse and 2DFT experiments are discussed in 


Chapter iV. 
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CHAPTER Li1 


THEORY AND IMPLEMENTATION OF THE HETERONUCLEAR 


NOISE DECOUPLING SCHEME 


thie  elntrocjuct1on 


fhe description of the simulation program jin “Chapter 
II, mentioned the heteronuclear noise decoupling scheme 
briefly. This chapter discusses in detail the density 
matrix description of heteronuclear noise decoupling. 
Some of the earlier work on the theoretical analysis of 
the noise decoupling process is due to Ernst (47), who has 
given the equations of motion of the density matrix in the 
presence of decoupler irradiation, with implications for 
CW experiments. In order to improve the broad-band 
characteristics of the decoupler, several modulation 
techniques have been introduced (47-50). A multiple pulse 
decoupling cycle called MLEV, which greatly enhances the 
broad-band characteristics, has been devised (51). Pegg 
et al. (52) have presented a simple theory which describes 
the behaviour of 13c magnetization components when 


subjected to ly decoupling. This theory assumes the 
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decoupling field to be spatially inhomogeneous and reduces 
theseffect of the decoupler to a randomization of the 
direction ofethe proton spin vectors in a plane 
perpendicular to the axis of the decoupling field. 
Recently Levitt et al. (53) have demonstrated that 
unanticipated effects can be observed in certain 
experiments employing spin decoupling. 

In the early parts of this chapter, the equations of 
motion of the density matrix in the presence of a noise 
modulated decoupler field are set up and the method of 
sodutioneis illustrated ain détai 17forVthess spin 
decoupling in the IS, system and the final results are 
given for IS and IS, systems. Another type of decoupling, 
called the spherical randomization decoupling is also 
considered and the solutions for the density matrix 
elements are derived. Finally the implementation of the 
decoupling procedure in the simulation program is 


described. 


Lito orheoretical Treatment 


The time independent Hamiltonian for the I,S, system 
in a frame rotating at the frequency of the decoupler, va, 


can be written as 
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where vy; is the resonance frequency for the I nuclei, and 
Vor, the resonance frequency for the § nuclei is taken to 

be identical to the decoupler transmitter frequency. Jrre 
JIgsg and Jyog are the spin-spin couplings in Hz, Li and Sr 
are the dimensionless spin angular momentum operators for 
Enemy thet wnucdeus and) the p th S nucleus wespectively, I 


A A 


2 a, and S = ESp Are Une etoOtalsoinsOperators storm ethics! 
U 


nuclei and S nuclei respectively. 

The equation of motion for the spin density operator, 
an in the presence of a noise modulated decoupling field 
applied along the x-axis of the rotating frame of the X 


nuclei, has been shown by Ernst (47), to follow, 


do 


at = 2nilp,Ho]) oad Kiser ley 7c) ’ [3-2] 


where Keis»the spectral density of the decoupling field at 
frequency vo, and is proportional to the applied decoupler 
power. The operators are referred to the coordinate 
system rotating at frequency vo, about the laboratory z 


axils. 
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The equation of motion [3-2] can be represented 


the basis of eigenfunctions of Ho: as 


> 
where p is a row supervector containing the 
Eneedensity Matrix in wan appropriate order, 
supermatrix whose real part is proportional 


an imaginary part which is diagonal and has 


involve differences in the eigenvalues of Ho- 


elements 
and V is 
to K and 


elements 


in 


iea34 


a 
with 


which 


Details of 


the solution of the equation of motion are now given for 


the IS5 system. 


a. IS5 Spin System 


The eigenfunctions of fey for the iS spin system 


are simple products of the eigenfunctions |I,My> of I? and 


Z 


system, the basis functions are: 


T, and the eigenfunction |S,Mo> of Sand oy 


: For the IS5 
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The eigenvalue of Ho (the Hamiltonian in the frame 
rotating at vc) associated with the function [I,M;>|S,Mo> 
is -(vy-vg)Mz + JyoMyMco + Jggl2S(S+1)-3]/4. The functions 
|¢;> to |¢¢> belong to the Sj total spin S = 1 (triplet) 
manifold, while |$7> and |¢g> form the S, total spin S=0 
(singlet) manifold. Since there is no mixing of the 
singlet and triplet S5 spin states by any rf pulses or 
decoupling, the 8 x 8 IS5 problem reduces to a © x 6 
problem associated with the triplet manifold anda 2 x 2 
problem associated with the singlet manifold. It should 
be noted that the order of eigenfunctions of Ho differs 
from the order used in the simulation scheme (Section 
II.4) since it is more convenient here to have the 
blocking with respect to total S=1l and S=0. 

The equations of motion for the triplet manifold, can 


be written as 


) Kw 


d i 
ae 6117 P221P331P13° P31) = (01410 99133°P13°P31 col 


+ 2 1 J15(0,0,0, 0131-031) ’ [3-5] 


) Kw 


d —_ 
ae (912121 7° 23/932? = (015191 °P23'P32) KWo 


+ 2ni Jz5(p12+-P21+P23"—P32) | [3-6] 
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d 
de P14 1P 251736 7P 167734) = (P41 P 251% 36 1P 16134) 


[KW)-2ni( vy-vg) 1] ate 2n1 JtgP14191-p3670,0) t [3-7] 


) 


d 
de! 157024 'P261P35) = (Py5%Po41Po61P 35) * 


[KW) = 12n(vq-vg) 1] + 274 Jyo(P 51> 9475P967P 35) ’ 


[3-8] 
So '.441055 1066 P46 "0 64) = ((44°P55°P66"P46°P64) KW) 
+ i2m Jyo(0,0,01-p 467064) [3-9] 
and 
Se (P45 1954 1056" 65) = (0451054705670 65) KWo 
+ 12m Jygl-pggrP54r-P56°P65) * (3]001 


where 1 is a unit matrix of appropriate dimension, 
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-l iI 0 =17.2 -1/2 
tl -2 1 1 i: 
Wt = 0 il -) = 2 -1/2 ’ [3-11] 
-1/2 1 -1/2 -l 0 
-1/2 i -1/2 0 —] 
and 
-3/2 I i -1/2 
} -3/2 -1/2 i 
W5 = e [3-12] 
a 1 -1/2 -3/2 1 
— 1/2 1 l Aa 


The equation of motion for an element below the diagonal 
which is not coupled to elements above the diagonal (e.g. 
P41) is just the complex conjugate of the equation of 
motion for the corresponding element above the diagonal 
(e.g. py4)- 

We are interested only in the case of strong 
decoupling, so we shall assume that the decoupler power is 
sufficiently high that K >> Jyc, so that the terms 
involving Jyg in Equations [3-5] - [B= Vor canebe 
neglected, and each of these equations reduces to the form 


>= 


Se = OEKW 42h oa 1) + p33] 
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where W is a real symmetric matrix of numbers and A is a 


frequency variable (which may be zero). The general 


SOduUclTON LCOPhquation [3—12)) 1s 


pqi(t) = ) Dig oe OR) a ok Ole CxXD LO KV) Coos 
narre bte: 

where p,j(0) is the kl element of the density matrix at 

the beginning of the decoupling period. U and A are the 

eigenvectors and eigenvalues of the supermatrix W and 


~ 


Satisfy the eigenvalue equation 


ei Week, [3-15] 


with A} a diagonal Supermatrix. The index a in Equation 
ip) 4)eruns over ali columns Of | U (i.e. all pairs sot 
Mmovces 17 OL Pi3) The matrices W will always have 
negative or zero eigenvalues. This means that the terms 
in Equation [3-14] for the negative eigenvalues will have 
exponentially damped oscillatory time-dependence, while 
those terms in Equation [3-14] associated with zero 
eigenvalues simply oscillate at frequency A with no 
damping. Thus, at long times, the transient effects which 
Gharacterize the initial part of the decoupling period 


will have died away and the density matrix elements will 


have obtained "steady state" values given by 
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mcMid= a0) <P exp (-2nists, (32169 


where the elements of the supermatrix P are defined by 


Deke Ge Uk en a [3-17] 
04 


with the sum over qa restricted to those columns of U 
associated with the zero eigenvalues of W. It is the 
density matrix given by Equation [3-16] which 
characterizes the spin system at the end of a period of 
decoupling. If the duration t of the decoupling period is 
too short, or the decoupler power is too low, this “steady 
state" of the spin system will not be attained during the 
decoupling period. 

bheishould besneted that seacheblockweof fdensity matmix 
elements which are connected by the decoupling interaction 
can be characterized by a constant value of AM, = Mri-Myy, 
where M;; is the eigenvalue of i for the basis function 
1o5> (see Equation [3-4]), for all members P45 within a 
given block. Hence the frequency A in Equations [3-13], 
[3-14] and [3-16] associated with a particular block of 
interconnected elements is equal to AM, (vy-vg)+ 

For the triplet Sj manifold of the IS, Spin system, 


Bquations sia-o1y, (3-7) and [3-9] reduce to, [3-13 )ewith Wo= 
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Whe and with A = 0, (Vrove) mand). 0, respectively. “The 
scteadya state solutions fon each of, the b locks, of ‘5 


connected elements {p11 10991P331013"P3)} 6 


{e147P257P36°P167P34} ANd {pygressrPEgrP4grPg4} are Given 
byoEquation [3-17] with Pp = Pir 


3/8 1/4 3/8 -1/8 -1/8 
1/4 ae 1/4 1/4 1/4 
Py = 3/78 ibys 378 ey -1/8 ' 13— 237) 
=1/8 1/4 -1/8 37.8 3/8 
-1/8 1/4 =i1.8 3/8 3/8 


and with A = 0, (vj-vc) and 0 respectively. Equations 
[3-6], [3-8] and [3-10] reduce to Equation [3-13] when K 
>> Szg, with W = Wo, and with A = 0, (vy-vg)s 0 
respectively. The steady state solutions for each of the 
blocks of 4 connected elements {P191P211P23°P32}8 
{9151P247026135} aNd {o451P54 056 1PG5) Are Given by 


Bquakion [8-17] withy P= Boy 


Po = (1/4) [3-19] 
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andy with Av'= 0; (vy-Vvo) and 0° réspective ly. 
The equation of motion for the density matrix 


elements associated with the singlet 59 Spin states are 


Sag SOG 


jue cere = Oe 

and [3-20] 
ooo 2 

de 2 YI 9s) P78 


with, solutions (Equation [3=17]'’ with® Pp = 1) 


p77(t) ae p77(0) ’ 
Pag) =| pg3(0) (P=007 


and py 7g(t) = p7g(0) exp[-i2n(vy-vo)t]- 


These elements are not influenced by the decoupler in that 
the singlet S5 state is non-magnetic so that the singlet 
manifold of IS, behaves exactly like an I, Spin system 
does under § decoupling. 

Equations [3-17], [3-19] and [3-21] give the density 
matrix elements in the frame rotating at frequency vo 
about the laboratory z-axis. It is more convenient to 
have the density matrix elements in the laboratory frame 


since they display the full frequency dependence of the 
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system. The transformation from the Vg rotating frame to 

the laboratory coordinates can be effected using the 

relationship in Equation [2-26]. Taking into account the 
eal ele : 

frequency dependence of p(t) in Equation [3-16], one 


obtains 

lap 
Puce pCO) Pa S(t), [3-22] 
where S(t) is a diagonal supermatrix with elements 
534715) = exp{2nilvy(My5-Myzq) + vg (Mg 4-Mgj) I]t} 5 [3-23] 


In writing Equation [3-23], it has been assumed that the 
laboratory and rotating frames coincide at time zero. The 
frequencies in the exponential arguments in Equation 
[3-23] are the frequencies for the transition from state j 
to state i for the case Jp cue 0, i.e. they are the 


"decoupled" transition frequencies. 


De ts end 1s, Spin Systems 
For the IS spin system, with the eigenbasis labelled 


as, 
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[3-24] 


the equations of motion (Equation [3-3]) can be decomposed 
into sets of coupled differential equations with the 
SIE) AE) Ne Maes oO UL Oi a Molle i 

{9331044} and {9341043} - ine the Pilmie of Kk Jtg, each 
of these coupled sets of 8 equations reduces to Equation 


[3-13] with W = W3, where 


~~ 


Nl) [3-25] 


and the frequency A = AMIN VY T= Vici with AMy characteristic 
Ore thempartiwcular palr Of Coup led edements. “The steady 
state values of the density matrix elements in the 
laboratory frame are given by Equation [3-23] with P = ab 


where 


Beal) [3-26] 


For the IS3 system, the eigenfunctions of Ho can be 
grouped into irreducible sets belonging to total S spin 


a7 27/72) and)1/2". “the basis functions used are then, 
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foes elleeel/2> 13/2,3/2> o lés> = |1/2,-1/2>|3/2,3/2> , 
\Orralen ail erl/2>13/20172> idee = l72,-1/2>13/2)1/2>5 3 
lo3> = |1/2,1/2>|3/2,-1/2> , lo4> = |1/2,-1/2>|3/2,-1/2> , 
lo4> = |1/2,1/2>|3/2,-3/2> , ldg> = |[1/2,-1/2>|3/2,-3/2> , 
[3-27] 
bOg >a =o/l/241/2>|1/2/1/2> », [o44> = |1/2,-1/2>|1/2,1/2> , 
loyq> = [1/2,1/2>|1/2,-1/2> , |,9> = |1/2,-1/2>|1/2,-1/2> , 
Oia =e17271/ 21/2" 1/2>, loy5> = |1/2,-1/2>|1/2',1/2> , 
loq4> = |1/2,1/2>|1/2',-1/2>, |6,¢> = |1/2:-1/2>|1/2',-1/2>, 


where the prime is used to distinguish the two irreducible 
sets Of functions with total spin 1/2. Since rf pulses 
and decoupling do not cause transitions between levels of 
different irreducible sets, the IS3 problem can be broken 
down into one 8 x 8 and two 4 x 4 problems. Each of the 4 
x 4 blocks of elements of the density matrix in IS3 behave 
in exactly the same way as the elements of the IS density 
matrix given above, so they are not considered further. 
messes ais; Dulocksof the 1S, density matrix associated with 
the total spin S = 3/2 manifold has equations of motion 
which can be grouped into coupled sets for the ordered 


groups of elements: 
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{0111022 4P33"P446P 13131 P 241 4a}e 
{214 1P23 1°32 "P41 10129034 1P 21 P43} 8 
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pnw Gae™ aml tO kK >> Jtg, Gach of these coupled sets of 


equations reduces to Equation [3-13] with W = War where 


3/2 es ee 0 -/3/2 -/3/2 0 0 
3/2. =7/2 2 0 /3 V3 -¥3/2. --73/2 
0 2 -7/2 3/2, -(3/2. -/3/2. 73 V3 
0 0 3/2 =3/2 «0 0 -/3/2 -/3/2 
aaa ee = , [3-28] 
e -/3/2 3 -/3/2 0 -5/2 0 8/20 
(3/2. V3 -/3/2 0 0 -5/2 0 3/2 
Oe /3872 3-8/2 372° 0 -5/2 0 
0 -/3/2. 73 -73/2 0 3/2 0 -5/2 


and the frequency A = AMy(vq-Vo)s with My characteristic 
of the particular set of eight coupled elements. The 
steady state density matrix in the laboratory frame is 


given by Equation [3-23] with P = Py, where 
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[16-29] 


c. Spherical Randomization Decoupling 

Spherical randomization decoupling, introduced by 
Pegguet aleu(52), differs from normal noise decoupling in 
that the phase of the decoupler is randomly switched among 
Ene -Orywo0°,mls0% and 270° ‘channels by the noise 
modulation rather than between the 0° and 180° phases 
normally used for decoupling. The equation of motion of 
the density operator during spherical randomization 


decoupling is 


= ilp,Ho] == rece ore 3 a8 eevee [3-30] 
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which differs from Equation [3-1] only by the presence of 


the double commutator involving Sy. There are no cross 
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terms of the form [Sy [8),0]} in Equation [3-30] because 
the cross correlation function for x- and y-components of 
the decoupling field vanishes. Equation [3-30] can be 
represented 1n) the form “of Bquation [3-3] and in the limit 
Of iK >> Jis it will break down into a number of 
collections of elements which satisfy Equation [3-13] and 
whose steady state values in the laboratory frame will be 
given by Equation [3-23]. 

For the IS system represented in the basis (Equation 
[3-24]) the sets {P119Po0}% {p13°P 74} amd {p337P4q4} will 
have associated with them Wa and .P2 > ust. as they did ‘for 
normal decoupling. The sets {pj9-p9 1} and {py4rP93} 6 


however, are associated with 


Ne ie : [3-31] 


Since We has no zero eigenvalues, the corresponding P 
Supermatrix is the null supermatrix, which implies that 
Pilar Paz" P14 aNd p93 vanish during the decoupling 
process. Thus the only off-diagonal elements of ° which 
remain finite during decoupling are pj3,r poq4r P31 and pyy 
which correspond to the perpendicular I magnetization. 


For the IS, system represented in the basis (Equation 


{3-4]) the sets of elements {P7170 227P337P13°P31)" 
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{61472510367 161P34) ANd {pggrP551P6E104G1PG4) have 


associated with them a W sSupermatrix: 


ay 2 0 0 0 
2 -4 2 0 0 
We = 0 2 =2 0 0 ' Los 
0 0 0 2 0 
0 0 0 0 =2 


Pee awl 3) cjele hl els yO m0 [3-33] 


Pe is quite different from the corresponding supermatrix 
Py for the normal decoupling case. The 4-member sets of 
elements for IS9: {1910211237932} 4 {P151P24"P26"P35) 4 
and {p45705470567P65) Satisfy Equation [3-13] with Waa 


where 
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which has no zero eigenvalues so that the corresponding Pp 
Supermatrix is the null supermatrix and all of these 
elements are destroyed by spherical randomization 
decoupling. The elements p77, p7g and pgg of the IS, 
density matrix are not affected by either type of S- 
decoupling. 

For the spin S = 3/2 block of the IS, system 


represented in Equation [3-27], the sets of elements 
{011 10221P331P44 P13 70311924 9P 42) 8 


{015 7P261P377P48°P179P357P 289046) ANG 


{9551P661P77708810571975'P68'PR6) have an associated W 


Supermatrix: 
23 5 0 0 0 0 0 0 
3 = 4 0 0 0 0 0 
0 4 =a, 3 0 0 0 0 
0 0 3 iS) 0 0 0 0 
Wee 5 ; , ay ge: ; . ab [3-35] 
0 0 0 0 0 =] 0 3 
0 0 0 0 3 0 =, 0 
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which has no zero eigenvalues so that the corresponding P 


Ssupermatrix is the null supermatrix and all elements in 


these collections vanish with spherical randomization 


decoupling. 


iis 


General Effects of Noise Decoupling on the Density 


Matrix 


The results for normal decoupling show a number of 


interesting features: 


al) 


Gen) 


(d) 


the density matrix elements at steady state are 
linear combinations of appropriate collections of the 
elements of the density matrix at the beginning of 
the decoupling period, with "decoupled" frequency 
dependence (see Equations [3-22] and [3-23]); 

the linear combinations of elements of p(0) which 
MakewUpealpatblLicular steady state element ,oL 35) fare 
not simple averages in most cases; 

the decoupling mixes diagonal elements of e(0), which 
have zero frequency dependence, with 2 quantum 
coherence elements, which oscillate at t2Vo in the 
laporaLlory trame-; and 

it mixes elements of p(0) associated With transwuvons 


of the nuclei (frequency tvy7) with multiple quantum 


coherences (frequencies +(vyz+2vg) and +(vy-2vg)). 
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These features are in marked contrast to the results 
for spherical randomization mode decoupling where the 
elements of o(t) are just simple averages of elements of 
p(0); there 1S no mixing of "diagonal elements with 2 
quantum coherences nor is there mixing of the elements 
associated with I nucleus transitions with any multiple 
quantum coherences. Those elements of o(t) which 
correspond to S nucleus transitions and to multiple 
quantum coherences are driven to zero by the spherical 
randomization decoupling field. The decoupling field in 
the usual noise decoupling experiment jumps randomly from 
the +x.) to ‘the’ —x “direction ain the Vo rotating frame so 
that the x-components of the S-magnetization and the real 
pacts Of 9(0) associated with multiple quantum coherences 


are unaffected by the decoupling field. In the spherical 


randomization mode, both real and imaginary parts of the 


off-diagonal elements of p are modulated by the decoupling 


~~ 


field and those off-diagonal elements for which AMg # 0 


are destroyed during the decoupling period. 


Tershould be noted that the analysis in «Section Ill. 2 


does not assume spatial inhomogeneity in the decoupling rf 


field, which is a requirement in the theory proposed by 
Pegg et al. (52). It is important that a theoretical 
description of decoupling be applicable in cases of low 


inhomogeneity as well as large inhomogeneity in the 
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decoupling rf field, since in modern spectrometers the 
transmitter coils are so designed to minimize spatial 


inhomogeneity of the pulse. 


III.4 Implementation of the Decoupling Scheme. 


The ordinary noise decoupling scheme discussed in 


Section III.2 is implemented as the subprogram DECOPL in 


the simulation program described in Chapter II. One of 


the major considerations in the implementation of the 


decoupling scheme is the identification of the collections 


of density matrix elements which are coupled by the 
equation of motion (Equation [3-2]). The double 
commutator in Equation [3-2] causes the coupling of the 
density matrix elements. The i,j matrix element of the 
double commutator is 


A A A ~ ? ; 
Eigse Ae SyJi5 = BCSy ) ik Pk 3 =2 Ne SBE oes als 


k k, 


= z 01x (Sx) [3-38] 


and shows that the element Pi3 is connected via the 

decoupling perturbation to those elements Pk for which 
. 2 

(onciae # 0, to p3, for which (S, ey # O and to those 


elements py, for which (Sx)i,(Sx)15 70 eel Order ato 
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generate the lists of connected elements, the matrix 
representations Sx and oe are first generated. The set 
of elements coupled to a particular density matrix element 
Pijr are found by scanning the density matrix for elements 
Pkjr Pk1 endiipsj LO themsersmor indices! k= witch give 
non-zero contributions to the right hand side of Equation 
[3-38]. The lists of coupled elements so generated are 
arranged in an array with appropriate labels, so that the 
number of elements in each coupled group can be easily 
ascertained. The group of connected density matrix el- 
ements, which occur in the inom eystems With 1) <i ne <3, 
can be classified according to the number of elements in 
the list: 1l-element, 2-element, 4-element, 5-element and 
8-element lists. During a period of spin decoupling, the 
density matrix elements in each set transform among 
themselves (Equation [3-16]), the transformation 
coefficients being given by the P matrices in Section 
ipso.) The D matrices, corresponding to” the different 
lists are stored in the program so that the appropriate 
linear combinations can be computed. When the lists of 
coupled elements of the density matrix have been 
generated, the subprogram scans each list and determines 
the appropriate P matrix to be used in the transformation 
Seecilements oLtthe particular ist... To carry out. the 


transformation, the coupled density matrix elements 
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multiplied by appropriate elements of the Pp Supermatrices, 
are copied into a work area in the “data array". 
Bookkeeping information concerning the density matrix 
elements thus computed is then updated. 

In the laboratory coordinate system, the off-diagonal 
elements exhibit "decoupled" frequency dependence, i.e. 
just the Zeeman frequencies of both nuclei. The 
"decoupled" frequency is represented in the simulation 
program by specifying the frequency dependence in a manner 
analogous to that for a free precession period (Chapter 
II) with an additional variable set to indicate that the 
Spin-spin coupling constant Jts is to be set to zero in 


all frequency variables for this decoupling interval. 


III.5 Summary 


A density matrix description of heteronuclear noise 
decoupling has been presented in this chapter. The 
equations of motion of the density matrix were set up and 
the method of solution described for the case of strong 
decoupler power. Two types of noise decoupling procedures 
were considered: (i) normal noise decoupling or planar 
randomization decoupling and (ii) spherical randomization 
decoupling. The implementation of the decoupling scheme 


in the simulation program was discussed. The next chapter 
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wilustrabes the application, of the simulation program by 
analyzing various pulse sequences on representative spin 
systems. Included in the applications is an analysis of 
the decoupling experiments in which unexpected features 


are observed. 
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CHAPTER IV 


APPLICATIONS 


DV neroduction 


In this chapter, the capabilities of the simulation 
program are explored by including the results for 
Simulations of several pulse sequences applied to a 
variety of spin systems. The examples have been chosen so 
as to illustrate various features of the program. The 
experiments simulated here include the Distortionless 
Enhancement by Polarization Transfer (DEPT) (31), a two- 
dimensional version of the DEPT experiment for selective 
ebservat tonyor individual CH, systems @(54-57)) , the 
multiple-quantum coherence experiment INADEQUATE (58,59) 
and versions of the 2D gated J-§ correlation experiment 
which illustrate certain unexpected effects of spin 
decoupling (53). Many of the experiments mentioned above 
are not easily understood in terms of the motion of the 
magnetization vectors, since they involve pulses of 
intermediate flip angles, excitation and evolution of 


multiple-quantum coherences and so on. Throughout this 
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Chapter, the notations: used for certain terms have been 

modified so as to correspond with the convention used in 
program output. The vy and vo terms shall be denoted by 
VI and VS respectively and the spin system I,S, shall be 


denoted by ImSn. 


IVe2 Simulation of Experiments 


a. DEPT Experiment 


The DEPT sequence, shown in Figure 3, is an ingenious 
method for transferring polarization from one spin system 
to a heteronuclear spin system without phase or intensity 
distortions (31). MDoddrell et al. (31) have described the 
behaviour of 13cu, {CH and ESCH fragments, when DEPT is 
applied, using Heisenberg vector diagrams, and have 
demonstrated some of the spectral editing capabilities of 
this sequence. 

The most popular application of the DEPT sequence is 
to CH, systems. The decoupling during the acquisition 
stage of the DEPT sequence (Figure 3) is optional. The 
relative intensities of coupled carbon signals from DEPT 
sequence for CH, CHy and CH3 groups, are the normal 
Characteristic ratios 1<1, 1:2:1 and P:3:331 
respectively. In contrast, the relative intensities from 


a refocussed-INEPT experiment (60,61) are distorted and 
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DEPT PULSE SEQUENCE -EOREG-He SYSTEM 


COMM ** DEPT SIMULATION FOR CH SYSTEM *» 
COMM ** THETA PULSE=45 «+ 

SPUNT S41 

SH CUES EC 45.0 

WRDM 
PULSSA+X 
WRDM 
PRECD1 
SUBSD Te 1/1) 20) 
WRDM 
PULSSB+X 
PULSIA+X 
WRDM 
PRECD1 
SUBSD1=1/( 2u) 
WRDM 
PUESSC=¥ 
PULSIB+X 

WRDM 

PRECD2 
SUBSD2=1/( 2d) 
WRDM 

DCPLS 2 
WRDM 

MEASIXYO 

STOP 


pop = - == 
=4==25-- 


Figure 3. Pulse sequence and input data for simulation of 
the CH DEPT experiment. 
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are srespectively Isl; :O0:ls lsl:ls1. The intensities of 
the signals detected after the DEPT sequence depend on the 


Flip angle of the proton 6(-Y) pulse. For CH carbon 


Signals, the maximum enhancement occurs at 6 = 7n/2, CH5 
carbon signals are largest at 6 = n/4 and CH3 carbon 
Signals have the maximum enhancement at 6 = 0.196n. The 


dependence of the signal intensities on the angle 9 for 
CH, CHs and CH, systems are (31): sin ©, Sin 2e67and 
(Sin 6 + Sin 36) respectively. Spectral editing is 
accomplished by generating individual CH, CHy and CH3 
Subspectra from appropriate combinations of spectra 
recorded at 6 = n/4, n/2 and 3n/4. Spectral editing using 
the DEPT sequence is more exact than that using the 
refocussed INEPT sequence, since, in DEPT, the signal 
intensity depends mainly on the pulse angle and the J 
dependence of the signal intensity is a second order 
effect much smaller than that for INEPT. DEPT has: another 
advantage over INEPT, in that it employs fewer pulses, so 
that DEPT is less vulnerable to intensity variations 
arising from inhomogeneous or imperfect rf pulses. 

The input to the simulation program to perform the 
DEPT sequence with 6 = 45° on the 11S] (13cH) system is 
shown in Figure 3 and the output in Figure 4. The 
diagonal elements of the density matrix are populations of 


the spin levels: pj and p34 are proportional to the 
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DEPT CALCULATION FOR CH SYSTEM : THETA PULSE=45 


Initial Density Matrix 


DG tne i =se 52. 5008 0.0 ) 
DMI 2 we ues ain SOOk O20) 
Die Sa Sua = 47 SOOn a} 10) 
DMG Seale = iliee='52 500) Om) 
Density Matrix After 90(Proton) Pulse 
OL t= tlk, WAN ae <atl 200k ON.00 
oh Fah Sh Ge se OO 5.0 0008) 
DMRS 2 2) eae me On OO} 
DME 254i) es <3 OOF 501000)) 
DME 35) oi) si 2.500, CFO 5) 
OM ( 4, 4) = +( -2.500, (a) ) 


Density Matrix After First Precession Period (D1i=1/(2*uU)) 


Wid IE thy a ee eal 2.500, 0-0} 
DMS resis tl 5.02 000) 0.000)*EXPI](2*Pl*( 1.00 VS )=D1) 
DMN 22) = 25007 On Ome) 
ahs Wh ey Gi) sat Ryo sOlale). 0.000)*EXPI(2*PI*( 1.00 VS )*D1) 
Yo ll eb ehh Se AM 2.500. 9) 
OMe aie Sica + (te 2 500) O05) 


Density Matrix After 180(Proton)-90(Carbon) Pulses 


yi Sh eal) Sal 0.0 2.500) 
DM ( 1, 4) = + 0.0 -50.000)*EXP1(2*P1={-1.00 VS )=D1) 
DMN Zio ieee tal 0.0 50.000)*EXP1(2*PI*(-1.00 VS )*D1) 
OM Sie Au =i 0.0 2.500) 


Density Matrix After Second Precession Period (D1=1/(2*J)) 


ley f ahg Ail eit Sel Gyoyol 2 0.000)*EXPI](2=*Pl=( 1.00 VI )=D1) 
OM (1, 407 = +4 0.0 , -50.000)*EXPI(2=*Pl1s( 1.00 VI )*D1) 
DM =f -2 Sih eat 0.0 , 50.000)*EXP1(2*Pl*( =12100) Vi )*D1) 
DM) 3), 4) = i Ze 500K 0.000) *EXPI(2*Pl*=( 1.00 VI )*D1) 
Density Matrix After Theta! Proton) -180'!Carbon! Pulses 
OM 1, 2I= +0 A768), =85 2355 )/* EXP 1(2=P i= -1.00 VI )=D1) 
DM iia ten =e 768 C5. SOSNm EX PH =P) lan ( =O OF Nil y*D1) 
DMI 25S) ee = OB ars Sis 55)) SEXP li 2= Pi si 1.00 VI )=D} 
BLN (thy CN) eM 1.768, 35.355)*EXP](2=Pl*( = 00) Wil )*D1 
Density Matrix After Third Precession Period (D2=1/(2*u)) 
i {i As A) Soca eI iS, Sal (Asreh OS Mize ea = 00) Vi} )=D1 
PEXPI(29P1*( 1.00 V1 )*D2 
DMS Min 4 eat OC So eSoOnse Keno s Pos 1 OO Vil )=D1 
PEXPI(2*°P]*( 1.00 VS 1.00 V1 )*D2) 
OME 25) 3) "=" +> = 1e768) s-sSn eb) EXP 2 spi= 1.00 VI )*D1) 
*EXP](2=*Pl*( 1.00 VS -1.00 VI )*D2 
DMs 4st SOMOS LOCUM EAP lire Plait -1.00 VI )*D1) 
PEXP](2*Pl* 1.00 VI )*D2) 
Density Matrix During Decoupled Acquisition 
my € i, PS ot Cake, oils Std rou TOR Vt )*T2) 
PEXPI](2*Pl*( =41 00 VI yeD1) 
PEXPI(2*Pl*( OO) Vil )*D2) 
OM ( 1, 4) = +( -0.684, 17. 678))*EXPI(2*Pr*( 1.00 VS 1.00 V1 )*T2) 
*EXP] (2=P1*( =i) 008 Nid yeD1) 
*EXP](2*Pl*( 1.00 VS 1.00 v] )*D2) 
+( -0.884, 17.678)*EXP1(2*Pl*( 1.00 VS 1.00 V1 )*T2) 
EXP] (29Pl1*( = 00) V1 )*D1) 
PEXPI](2*P]*(-1.00 VS 1.00 V1 )*D2) 
DM ( 2, 3) = +( -0.884, -17.678)*EXP1(2*Pl*( 1.00 VS -1.00 V1 )*T2) 
*EXP1(2*Pl]*( TOO RVa )*D1) 
*EXP](2*eP]*(-1.00 VS -1.00 VI )*D2) 
+( -0.884, -17.67B8)*eEXPI1(2*Pl*( 1.00 VS -1.00 V1 }eT2) 
SEXP] (2*P1*( 1.00 V1 4*D1) 
PEXPI(2*#P}e( 1.00 VS -1.00 VI eer 
=) +f .355, -1.76B8)%EXP1(2*P1*< 4.00 V1 ° 
ee eg ®EXP1(2=P l= ( 1.00 VI )=D1) 
*EXPI](2eP1*( 1.00 V1 )*D2) 
Carbon-13 Signals During Acquisition 
: +( 70.710, -3.536)eEXPI(2*Pie{ 1.00 V1 )*T2) 
are ot 4 : *EXPI](2*Pl*( -1.00 VJ beD1) 
EXP] (2*Pl*( 1.00 VI )eD2) 


Eiguse 4. (Output of the DEPT simulation for CH system. 
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complex magnetizations (M, + iM.) associated with the I 
resonances at VI + J/2 and VI - J/2; p 3 and poq4 are 
proportional to the complex magnetizations associated with 
the S resonances at VS 4+.9/2 and vS = J/72: P14 Bepresents 
a two quantum coherence at frequency VI + VS; and P93 
represents a zero quantum coherence at frequency VI - VS. 
As shown below, the density matrix information output by 
the program can be translated into a collection of 
magnetization vector diagrams for visualization of the 
magnetization transfer process. The magnetization 
diagrams for the DEPT experiment on the 13cy system are 
Givenwin Eiqube (5. In this figure, the axes ore the 

doub le-quantum and zero-quantum rotating frames represent 
the real (X) and imaginary (Y) parts of the corresponding 
density matrix elements. 

The step-by-step time evolution of the density matrix 
can be followed by translating the output in Figure 4 to 
conventional representation. The initial thermal 
equilibrium density matrix in the simple product basis 
fO/ 2/2 eA 22g l/ 241) 2> feo =e 
hyp 72> |e ly 2omand |/1/2,—1/2> | 271 2>) eo avena in 


Figure 4, can be represented as 
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6+A 0 0 0 

0 —-h+A 0 0 
OO) mea Gl 722 ) ’ [4-1] 
% 0 0 5-A 0 

0 0 0 -6-A 


where A = YqNH G/8kT, § = Y¥chHo/8kT and yur yor Ny Hor k 
and T have their usual meanings. A and 6 have been 
assigned numerical values 100 and 5 respectively in the 
Simulation program. The proton magnetizations are put 
into the x-y plane by the proton 90° pulse (Figure 5A), 
and they precess at frequencies VS + J/2 in the frame 
rotating at the proton transmitter frequency. At time 
Dl = 1/(2Joq), the phase difference between the two 
Gomponents 1S 1e0e (rigure 5B).)) At this point, the 


density matrix is given by 


6 0 a 0 

0 -6§ 0 -a 
eCDi) =" (172) ' a2 
5 a* 0 6 0 

0 -a* 0 -§ 
Were ans A eExpiena VS Dit. The carbon 90° pulse 


transfers the proton magnetization into the zero-quantum 
and double-quantum coherences, and puts the carbon thermal 


equilibrium z-magnetizations into the perpendicular plane 
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where they rotate at frequencies VI + J/2 in the frame 
rotating at the carbon transmitter frequency (Figure 

5C). During the second precession period, Dl = 1/(25 oy), 
the carbon magnetizations acquire a phase difference of 
180°. In this period, the double-quantum coherence has 
acquired a further phase angle 2n(VI + VS)Dl in the frame 
rotating at the sum of the frequencies of the proton and 
carbon transmitters (Figure 5D). During this period, the 
zero-quantum coherence acquires an additional phase angle 
of 2n(VI - VS)D1l in the frame rotating at the frequency 
difference between the proton and carbon transmitters. It 
should be noted in the density matrix, that the purpose of 
the proton 180° pulse is to have the 2n VS Dl phase angles 
in the proton magnetizations refocussed in the zero- and 
double-quantum coherences in the second precession 

perroc. | The density matrix at the end of “the second DI 


period is given by 


0) Cc 8) b 

c* @) -b* @) 
o.(2D1) = ’ [4-3] 
as O -b @) -c 

jens O -c* @) 


WieresDl=r—-IN expy2ni Vi DI} and ¢°= 6 expy2nrisVvl Dig. 
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The proton 45° pulse transfers some of the zero- 
quantum and doub le-quantum coherences into carbon 
magnetization components (p19 and p34) which have a phase 
diusberence! Of 180° (Figure 58)i.,4.these: pulses. also 
transfer some of the transverse carbon magnetization 
components into the zero- and double-quantum coherences 
(p53 and py4) where they remain during the detection 
period. The purpose of the 180° carbon pulse applied 
along with the proton 6 pulse is to effect the refocussing 
of the carbon resonance offsets. During the third D2 = 
1/(2J3ey) precession period, the transferred carbon 
Magnetization components are completely refocussed (Figure 
5F). The phase angles of + 2n VS Di built up on the zero- 
and doub le-quantum coherences during the second precession 
period srandeconverted sinto wphasewang les) of —27, VSeD1 “tor 
each of the carbon magnetization components are exactly 
cancelled by the 2n VS D2 phase angle acquired during the 
third precession period. The density matrix just before 


the signal acquisition stage is given by 


0 d 0 e 

a* 0 e 0 
O20 3D2)) = ' [4-4] 
~ 0 e* 0 d 
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where d = (1//2)(A-i6) and e = (1//2)(-8+tiA) 

exp{2ni VS D2}, since the time period Dl is identical to 
D2. The output in Figure 4 does not show the complete 
distortionless magnetization transfer as clearly as one 
might expect because of the program limitation which 
restricts the appearance of any time variable more than 
twice in a pulse sequence. 

When the decoupler is turned on during the 
acquisition period, both p}5 and p34 precess at the 
decoupled frequency VI and the signal measured with a 
quadrature detector corresponds to a single resonance at 
frequency offset VI from the carbon transmitter with an 
intensity which corresponds to a mixture of transferred 
proton magnetization and the normal Boltzmann carbon 
magnetization. The signal expression output by the 


program can be written as 


Signal = (1//2)(A-i6s) exp{2ni VI (D1-D2)} exp{2ni VI T2} 
= (1//2)(A-i6) exp{2ni VI T2} [4-5] 
since Dl = D2. The Boltzmann carbon magnetization, which 


correspond to the 6 term in the output can be removed by 
application of a carbon 90° pulse at the beginning of the 


DEPT sequence or by alternating the phase of the proton 6 
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pulse by 180° on successive acquisitions. The signals are 
independent of the proton and carbon transmitter offsets 
during the free precession periods Dl and D2. This offset 
independence feature gives rise to distortionless 
polarization transfer, since the chemical shifts are all 
rerocussed, regardless of their offset from the 
transmitter frequencies. 

The coupled spectrum using the DEPT sequence has the 


carbon signal given by (Figure 6) 


Signal = (1/2¥2)(A-i6d) [exp{2ni(VI+3/2)T2} + 


eExpie na (Vl) 2) Toh ley [4-6] 


which corresponds to two lines at frequencies VI + J/2 and 
ViEemuyzyewith bs lere lative: intensity, 

The dependence of the signal intensity on the pulse 
angle 6 has been studied by carrying out the density 
matrix calculations with various values for the 9 pulse. 
Thewsimulation output) is given in higure 7. ror 
transferred polarization from protons, the Signal 
intensity dependence on the @ value corresponds to the 


general expression 


Signale=wA sing expt2ni. VI 12} [4-7] 
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INPUT FOR COUPLED DEPT EXPT ON CH SYSTEM 


—— 


COMM COUPLED DEPT SIMULATION FOR C-H SYSTEM 
COMM THETA PULSE=45 

SPINI1S1 

SET FPULESC 45.0 

PULSSA+X 


UESSCSY 
ULSIB+X 

RECD2 
UBSD2=1/( 2u) 
PREC HZ 
MEASIXYO 

STOP 


P 
P 
Pp 
SUBSD1=1/( 2u) 
Pp 
P 
P 
S 


OUTPUT FOR COUPLED DEPT EXPT ON CH SYSTEM 


Carbon-13 Signals During Coupled Acquisition 


SIGNAL 


Figure 6. 


= al Sey Ss seA Nee i 
*EXP1(2*P I ( > il oe) Wat 
ENP 2 PT 1 


il Sane = ie MaRS ES ey Ae esc Wc 
*EXPI1(2*P1*( =A Oe Wat 
*EXP1(2*PI*( 1 


Input data and output from the simulation of CH 


coupled DEPT experiment. 
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SAMPLE INPUT FOR DEPT EXPT ON CH SYSTEM 


COMM DEPT SEQUENCE FOR C-H ; THETA PULSE ANGLE = 
SPINI1S1 

SEs PULSE 45.0 
PULSSA+X 

PREG 

PULSSB-X 

PULSIA+X 

PRECT1 

SUBST1=1/( 2u) 
BUIESS Gay 

PULSIB-X 

PRECD 1 

SUBSD1=1/( 2u) 
MEASIXYO 

STOP 


OUTPUT FOR DEPT EXPT ON CH SYSTEM FOR VARIOUS THETA VALUES 


Carbon-13 Signals During Decoupled Acquisition 
Theta Pulse = 45 Degrees 


45 


SIGNAL = +( 70.711, -3.536)*EXP1(2*P] = ( =e Om Vil )*T1 
*EXPI (2*PI *( 1.00 VI )*D1 

Theta Pulse = 90 Degrees 

SIGNAL = +( 100.000, -0.000)*EXPI1(2*P]*( =i} (OO) Wat )*T1 
*EXP] (2*PI * ( 1.00 V1 )*D1 

Theta Pulse = 135 Degrees 

SIGNAL = ch AO PN. 3.536)*EXPI (2*PI *( =| (Oe) Yt )*T1 
*EXPI] (2*PI * | 1. WO Wit )*D1 

Theta Pulse = 180 Deqrees 

SIGNAL = + ( OC 5.000) *EXP1 (2*PI*( a0 Oma )*T 4 
*EXP] (2*PI * ( TOO Val )*D1 


Migure) 7. 


Simulation of the 6-pulse angle dependence in 
CH DEPT experiment. Carbon-13 signals 
COrrespondaing too = =45, 37 90.7. 155 ssance tom 
are given. 
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as *Leported by Doddre ll et val. (31). 

Simulations of the behaviour of the P2eHe and chs 
fragments with the DEPT sequence have been performed, and 
the output is given in Figure 8 and Figure 9 
respectively. The results are shown in pictorial form in 
Figures 10 and 11 where the transfer of proton 
magnetization to the multiple quantum coherences and on to 
carbon magnetization is shown for CH» and CH3 systems 
uSing the same scheme used in Figure 5 for the CH 
system. For clarity, the Boltzmann carbon magnetization 
components are not shown in Figures 10 and 11, since these 
may be removed by an initial carbon 90° pulse. The 
polarization transfer in all CH, cases follows the same 
route: (i) the proton magnetization is transferred to 
zero- and double-quantum coherences by application of a 
carbon 90° pulse when the two components of the proton 
magnetization are 180° out of phase (Figures 10B-C and 
11B-C); (ii) some of the zero- and doub le-quantum 
coherence is transferred to carbon magnetization by the 
proton 6(-y) pulse (Figures 10D-E and 11D-E); and (iii) 
the components of the transferred magnetization are 
allowed to precess until they are in phase and a net 
polarization transfer from the proton to carbon spin 


systems has taken place. In the CHy and CH3 cases, there 
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DEPT CALCULATION ON C-H2 SYSTEM ; THETA PULSE = 45 


Initial Density Matrix 


OM Hf, AW) Sse HO E00, Oe 

yy) W 2 2 Se sel tly isola) i) 4} 

Bi ( Bi) eS sel 2e5.00R C0a 

DMN 4a) Aa a+ (eo S00) ) } 

Om 4 Sy, BS) Bs Sty O., (8) )) 

DMM Ghe Gils = = 1)02 50.0) 0.0 

Oh Wk Se eal PX tS \0}0)\, 0.0 

ON feist Saal Sek Bal, 0.0 

Density Matrix After 90(Proton) Pulse 

BW aba AO Be PL 1S} 0)(0),, 0.0 

OM) i, Si S se ORO ane Olen iaial 

OW ( 2, Pie & =2 7,5 0/0 0.0 

om ( 2, 24) B ORG. —s OF /alntn 

DIN (Ch eh) ce PE RSI O0)|, soy} 

OM iy Sy SB set OG O71 

DM 454) = Oe SOO 0.0 ) 

OY 4, GS) Be O.@ » svn! 

DMS Sian 255007 OROn es) 

Ol WS) eS P25 S0}(0) 0.0 

pw ( 9, Ph a ss PE s)(0)(0), OO } 

Ot teh, te} se =A), S00) ORO) 

Density Matrix After 180(Proton)-90(Carbon) Pulses 

DM i, 2 ss 0.0 Ph KONO) 

Om (( a, 2h) Ss ORO 7 Oe tal SEXP G2 Pil*( i OOM VIS )*D1 ) 
DM AL, hi) Se oP OPO Ona BXPT (25 Pix San OO VS )*D1 ) 
DM Si, Zl) = ee On0_ | 2.500) 

DM 6, ta} & sri QO 5 S705 METS (Bea ea Sal stole) WS )*D1 ) 
i «4 Sy) s3t OO HO AEM TIE i 100) MS )*D1 ) 
i | Sy Gi & a On Oma, 2.500 ; 

DM @,, i) Se Oi . 25500 


DM 1, 2 2 se sy), Otol . 1.250)*EXPI (2*PI] x | =H (lO, Wat )*D1 
DM i Ch & dal 0.000, 1.768)*EXPI(2*PI* | =) OO) Wine )*D1 ) 
Ow of, We) Ss 50.000, 1.250) *EXPI(2*P] x ( =} (OO), Wa )*xD1 ) 
DM A Si) Ste 0.000, -1.768)*EXPI(2*P1*( Wa hOley iat )*D4 
OM 2, 2a)) Be 50.000, -1.250)*EXP1(2*PI*( GO) VI )*D4 
ol i, AE) Rail GOO Soe, -0.000)*EXPI(2*PI*( = OOM Vil )*D4 
OM) ( S ) Boo -0.000, -1.768)*EXPI(2*PI*( Saf {Olay Wi )*D1 ) 
DM Al (i) 2 (0), (OON0), 1.768) *EXPI(2*P1*( WOO Ware )*D1 ) 
om Nt, GH) SB aa Sis) Xeltee);, 1,250)*EXPI (2*P1*( =1,00 VI )*D1 
OM ( 7, Bi) ae s 0.0 , -2.500)*EXP1(2*P1*( =, 00 VI )*D1 
Carbon-13 Signals During Decoupled Acquisition 
SIGNAL = +( 199.999, -5,000)*EXPI(2*PI*( 1.00 VI )*T2 ) 
*EXPI (2*P 1 * Hi OOmNT )*D1 ) 
*EXPI(2*P I * ( 1.00 VI )*D2 ) 


Figure 8. Density matrix and output from the simulation 
of CH» DEPT experiment. 
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Figure 9. 


EPT CALCULATION ON C-H3 SYSTEM ; THETA PULSE © 45 
nities) Density Matrix 

DM il) eae 1525500 OF:0) >) 
DM ( 2, 2) # +( 147.500 0.0 ) 
DOMES 3) =a 52.500 0.0 ) 
Dm ( 4, 4) = +( 47.500, 0.0 ) 
OM ( 5, 5) © +( -47.500 00>) 
OM ( 6, 6) & +( -52.500, 0:10) )) 
OM ( 7, 7) & #(-147.500 0.0 ) 
DM ( 6, B) = +(-152.500, 0.0 ) 
OM ( 9, 9) & +( 52.500 0.0 ) 
OM (10,10) = +( 47.500 0.0 ) 
DM (11,11) & +( -47.500 0.0 ) 
DM (12,12) © -52.500 0.0 ) 
OM (13,13) = 52.500 0.0 ) 
DM (14,14) © +( 47.500 O70) 
OM (15,15) © +( -47.500 ORC) 
DM (16,16) © +( -52.500 0.0 ) 


Density Matrix After 90(Proton) Pulse 


DM ( 1, 1) & 4( 2.500, OF 0m) 

DM le ey es cot 0.0 , 66.602) 

DMI 2s (2) Sano re 2500), CxO s) 

DM ( 2, 4) & +( 0.0 , 86.602) 

DM SS) = a 2.500, ety) 

DMS See +i 0.0 , 100.000) 

OM ( 4, 4) & +( -2.500, OO) 

OM ( 4, 6) = +( 0.0 , 100.000) 

DM 5 Se + 2.500, 00) 7) 

DM ( 5, 7) = +( 0.0 86.602) 

DM ( 6, 6) = +( -2.500, (he) 

OW, 6, B) = / 0.0 , 66.602) 

OM tha, PN ee soll 2.500, 0.0 

pM 8. 8) = +( ~-2'2500 OOF 

OM, 9, 9) = +( 2.500, 0.0 

OM ( 9,11) & +( 0.0 -50.000) 

DMO nsO)me n+ (eer DO OF 0.0 

DM (10,12) = +( On0 ey, -50000) 

Iara) G call 2.500, ORO) 

DM (12,12) = +( -2./500), 0.0 

DM (13,13) = +{ 2.500, 0.0 ) 

DM 0135 tS); esa 0,0, , 50.000) 

DM (14,14) = +( -2.500, ORO) 

DM (14,16) = +( 0.0 , 50.000) 

DM (15,15) = +( 2.500, 0.0 

DM (16,16) = +( -2.500, 0.0 ) 

ensity Matrix After Theta(Proton!-180/Carbon) Pul 
DM tle, Co) ce oH 0.684, 53.033)*EXP1(2°P19( 
DMI eran earl 1.531, 30.618)*EXP1(2*P1*( 
OW ls Gi to aa 1.531, -30.619)*EXP](29PI=( 
DM ( 1, 8) & 4( 0.684, -53.033)*EXP1(2=Pl=| 
DME 2 ecu eies( 1.531, -30.61B)*EXP1(2ePl=( 
i (A Sh) ha cal 1.531, 30.619)*EXPI](2*Ple( 
yl (C2 A ae atl 0.884, 53.033)*EXPI](2*Ple( 
OM 3, 4) & +( 0.884, -88.388)*EXP1(2=Pl*( 
DW 3, 6) = +( -0.8684, -88.388)*EXP1(2*Pl*( 
DM 3, 8) = +( -1.531, 30.619)eEXP](2eP19( 
OM ( 4, 5) = +( -0.884, 68.388)*EXPI](2*P1*#( 
DM ( 4, 7) © +( -1.531, -30.619)*EXP1(2eP]s( 
DM ( 5, 6) © +( -0.884, 68.388)*EXP1](2*P1=( 
OM yy Weis st 1.531, 30.618)*EXP1(2ePl*( 
Dm ( 6, 7) #& +( 1.531, -30.61B)*EXP1(2°Pl*( 
DM ( 7, 6) & +( -0.884, -53.033)*EXP1(2*P1*( 
OM ( 9,10) © +( +-1.768, -35.355)*EXP1(2*Ple( 
DM ( 9,12) ® +( 1.768, -35.355)*EXP1(2*Ple( 
DM (10,11) = +( 1.768, 35.355)*EXP1(2*Pl*( 
DM (11,12) & +( 1.768, 35. 355)*EXP](2ePle( 
DM (13,14) © +( -1,768, -35.355)*EXP](29Pl#( 
DM (13,16) © +( -1.768, 35.355)eEXPI(2ePle( 
DOM (14,15) © +( -1,768, -35.355)*eEXPI(2ePl( 
DM (15,16) & +( 1.768, 35.355)*EXP1(2*Ple( 

r | ignal rd i} A igti 

SIGNAL «® +( 424.252, -7,.071)eEXPI(2ePle( 


Density matrix and output from the simulation 


Of CH, DEPT experiment. 


SEXP] (2ePl]*( 
EXP1(2ePI9( 
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are several zero-quantum and double-quantum coherence 
components, and parts of these coherences are converted by 
the @ proton pulse into multiple-quantum coherences at 
Frequencress2VSi+ Vi (for both CH, and CH3), and at 3VS + 
VilleCLor CHa Jem these multip Je-quantum coherences -are: not 
shown in Figures 10 and 11 since they are not detected in 
the DEPT experiment. 


For the CH, system, the computed signal intensity 


(Figure 8) can be written as 


srgnal = (A—-15) exp{ 211i Vi T2) >, [4-8] 


in the case of decoupled acquisition. When the carbon 
Spectrum is obtained in the absence of the proton 


decoupler, the signal intensities are given by (Figure 12) 
Signal = (A/ 2—16/4) exp{2n1t (Vitd)T2) + 
(A-i8/2) exp{2nxi VI T2} + [4-9] 
(h/ 2-167 4) expt 274. CVI=3 i025 9, 
whieh corresponds to the characteristic 1:2:1) triuplee. 
The dependence of the transferred polarization signal on 
the duration of the 6 pulse has been studied (Figure 


13). The transterred polarization follows the general 


expression 
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INPUT FOR COUPLED DEPT EXPT ON C-H2 SYSTEM 


COMM COUPLED DEPT EXPT FOR C-H2 SYSTEM ; THETA = 45 
SPINI1S2 

SET (PUESG 45.0 
PULSSA+X 

PRECD 1 

SUBSD1=1/( 2u) 
PULSSB+X 

PULSIA+X 

PRECD1 

SUBSD1=1/( 2d) 
PUSS Cay 

PULSIB+X 

PRECD2 

SUBSD2=1/( 2u) 
PRECT2 

MEASIXYO 

STOP 


OUNEUINEORECOUPRREDEDERINEXRIRONECSH2ESYSIEM 


Carbon-13 Signals During Coupled Acquisition 


SIGNAL = +(_ 50.000), =1).250)i*EXPI (2*PI>( 1.00 VI HoWO gee jj 
*EXPI(2*P I x ( = OO) Wa )*D1 ) 

*EXPI(2*P I ( 1.00 V1 )*D2 } 

S(T MOOR OCO FS S250 0) EXP In 2 Pi i 1.00 VI )*T2 ) 

*EXPI(2*PI* ( SOO) Val )*D1 ) 

*EXPI (2*P 1 ( 1.00 VI )*D2 ) 

sel SOS OOO, =a ASO reo eae | 1oOO Wi Si sOO W VE) 

*EXPI(2*P I * ( =e OOM Vel )*D1 ) 

*EXPI (2*P I ( 1.00 VI )*D2 } 


Figure 12. Output £rom ‘the simulation of CH» coupled DEPT 
experiment. 
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OUTPUT FOR DEPT EXPT ON C-H2 SYSTEM FOR VARIOUS THETA VALUES 


Carbon-13 Signals During Decoupled Acquisition 
Theta Pulse = 45 Degrees 


SIGNAL = +( 200.000, -5.000)*EXPI(2*PI = ( 


*EXPI(2*P I] * ( 


Theta Pulse = 90 Degrees 


SIGNAL = FA 0.000, -0.000)*EXPI(2*PI1( 


FEXPH (2*P Ie 


Theta Pulse = 135 Degrees 


SIGNAL = +(-200.000, -5.000)*EXPI(2*PI*( 


*EXPI(2*PI*( 


Theta Pulse = 180 Degrees 


SIGNAL = + | O50) =c10. 000) *EXPT(2* Pt =| 
*EXPI(2*PI1* | 


Figure 13: 
DEPT experiment. 


Simulation of §6-pulse angle 


SOO VA )*T4 
Ayn OOM Vel )*D1 


=} SOLO) Wat )*T1 
1.00 VI )*D1 


ail OO MVE )*T 1 
1.00 VI )*D1 


Se OOe Vil )*T1 
COM Vil )*D1 
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Signal = 2A sin(20) exp{2ni VI T2}) , [A=-10] 


which agrees with the results of Doddrell et al. (31). 
Results for the CH3 system are shown in Figure 9. 
The simulations of the DEPT experiment on the CH3 system, 
with carbon signal acquisition in the absence of the 
proton decoupling shows the characteristic 1:3:3:1 quartet 
pattern. This result is shown in Figure 14. The 
simulations (Figure 15) for investigation of the 
e-dependence lof the CH, signal antensities agree with the 


(sin 6 + sin 39) dependence given by Doddrell et al. (31). 


b. Selective 2D Heteronuclear Shift Correlation 

Experiment 

The spectrum editing capabilities of the DEPT 
sequence have been utilized by Levitt et al. (54), and 
Bendall et al. (55), to produce selective two dimensional 
heteronuclear shift correlation maps. A selective shift 
correlation map displays peaks arising from the shift 
correlation of only a particular type of CH, group, thus 
making it possible to differentiate between signals from 
CH, CHy and CH3 fragments. The pulse sequences used in 
Both investigations (54,55) yield similar correlation 
maps, and the pulse sequence proposed by Levitt et al. 


(54) will be analyzed using the simulation program. 
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INPUT FOR COUPLED DEPT EXPT ON C-H3 SYSTEM 


COMM DEPT SIMULATION FOR C-H3 SYSTEM 
SPINI1S3 

SET PULSC 45.0 
PULSSA+X 

PRECD1 

SUBSD1=1/( 2u) 
PULSSB+X 

PULSIA+X 

PRECD1 

SUBSD1=1/( 2u) 
PULSSC-Y 

PULSIB+X 

PRECD2 

SUBSD2=1/( 2u) 
PRECT2 

MEASIXYO 

STOP 


; THETA PULSE=45 


OUTPUT FOR COUPLED DEPT EXPT ON C-H3 SYSTEM 


Carbon-13 Signals During Coupled Acquisition 


SIGNAL = sil SOs WSs SO) Shey este Tose eae esl 1.00 V1 Wee UW yee) 
® Pea 2s Pale Sal OOMVa )*D1 ) 

*EXPI(2*P] = ( 1.00 V1 )*D2 ) 

+ (159,089), =2652) *EXP1(2*P1*( 1.00 VI Dos) od) era), 
*EXPI(2*P]*( OOM VA )*D1 ) 

*EXPI] (2*P I * | 1.00 VI y= D257) 

+ SOROSS)  =2 4652) SE AP IN 2 Pie) LOO Woo oe .'s@ wl yew | 
ar OPM 2.76) tea SH} ol. Wat EO 

*EXPI(2*PI] «| 1.00 VI )*D2 ) 

(SS nOCshe a OnOe4) SEXP lee Pil #i Toe) Wi Sj) eae | 
*EXPI (Qe PI «| Al (Ole) Wx )*D1 ) 

*EXPI(2*P] * ( 1.00 V1 )*D2 ) 


Figure 14. Output from the simulation of CH3 coupled DEPT 
experiment. 
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OUTPUT FOR DEPT EXPT ON C-H3 SYSTEM FOR VARIOUS THETA VALUES 
Carbon-13 Signals During Decoupled Acquisition 
Theta Pulse = 45 Degrees 


SIGNAL = +( 424.260, -7.071)*EXPI(2*PI*( ne. Wal Heal 
*EXPI (2*PI1*( je OONVAL )*D1 ) 


Theta Pulse = 90 Dgerees 


SIGNAL = + ( 0.000, -0.000)*EXP1(2*PI#( Solo) Wel )*T1 ) 
*EXPI(2*PI1*( 1.00 V1 )*D1 ) 


Theta pulse = 135 Degrees 


SIGNAL = +( 424.257, 7.071)*EXPI(2*PI*( -1.00 VI eat) 
EXP] (2eP1*( 1.00 V1 )*D1 ) 


Theta Pulse = 180 Degrees 


SIGNAL = = ( 0.000, 20.000)*EXPI(2*PI*( =p ooo) Ya y*T1 ) 
*EXPI (2*PI*( 100M VA )*D1 ) 


Figure 15. Simulation of 6-pulse angle dependence in CH3 
DEPT experiment. 
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The selective 2D shift correlation pulse sequence is 
shown in Figure 16. The original DEPT sequence (31) was 
modified to include a variable time period Tl, during 
which the proton magnetizations are frequency labelled. 
Thus the proton magnetization at the end of the T]) time 
period is modulated by proton resonance frequencies. 
Arter the introduction of this modulation, the DEPT 
sequence was used to transfer the modulated proton 
magnetization into carbon magnetization. The carbon 
Signals acquired in the presence of proton decoupling is 
them a £unction of two time variables Tl “and 172. A two 
dimensional data set is acquired by collecting the FIDs 
GEunctrons Of “T2) “for avtseries of values sor Tl. ‘The two 
dimensional data is then Fourier transformed twice, with 
respect to the time variables T2 and Tl, to produce the 
shift correlation spectrum. The 2D spectrum has proton 
frequencies along the Vy dimension (Fourier conjugate of 
time T1) and the carbon frequencies appear along the v5 
dimension (Fourier conjugate of time T2). 

AS mentioned in Section IV.2.a, the signals acquired 
after the DEPT pulse sequence depend on the 6-pulse 
duration, with the intensities proportional to sin 6, 
sin 26 and (sin @ + sin 36) for CH, CHy and CH3 systems 
respectively. The linear combinations of FID data at 


three different 9 values have been used to generate the 
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individual subspectra of the CH, systems. The 6 values 
used in the 2D DEPT work (55,56) were 45°, 90° and 135°. 
Anotner set of @ values 30°, 90° and 150° faci litates 
easier manipulation of the 2D data (57). 

inworGgeaeto illustrate the application of the 
Simulation system to the analysis of a selective 2D shift 
corre lation €xperiment, computations were carried out for 
CH, CHy and CH3 systems. The input pulse sequences and 
outputs are given in Figures 17-19. The CH case will be 
used to illustrate the mechanism of shift correlation. 

The calculations have been carried out for 9 = 45°. 
The density operator for the CH system is represented in 
the simple product basis, as in Equation [4-1]. The 
initial proton) 902 spulsescreates transverse proton 
magnetizations which precess during the evolution period 
Tl. This introduces a frequency labelling of the proton 
Magnetization. The carbon 180° pulse at the midpoint of 
the Tl period serves to eliminate the effects of carbon- 
proton coupling on the proton magnetizations. The density 
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INPUT FOR 2D SELECTIVE DEPT EXPT ON C-H SYSTEM 


ee 


COMM 2D SELECTIVE DEPT EXPERIMENT ON C-H SYSTEM 


SPINI1S1 
SETMPUESE 45.0 
PULSSA+X 

PRECT1 

PULSIB+X 

PRECT1 


WRDM 


PRECD1 
SUBSD1=1/( 2d) 
PULSSB+X 
PULSIA+X 
PRECD 1 
SUBSD1=1/( 2u) 
PULSSC+Y 
PULSIB+X 
PRECD2 
SUBSD2=1/( 2u) 


DCPLS 


T2 


MEASIXYO 
STOP 


OUR URE BOR 2 DS EMEC IHR D ERiteX Pi ONSGS Hm SViomeEM 


Density Matrix After 11 Evolution Period 


DN, 
DM 
DN 
DN 
DM 
DM 


( 
( 
( 
( 
( 
( 


’ 
’ 
’ 


) 
) 
) 
) 
) 
) 


Te | 


aril Pa OLD, OOM) 
+ | 0.0 , 50.000)*EXPI(2*PI*( 2.00 VS 


+ ( Pata), Os® |; 
+ ( 0.0 , 50.000)*EXPI(2*PI*( 2.00 VS 
) 


+ 229.500.) 056 
Hi) e500 m0 Om) 


Carbon-13 Signals During Decoupled Acquisition 
SIGNAL = 


EUOURe TL 7. 


+ ( 0.000, 3.536)*EXPI (2*PI*( 
*EXPI(2*PI* ( 


*EXPI(2*P1*( 


+( -35.355, -0.000)*EXPI(2*PI*( 2.00 VS 
*EXPI(2"P] » ( 
EXP] (2ePI*( 
*EXPI] (2*P] »( 


+( -35.355, -0.000)*EXPI(2*PI*( -2.00 VS 
*EXPI (2*PI *( 
*EXPI (2*PI1*( 
*EXPI (2*PI*( 


Input. Gata and olitpue from the simulation cf 
CH selective shift correlation experiment. 
The 6 pulse (pulse C) has been set to 45°. 
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INPUT FOR 2D SELECTIVE DEPT EXPT ON C-H2 SYSTEM 


COMM 
SPINI1S2 
Sail TUsSe 
PULSSA+X 
PRECT 1 
PULSIB+X 
PRECT 1 
WRON, 
PRECD1 
SUBSD1=1/( 2u) 
PULSSB+X 
PULSIA+X 
PRECD1 
SUBSD1=1/( 2u) 
PUESSEAY 
PULSIB+X 
PRECD2 
SUBSD2=1/( 2) 
DCPLS 2 
MEASIXYO 

STOP 


45.0 


CUTLPWiEORMZDESEDECT IVES DERIMEXPIAONEC=H2ES SIEM 


Density Matrix After T1 Evolution Period 


( 
( 
( 
( 
( 
DM ( 
( 
( 
( 
( 
( 
( 


Carbon-13 Signals During Decoupled Acquisition 


SIGNAL 


Figure 18. 


> 


OnrNIMDMNH PWWNHNM — = 


= Son 5008 C207) 

3) = +{ OO Old TA EXPT (2 Pulse 
Dy ee sah 225.00), On Oe) 

4) = +( O20) 700711 )*EXPT(2=PI=( 
Si) Beh 2 SOO), 0.0 

5) 2 sx 0.0 70.711)*EXPI(2*PI*( 
Li ES all 2.500, ORC) 

6) = +( OS ON OL IPE NP Ta2% Palen 
Bh) Saal See Eyl): 0.0 ) 

6) = +( ZOO; 0.0 ) 

i) See Asse). ORCaay 

3h) SS Sel| 2.500, GRC) 


5.000)*EXPI(2*PI*( 
ede \( el DES 
*EXPI (2*P1*( 


= + ( 


0.000, 


-0.000)*EXPI(2*PI*( 
*EXPI (2*PI *( 
*EXPI] (2*PI*( 
*EXPI(2*PI «| 


+(-100.000, 


-0.000)*EXPI(2*PI*( 
*EXPI(2*PI*( 
*EXPI (2*P1 *( 
*EXPI(2*P1*( 


+(-100.000, 


imputsdata and outputetromathessimubeationtor 
CHs, selective shift correlation experiment. 
(pulse C) has been set to 45°. 


The @ pulse 


2.00 VS 
2.00 VS 
2.00 VS 
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2D SELECTIVE DEPT EXPERIMENT ON C-H2 SYSTEM 
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JNPUT FOR 2D SELECTIVE DEPT EXPT ON C-H3 SYSTEM 


COMM 2D SELECTIVE DEPT EXPERIMENT ON C-H3 SYSTEM 
SPIN1 183 

SET PULSC 45.0 
PULSSA+X 
PRECT1 
PULSIB+X 
PRECT) 

WRDM 

PRECD1 
SUBSD1=1/( 2d) 
PULSSB+x 
PULSIA+X 
PRECD1 
SUBSD1=£1/( 2u) 
PULSSC+Y 
PULSIB+X 
PRECD2 
SUBSD2=1/( 2u) 
DCPLS T2 
MELSIXYO 

STOP 


QUTPUT FOR 2D SELECTIVE EXPT ON C-H3 SYSTEM 


ensity Matrix After Ti Evoluti Peri 


DN ieee = 25500 OF Cie) 

DM ( 1, 3) = +( 0.0 , 66.602)*EXPI(2=PIl*( 2.00 VS jeT1 ) 
DM (252 =) + 25500) OOS) 

Dm ( 2, 4) & +( 0.0 , 86.602)*EXPI(2=PI*( 2.00 VS tain gl 
DVM Sees!) cease 255005 0.004 

Dione sino) eat 0.0 , 100.000)*EXPI(2*Pl*( 2.00 VS yeti ) 
om { 4, 4) & «( 2.500, (et), ) 

DM ( 4, 6) & + 0.0 , 100.000)*EXPI(2=PI*( 2.00 VS jeT1 ) 
OM ( (5, 5) # +( ~=2-500; (ora). 

ol AC. WO) 3 sal 0.0 , 86.602)*EXP1(2=PIl*( 2.00 VS yeT1 ) 
DM ( 6, 6) = + 2500), ORO) 

DM (16, 6) = + 0.0 , 66.602)*EXPI(2ePl*( 2.00 VS )eT1 ) 
DM 74 7) ee 27/500; 0.0 ) 
DM ( 6, B) © + 2.500, ORO me) 
DM ( 9, 9) & + =2. 500) ORCI) 

DM ( 9,11) & + 0.0 , -S0.000)*EXPI(2*Pl*( 2.00 VS yeT1 ) 
DM 110,10) © + 2.500, OF 0) 

DM (10,12) © + 0.0 , -50.000)*EXPI(2=°Pl*( 2.00 VS )eT1 ) 
(oy (ESIGN alta) eS 2 -2.500, 0.0 ) 
DM (12,12) = +( 20500) 0.0 ) 
OM (13,13) © +( -2.500, OO) 

OM (13,15) = +( 0.0 , 50.000)*EXP1(2=Pl*( 2.00 VS yeT1 ) 
DM (14,14) = +( 2.500, O01) 

DM (14,16) © +( 0.0 , 50.000)*EXPI(2=Pl*( 2.00 VS )eT1 ) 
DM (15,15) © +( -2.500, O70) } 
DM (16,16) = +( 2.500, OO) 

rbon-13 Signal r4 i A isiti 

SIGNAL = + OT 2 ei olig O.000)*EXP1(2=Pl*( 2.00 VS )eT1 ) 

*EXP] (29Pl#( 1.00 VI Meas )} 

@EXP] (2ePle( =1,00 V3 )*D1 ) 

PEXP](2°P]=( 1.00 V1 )*D2 ) 

+(-212.131, O.000)*EXPI(2eP]e/ -2.00 VS )=Ti) 

PEXPL(2eP lo ( 1,00 VJ )eT2 ) 

EXP] (2=P]e( =1,00 VI )eD1 ) 

PEXP] (2ePle( 1.00 VI yeD2 ) 

+( 0.000, 7.071)*EXP1(2*Ple( 1.00 VI )eT2 ) 

EXP] (2ePle( -1,00 V1 )eD1 ) 

EXP] (2*Ple( 1.00 V1 )eD2 ) 


Figure 19. Input data and output from the simulation of 
Cig selective shift correlation experiment. 
The @ pulse (pulse C) has been set to 45°. 
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noted that the density matrix elements P73 and p94 which 
correspond to proton transverse magnetizations are 
frequency labelled’ in Tl. “Further, the evolution of the 
density matrix is similar to that in the DEPT pulse 
sequence described in Section IV.2.a. The proton-shift 
labelled magnetizations are transferred to the multiple- 
quantum coherence levels by the carbon 90° pulse. They 
are then converted to observable carbon magnetization by 
the proton 9 pulse. The chemical shift offsets during the 
precession periods Dl, D2, are refocussed to yield carbon 
signals free of intensity and phase distortions. The 
carbon signals observed under proton noise decoupling is 


given by (Figure 17) 


Signal = {A/v2) cos{2n VS Tl} expt2ni VI T2} + 


(S/72) expi2n. Vi T2) 7; [4-12] 
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after recognizing that Dl = D2. Double Fourier 
transformation of this signal with respect to Tl and T2 
yields the 2D spectrum. 

It should be noted that in Equation [4-12], there is 
no quadrature detection in the Tl time period and hence in 
its frequency counterpart, the v,; dimension. Lack of 
quadrature detection in Tl, dictates that the proton 
carrier frequency be placed at one end of the proton 
Spectral width for distinction of plus and minus 
EVequencies. = This 1s inefficient. in that the v, Spectral 
width is doubled, leading to the situation where twice the 
number of Tl values are required to achieve the same 
resolution achieved using a quadrature detection method. 

A phase cycling sequence for introducing quadrature 
detection in Tl has been worked out using the simulation 
program and is given in Figure 20. The carbon signal at 
the end of this 4-phase cycle sequence is given by 


(Figure 20) 


evonaimear (N/ 272) Expl 2n1, VS Il) expi2nl Vi wi) =. [4-13] 


Another interesting feature of the signal in Equation 
[4-13] is the absence of carbon Boltzmann signals. The 
carbon Boltzmann signals are nulled by cycling the first 


90° carbon pulse through all four phases of the rotating 
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INPUT FOR 2D SELECTIVE DEPT EXPT WITH QUADRATURE DETECTION 


COMM QUADRATURE DETECTION IN T1 DOMAIN 
SPINI1S1 

SET SRUESC 45.0 

PHASECYCLE 4 


PULSSA+XSA+XSA+XSA+X 
PRECT1 
PULSIB+XIB+XIB+X1B+X 
PRECT 1 

PRECD 1 

SUBSD1=1/( 2u) 
PULSSB+XSB+XSB+XSB+xX 
PULSIA+XIA+YIA-XIA-Y 
PRECD 1 


SUBSD1=1/( 2u) 
BUESSC+YSCSNSG=VSEXX 
PULSIB+XIB+X1B+X1B+X 
PRECD2 

SUBSD2=1/( 2d) 

DCPLS T2 
MEASIXYOIXYOIXYOIXYO 
PHASECYCLE 4 

STOP 


OUTPUT FOR 2D DEPT EXPT WITH QUADRATURE DETECTION IN 11 
Carbon-13 Signals (Decoupled) After 4-Phasecycles 


SIGNAL = +(-141.421, -0.000)*EXPI(2*PI*( 2.00 VS Very | 
*EXPI (2*PI*( 1.00 VI )*T2 ) 
*EXPI(2*P1*( Sale COVA } =D) 
*EXPI(2*PI1*( 1.00 VI )*D2 ) 


Figure 20. Input data and output from the 4-phase cycle 
CH selective shift correlation experiment. 


The signal exhibits quadrature phase in Tl 
domain. 
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frame. The phase of the second 90° proton pulse is 
rotated in the same direction but with a 90° relative 
phase shift from the carbon pulse. This ensures that only 
the transferred polarization signals are detected by the 
carbon receiver. This phase cycling also effects 


quadrature detection in v, (Figure 20). 


Cc. INADEQUATE Experiment 

The INADEQUATE experiment has been developed to 
observe natural abundance l3c-1l3c couplings without 
interference from the much more intense resonance due to 
the corresponding 13¢e-12¢ fragments (58). This sequence 
takes advantage of the phase properties of multiple- 
quantum coherences. Variations of the basic INADEQUATE 
sequence have been employed in the observation of 13¢-l3¢ 
couplings with enhanced sensitivity (62); in a two 
dimensional version for easier interpretation of spectra 
(63); for observation of long range 13¢-13¢ couplings 
(64) and form tracing, out the carbon skeleton inva 
mMovecule (59). Lheamain ditticudty in 1dentifying (the 
weak 13c-13¢c satellite lines in normal carbon spectra 
stems from the presence of many other weak signals due to 
Spinning side bands, incomplete proton decoupling and 
samp le impurities (58). The basic INADEQUATE pulse 


sequence is shown in Figure 21 and contains the minimum 
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90, 180, 90, 90) 


| 1 D1 acquisition T2 
Read Pulse (p ) 


+X 


Receiver 


INADEQUATE PULSE INPUT FOR C-C SYSTEM 


HOMO 

SPENT IS 
WRDM 
PULSSA+X 
WROM 

PRECD1 
PULSSB+X 
PRECD 1 
SUBSD1=1/( 4J) 
WRDM 
PULSSA+X 
WROM 

PRECT 1 

WRDM 
PHASECYCLE 2 
PULSSA+XSA+Y 
PRECT2 

WRDM 
MEASSXYOSXY 1 
PHASECYCLE 2 
STOP 


Pigure Zl. Pulse sequence and input data for simulation 
of the INADEQUATE experiment on 13¢_13¢ 
fragment. 
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phase cycling required to remove the strong resonances 

from molecules which contain only a single 13c nucleus. 
More complicated phase cycles are often used to correct 
for pulse imperfections and other experimental problems 
(62). 

In this section, the simulation of the INADEQUATE 
sequence applied to a homonuclear I1S1 system (a 13¢_13¢ 
fragment), to the S spins of a heteronuclear I1S1 system 
oe 


(a 12, fragment), and to a homonuclear I1S2 system (a 


13¢_13¢_13¢ fragment in an isotopically enriched molecule) 


are described. 

The input to the simulation program to perform the 
INADEQUATE sequence on the I1Sl1 homonuclear system is 
shown in Figure 21. The output from the simulation 
program is given in Figure 22. The initial density matrix 


for the homonuclear I1S1 system (Figure 22) can be written 


as 

2h @) 0] O 

O ) 0 O 
ClO )ee= Cle) 2)) : (A-14 ] 
a O (0) O O 

@) @) @) —2A 


The initial 90°(+x) pulse puts all of the magnetization 
into the transverse plane and the four components pj5, 


P34- P13 aNd pp, are proportional to the complex 
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Figure 22. 


INADEQUATE SEQUENCE FOR C-C SYSTEM 
Initial Density Matrix 


omon on 


+( 100.000, ORC) 
+( 00) 0.0 ) 
+( OO O70) 
+(-100.000, 070) ) 


Density Matrix After First 90 Pulse 


Densit 


DM 
DM 
DM 
DM 


( 
( 
( 
( 


( 
( 


( 
( 


1, 
ihe 
2 
35 


1, 
ily; 
2m 
35 


2) 
) 
) 
) 


Matrix After Second 1 


2) 
3) 
4) 
4) 


oon 


cA 1330), 
call Ee 
= (S50) 
er\( S00 


0 
0. 
+( 0. 
0 


0 50.000) 
0 6,6 | 80). 000) 
Oi 0-000) 
0 , 50.000) 
(4JU) Precession Period 
000, -0.000) 
000, -0.000) 
000, -0.000) 
000, -0.000) 


Density Matrix After Second 90 Pulse 


DM 


( 


1, 


4) 


se 0. 


OF 100710003) 


Density Matrix After Precession Period 11 


DM 


OM 
DM 


DM 
DM 


DM 
OW 


OM 


( 


( 
( 


( 


( 


DM ( 


DM 
DW 


DM 
DM 


DM 
DM 


OM 
DM 


DM 


OM 
DM 


DOM 


{ 
{ 


{ 
{ 


( 


( 


1, 


1, 


4) 


1) 
2) 


3) 
4) 


2) 
3) 


4) 


1) 
2) 


we 


Cm 


= 0. 


+( 50 
a= OO 
call Ay), 
+{/ =50 
+{ =50 
(0) 
+( 50 
SAN SAO) 
+( 50 
pall cal0}e 


molt) ASO) 5 
+ | le 
+/ 0. 
pall) “aaca0)y 
th 50 
+( 50 
+( 0 
well Aa), 
+( 0: 
G3 tall osioh3 


Carbon-13 Signals During Acquisition 
SIGNAL 


+(-100. 


+( 100. 


+( 1100". 


+(-100. 


0 , 100.000)*EXP1(2=PJ#( 1.00 VS 1.00 VI )eT1) 
O.M. During Acquisition, After Third 90 Pulse (+x ; Phasecycle 1) 
000, 0.0 )* SIN(2*Pl*( 1.00 VS 1.00 VI )*T1) 
000, 0.0 )* COS(2*PI*( 1.00 VS 1.00 VI )eT1) 
*EXPI](2=P]*( 1.00 VI 0.50 uU)*T2) 
000, 0.0 )* COS(2*PI*( 1.00 VS 1.00 V1 )eT1) 
*EXPI(2*Pl]*( 1.00 VS 0.50 J)*T2) 
.000, 0.0 )* SIN(2*Pl*( 1.00 VS 1.00 VI )=T1) 
*EXP](2*Pl*( 1,00 VS 1.00 V1 )*T2) 
000, 0.0 )* SIN(2*Pl*( 1.00 VS 1.00 V1 )=T1) 
.000, 0.0 )* SIN(2*PI*i 1.00 VS 1.00 V1 )=T1) 
*EXP1(2*PI*( 1.00 VS -1.00 VI )*T2) 
000), 0.0 )* COS(2*P]*( 1.00 VS 1.00 VI )*T1) 
*EXPI](2*Pl*( 1.00 VS =O 50) ve te) 
.000, 0.0 )* SIN(2*PI*( 1,00 VS 1.00 VI )*T1) 
.000, 0.0 )* COS(2*PI*( 1.00 VS 1.00 V1 )eT1) 
*EXP1(2*PI] =| TOOF Vie = O050) Us) 23) 
000, 0.0 )* SIN(2=*Pl*( 1.00 VS 1.00 VI )eT1) 
D.M. During Acquisiton, After Third 90 Pulse (+y ; Phasecycle 2! 
000, 0.0 )* SIN(2*PI*({ 1.00 VS 1.00 VI )=T1) 
0 , 50.000i* COS'2=PI*( 1.00 VS 1.00 V1 )=T1) 
*EXPI](2*PI] =| 1.00 V1 0.50 J)*T2) 
0 , 50.000!* COSi2=PI*( 1.00 VS 1.00 VI )*T1) 
EXP] (2*PI*( 1.00 VS 0.50 JU)*T2) 
000, OO) )* SIN(2=P T= 17.100 VS OOF Vi )eT1) 
sEXP1](2=PI*({ 1.00 VS 1.00 VI )*T2) 
.000, 0.0 )d* SIN(2*Pl*( 1.00 VS 1.00 V1 )*T1) 
.000, 0.0 J= SIN(2*PI*( 1.00 VS 1.00 V1 )=T1) 
*EXP](2=PI*( 1.00 VS -1.00 VI )=T2 
.0 , -50.000)* COS(2*PI*( 1.00 VS 1.00 V1 )*T1) 
*EXP1(2*P1*( 1.00 VS -0.50 J)*T2 
000, 0.0 )* SIN(2*PI*( 1.00 VS 1.00 VI jet 
0 , -50.000)* COS(2*Pl*( 1.00 VS 1.00 V1 )eT1) 
EXP] (2*Pl*( 100 Vi =0.50 Jhet2 
000, 0.0 )* SIN(2*PI*( 1.00 VS 1.00 V1 y=T1) 
000, 0.0 )* COS(2*Pl*( 1.00 VS 1.00 V1 )=eT1 
SEXP] (2*PI]*( 1.00 VI 0.50 J)=T2 
000, 0.0 )* COS(2*PIl*{ 1.00 VS 1.00 ¥1 yeT1) 
*EXP] (2*PI*( 1.00 VI -0.50 J)#T2 
000, 0.0 J* COS(2«P}*{ 1.00 ¥S 1.00 VI )eT1 
PEXPI(2=Pi*( 1.00 VS -0.50 J)#T2 
000, 0.0 )* COS(2*PI*( 1.00 VS 1.00 V1 )*T1) 
*EXP](2*Pl*( 1.00 VS 0.50 J)*T2 
INADEQUATE 


Output data from simulation of 
experiment on 13¢_ 


8) 


C fragment. 
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magnetizations associated with the resonances at VI + J/2, 
Vir= J/2, 2VS +uJ3/2 .<andivSi-"J/2 réspectivelys ‘During the 
first precession period, these magnetization components 
dephase and if the duration of the precession is exactly 
a7 4s); phase differencessoh190° between itheryic— J/2 and 
VI + J/2 magnetization vectors and between the vS - J/2 
and VS.+ J/2 vectors are built up.” The 180° pulse is 
applied to cancel the phase shifts due to resonance 
offsets which build up during the first precession period 
by refocussing them during the second precession period. 
Thus, at the end of the second precession, the 
magnetizations associated with the VI - J/2 and VS - J/2 
resonances will lie along the -x axis of the rotating 
frame and the magnetizations associated with the VI + J/2 
and vS + J/2 resonances will lie along the +x axes. The 
90°(+x) pulse causes all of the magnetizations to be 
transformed into the double-quantum coherence pj, which is 
modulated at frequency VI + VS. During the Tl delay, the 
double quantum coherence acquires a phase shift of 

a (VEIFVSNEl. osthe density matrix abethis polnt;  ismgiven 


pyn( Bigure=22) 


0) 0 0 a 

0 0 0 0 
BU20IeT 1) = (1/72) ’ [4-15] 
os 0 0 0 0 
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where a = 21A exp{2ni(VS+VI)T1}. The 90°(6) pulse 
converts most of the double-quantum coherence back to 
observable magnetization. For 6 = +x, the VI + J/2 anda 
VS + J/2 magnetizations lie along the -x direction 
(ignoring the small phase shifts 2n(VS+VI)Tl) and the 
VI - J/2 and vS - J/2 magnetizations lie along the +x 


direction. For s$ = +y, the Vi +.3/2 and VS 4.3/2 


magnetizations lie along the +y direction and the VI - J/2 


and VS - J/2 magnetizations lie along the -y direction. 


The phase of the detector is chosen to measure M, + iM, 


£OF t.= +x and “My cai 


intensity patterns at VI + J/2 and at VS + J/2 in the 


for ¢ = -y, which yields 1:-1l 


Fourier transformed spectrum. 

The input to the simulation program for the 
INADEQUATE sequence on the S spins of an I1S1 system is 
Similar to that for the homonuclear I1S1 system given in 
Figure 21, but the HOMO and SUBS commands are deleted. 
The simulation results for this case are given in Figure 
23) in, this case; all pulses ane applied only to thers 
spin so that first 90°(+x) pulse simply generates 


transverse magnetization components at frequencies 


vs + J/2 and vS - J/2, associated with the density matrix 
elements p}3 and 274 respectively. Here J represents the 


coupling of the S spin (13c) to a heteronucleus (1H), not 
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INADEQUATE SEQUENCE ON HETERO NUCLEAR 1 S SYSTEM 
Initial Density Matrix 


i ty aD) Sse yt, syeie) O@ } 
Wi te ZA S salt 27 also) Oo } 
oi | Oy Si) es Sel Sole Do 
DM 4 4y =H =5 22,500 0.0 ) 


Density Matrix After First 90 Pulse 


DN yee) eae PL 00), OR Oe) 
DMT ce =a 0.0 50.000) 
DN Nar eee eae Nae O10 Da ) 
yn A, Gil eS api 0.0 50.000) 
aM Gy SHS sal 200), OR0m) 
Bil MC 2th Gh Ee cal Pel TiO0) st) 


Be te diese ca 2.500, OF Oi) 
ON (ah, Sp Se Sk Or, =e (NOLO) 
OM ( 2. Bi sal 2, S00), j. )) 
om t 2, A eB eel OG 5 SS), Wolo) 
OW Ml Sy, SO) S sei 25 500%, 0.0 

DIG aire (ee Oe G0  } 


Density Matrix After Second 90 Pulse 


experiment on 


3e-12. fragment. 


Bey NW oh, TU) ee sell Sel), 0.0 } 

ON ( By 2a) ea cal Ai SS0). 0.0 

DM 6, 8) 2 ol =47 S500, 0.0 

DM ii AN) ee catt ata) Ia(00),. 0.0 

Density Matrix After Precession Period 11 

DM ty a) €& of 62.500) 0.0 

DM 2 BO) Eis 7477 Sao OF0" *) 

DM S, oi) 2 sl 47 S00) Omo 

DM A 7) cell lay) Si(0)0) 0 }) 

D.M. During Acquisition, After 90(+X) Pulse (Phasecycle 1) 

DM Hig A pal PL M00, 0.0 } 

DM ty Oy tee 0.0 , 50.000)*EXP1(2*PI*( 1.00 VS 0.50 J)*T2) 
Bi ( 2, 2) SB set S200), O.@ 4 

Oy tl By Ay ssa 0.0 50.000)*EXP1(2*PI*( 1.00 VS -0.50 J)*T2) 
Die We Sh 3S)! eck | Ph NON), 0.0 

By) Zk UN) eS SEC et iO) 0.0 

D.M. During Acquisition, After 90(+Y) Pulse (Phasecycle 2) 
yi hy ad) ee SE Ps S08), QO. }} 
OM ( i, Ch) S cat Seie) OW, 0.0 )*EXPI(2*PI*( 1.00 VS 0.50 J)*T2) 
Bi By By) B sai Saree). O)( )) 
DM i 2h 4 ae 5 ORO OOF 0.0 )*EXPI(2*PI*( 1.00 VS -0.50 J)*T2) 
oy tt Si, ey Se Sell 2.500 0 3 
eM 2h 2) Be mel Se IS i0N0),, On(0)  } 

Carbon-13 Signals During Acquisition 

RESULTANT SIGNAL IS ZERO 

Eiqure 23.) Output data from simulation of INADEQUATE 


a l3c-l3¢ coupling. The 180° pulse causes the phase 
shifts which build up on P13 and py, during the first 
precession period to be refocussed in the second period 
and the magnetization vectors for the two components both 
lie along the -y axis. Hence the second 90°(+x) pulse 
rotates these components to the +z axis and the density 
matrix is identical to its initial thermal equilibrium 
values “Lt does not change during the fT] intervals The 
90°(6) pulse converts the carbon z magnetization into 
components along the +y (for » = +x) or -x (for 6 = +y) 
directions. With detector phases set to measure Mote iM 


Y 


fOr 6 = +X and;=M., + 1M) for 4 = +y, at is clear that ‘the 


y 


sum of the two acquisitions will be zero as the output in 


x 


Figure 23 shows. 

It is interesting to consider the simulation of the 
INADEQUATE sequence on an I11S2 homonuclear system in order 
to determine the behaviour of a molecule containing three 
coup led 13¢ nuclei with this pulse sequence. This example 
is relevant in 13c studies on isotopically enriched 
molecules. The output of the simulation for this case is 
given in Figure 24. The components of the carbon FID with 
frequencies VI —- J, VI, VI + J, VS = J/2 and VS + J/2 have 
different phases so the Fourier transformed spectrum will 
Consist Of both dispersive and absorptive components. If 


we assume that the delay Tl is very short (in practice ~10 
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INADEQUATE SEQUENCE ON THREE COUPLED CARBON (C-C-C) SYSTEM 


Carbon-13 Signals During Acquisition (from 2 phasecycles) 


SIGNAL 


Figure 24. 


+( 


Tal 


Output data from simulation of the INADEQUATE 


experiment on 3¢-13¢_-13¢ fragment. 
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uS), the phase factors involving double quantum 
brequencies VS" + vi +e /2eand 2VS  Jycan be et to 

unity. Then, if the spectrum is phased so that the lines 
aiaavs = Vandy VSi+ J “arelabsorptions with ire lative 
bntensities"2:=2, the line at VI will be purely dispersive 
and have relative intensity 1; while the lines at VI - J 
and VI + J will have absorptive components of intensity 
l:-1, and dispersive components of magnitude -1/2:-1/2. 

It should be noted that the multiple-quantum signals arise 
from even-order coherence, i.e. in the case of 1182 
system, from zero- and double-quantum coherences. The 
INADEQUATE signals for a molecule containing three 13¢ 
nuclei are rather different from those for molecules with 
only two 13¢ nuclei in that the phases of the various 


lines are Significantly different. 


d. Unexpected Effects of Spin Decoupling 

Recently Levitt et al. (53) have demonstrated that 
unanticipated effects can be observed in certain 
experiments employing spin decoupling. They have studied 
the effect of heteronuclear decoupling applied to the §S 
spins, on the resonance of another nuclear species I, in 
IS, Spin systems which are in arbitrary states at the 
Start of the decoupling period. It was concluded that, 


even under strong decoupling the two spin systems, I and 
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S, do not behave independently. The multiplet components 
of I spin magnetization are affected by the rf 
perturbation to the S spin when the I multiplet components 
differ in magnitude or phase at the beginning of the 
decoupling period. 

An analysis of the time evolution of the density 
operator o(t) during the decoupling period, was carried 
Out by hevitt etal, (53), "and has shown that the 
necessary conditions for the S spin decoupling field to 


affect the I spin expectation value <Q> are: 


(etier, 600) lar oO 
[4-16] 
(2) gitar Odets0 


where He is the S spin Hamiltonian during the decoupling 
period and 0(0) is the density operator at the start of 
the decoupling period. The above conditions have been 
used to explain the effect of S spin decoupling on the I 
spin multiplet components (53). Consider the example of 
an IS system. When the I spin multiplets have the same 
phase, for example, a state obtained after a 905 (+4y) I 
pulse acting on an equilibrium system, 0 (0) is 
pEOportional to Lite If the system is then subjected to 


strong S spin decoupling, the effective S spin Hamiltonian 


Hg in the S spin rotating frame is given by 
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Hg does not affect the expectation value <Q> because 0 (0) 
commutes with Hee However, if the I spin multiplets are 
allowed to dephase for a time period, before the start of 
the decoupling, the magnetization components will develop 
a phase difference. At the start of the decoupling period 
(time t = 0), the density operator of the system will have 
a component which is proportional to TS Now 0 (0) does 
not commute with Hs and, aS a consequence, the I spin 
multiplets are affected by the S spin decoupler. 

This unexpected effect of spin decoupling has been 
demonstrated Dbygievitt elsal. (53) using two varlantsror Ja 
2D gated decoupling experiment. The pulse sequences are 
shown in Figure 25. The gated decoupler carbon J resolved 
experiment without proton decoupling during the 
observation period is shown in Figure 25A. The proton 
decoupler is turned on before the carbon 180° refocussing 
pulse in the? evodve1ony period? “The sother varrant of this 
experiment, shown in Figure 25B, has the proton decoupler 
tumned on after the: carbon 180° pulse in, the evolution 
period. It might, at first, be expected that both 
experiments would yield identical results. The pulse 
sequence in #£igure 25A gives Pri," 12211 so 23e1 
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Figure 25. 2D gated decoupler 136 g-§ resolved pulse 


sequences: (A) the decoupler is gated on 
during the first half of the evolution period 
tae (B) the decoupler is gated on during the 
second half; and (C) the decoupler is gated on 
in the middle of the evolution period t,. 


ie 


to the v, axis for CH, CH» and CH3 groups respectively. 
With the sequence in Figure 25B, one obtains an (n+l) x 
(n+1) rectangular array of peaks for each CH, group with 
v, cross Seceion@ritensttres 1:1 and 1s) tor CH, Lslst; 

ie Arterand hele sor CH5, FebeKol, ha eal bess Guile | lesen oak eis Kel 
l:l:1:1 for CH3 with spherical decoupling. 

In the following, the results of a density matrix 
analysis of the behaviour of all three spin systems using 
the pulse sequences in Figures 25A and 25B, carried out 
using the simulation program are reported. The simulation 
results for CH, CH> and CH3 systems are given in Figures 
26-29. A detailed description isypresented for the sch, 
system below. The simulations were carried out with a 
special version of the decoupling subprogram which employs 
the spherical randomization decoupling procedure instead 
of the normal decoupling process. 

For the CH, spin system at thermal equilibrium, the 
density matrix in the basis vdeftined' in Section ~1.4,.1s 


(Figures 27 and 28) 
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INPUT PULSE SEQUENCE FOR GATED-DECOUPLER 2D EXPERIMENTS 


COMM ILLUSIONS OF SPIN DECOUPLING SEQUENCE A 
SPINI1S1 

PULSIA+X 

DCPLS 11 

PULSIB+X 

PRECT1 

PRECT2 

MEASIXYO 

NEXT 

COMM ILLUSIONS OF SPIN DECOUPLING SEQUENCE B 
SPINI1S1 

PULSIA+X 

PRECT1 

PULSIB+X 

DCPLS 11 

PRECT2 

MEASIXYO 

STOP 


QUTPUT FOR PULSE SEQUENCE A ON CH SYSTEM 
Carbon-13 Signals During Acquisition 


SIGNAL = 


+( 0.0 , -2.500)*EXPI(2*PI#( 
*EXPI (2*P1*( 


+( ONO) =25500)*EXPI(2*P 1 | 
*EXP] (2*PI*( 


OUTPUT FOR PULSE SEQUENCE B ON CH SYSTEM 
Carbon-13 Signals During Acquisition 


SIGNAL = 


Figure 26. Input data and output 
decoupler pulse sequences A and B on CH 
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OUTPUT FOR GATED DECOUPLER SEQUENCE A ON CH2 SYSTEM 


Density Matrix After 90(+X) Carbon Pulse 


Oo 

= 
SO ONUIW A= = 
DOD SMH — 


= +( 100.000, 0) |) 
= +| 0. ; 2.500) 
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z= +{ OF ; 2.500) 
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= +( OneT 2.500) 


Density Matrix After the Decoupling period 
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Carbon Signals During Acquisition 


SIGNAL = 


Figure 27. 


+ ( 0.0 , -2.500)*EXPI(2*P1*( 
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+( OO eee 21 U0) ENP (2% Pi: 
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Output data from simulation of gated decoupler 


pulse sequence A on CHy system. 
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OUTPUT FOR GATED DECOUPLER SEQUENCE B ON CH2 SYSTEM 


Density Matrix Before the Decoupling Period 
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Density Matrix After the Decoupling Period 
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Carbon Signals During Acquisition 


SIGNAL = 


Figure 28. 
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Output data from simulation of gated decoupler 
pulse sequence B on CH»s system. 
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OUTPUT FOR GATED DECOUPLER SEQUENCE A ON CH3 SYSTEM 
Carbon Signals During Acquisition 


SIGNAL = a 0.0 , -2.500)*EXPI](2*PI = ( Vee) od Vea) 
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OUTPUT FOR GATED DECOUPLER SEQUENCE B ON CH3 SYSTEM 
Carbon Signals During Acquisition 
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Figure 29. Output data from simulation of gated decoupler 
pulse sequences A and B on CH3 system. 
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where A and 6 have been set to 100 and 5 respectively in 
the simulation program. Here the I spin corresponds to 
the carbon and S spins correspond to protons. In both 
sequences A and B (shown in Figures 25A and 25B 
respectively), a 90°(+x) carbon pulse is applied first, 


and the density matrix becomes (Figurese27 sang 92a) 
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where the superscripts A and B refer to the pulse 
sequences. In Sequence A, the proton decoupler is turned 
On immediately after the carbon pulse and, since all of 
the carbon multiplets are in phase at this point, it 
decouples the proton-carbon interactions during the period 


Tie. sence, atetime Tl, thesdensity matrix (Pigure 27). 1s 


0 (71) = (1/2) , [4-19] 


where a = if exp{2ni VI Tl}. The 180° carbon pulse and 
subsequent delay period Tl, refocusses the carbon chemical 
shifts but retains the modulation due to heteronuclear 
coupling to»the protons. Carbon signals are then acquired 
in the absence of proton decoupling. The signals computed 
by the simulation program (Figure 27) yields the 
conventional 1:2:1 triplet at positions (J, Vit+J), (0, 
VI), (-J, VI-J), where J represents the carbon-proton 
coupling constant. This result is in complete agreement 


with the reported experimental spectrum (53). 
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In the case of pulse sequence B (Figure 28), the 
carbon magnetization components develop phase differences 
during the evolution period Ti. Thej180°% carbon pulse 
following the period T]"interchangesithe labels of the 
carbon magnetization components and the density matrix 


after the 180°s pulse is givensby (Figures28) 
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0 0 0 0 0 0 ety EO 
where b = -i§& exp{2ni(-VI-J)T1}, ¢ = -16 exp{2n1 (-VI)T1}, 
and d= —16expi{271(-VitJ) Tl). it shouid be noted that 


the carbon magnetization components, corresponding to 
density matrix elements pj5, p347 P56 and p7g have 
different phases. proton decoupling is then applied 
during the next evolution period Tl. The density matrix 
at the end of the decoupling period Tl is given by 


(Figure 28) 
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wherere = (—16/3) fexpi—291, J Tl} + 1+ explent a) Tl) and 
£ = -i6. Carbon signals in the absence of proton 
decoupler has been computed (Figure 28), which leads to a 
Ba eaLLay -OL speaks vat. positions (a5 Vito) peur VO 
CiyeViad eC Ope vised) > COs VL); CO,= Vi-J) sade=d = Vide 

(-J, VI), (-3, VI-J) with relative intensities 

Pile eet 4a els lsie Ths 2p Spectrum is vagain: Im agreement 
with experimental results reported (53). 

Tf the proton decoupler is gated on during the 
acquisition period of either sequence A or B, the 2D 
Correclattonamap wit CONSISE, Obs slic ll randaeless cate) 
patterns of peaks at the appropriate carbon chemical shift 
in ¥5 for CH, CH5 and CH3 groups respectively. The 
equivalence of sequences A and B when the decoupler is 
gated on during acquisition, might tempt one to conclude 
(incorrectly) that "illusions" of spin decoupling are 


observed only in coupled spectra. 
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If the decoupling period is placed in the center of 
the evolution period (Figure 25C) rather than at the 
beginning (sequence A) or end (sequence B), carbon signals 
acquired in the presence of proton decoupling will show 
effects of decoupling during the evolution period. Carbon 
Signals are acquired in the presence of proton decoupler 
irradiation. thes 2D contour map and the v1 cross-sections 
obtained from 2-butanol (90% in CDC 13) with the pulse 
sequence C are shown in Figure 30. The 2D map consists of 
aceGipJlet Lor athe. CH group, a quintet, for the CH group 
and septets for the CH, group, with spacings between the 
lines equal to Joqy/4. 

Detisity Matrix analysis Of (thes pu lsemeequcnces: 
app ilzedvon CH (11S); CH> (11S2) and CH3 (I11S3) systems 
are, show in Figuyes 31, 32 and 33 ‘respectively. he 
spherical randomization decoupling procedure has been used 
forvthe Simuliation.  Asbriet description of the CHF results 
(from Figure 31) are given in the following. 

For the CH spin system at thermal equilibrium, the 


density matrix in the basis defined in Section II.4, is 


(Figure 31) 
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INPUT FOR SEQUENCE C (DECOUPLED) ON CH SYSTEM 


ILLUSIONS OF SPIN DECOUPLING SEQUENCE C 


OUTPUT FOR SEQUENCE C (DECOUPLED) ON CH SYSTEM 
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Density Matrix During Decoupled Acquisition Period 
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Carbon Signals During Decoupled Acquisition 
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Input data and output from simulation of gated 
decoupler pulse sequence C (with decoupled 
acquisition) on CH system. 
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OUTPUT FOR SEQUENCE C_ (DECOUPLED) ON CH2 SYSTEM 


Carbon Signals During Decoupled Acquisition 
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Output data from simulation of gated decoupler 
pulse sequence C (with decoupled acquisition) 
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OUTPUT FOR SEQUENCE C (DECOUPLED) ON CH3 SYSTEM 


Carbon Signals During Decoupled Acquisition 
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Output data from simulation of gated decoupler 
pulse sequence C (with decoupled acquisition) 
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6+A 0 0 0 
0 —-6+A 0 0 
oO = a2)) ; [4-22] 
0 0 6-A 0 
0 0 0 =6=A 


where A and 6 have been set at 100 and 5 respectively in 
the simulation program. The application of a carbon 90° 
pulse followed by the free precession period Dl transforms 


the density matrix to 


A g 0 0 

g* A 0 0 
e(D1) = (1/2) ' [4-23] 
- 0 0 -A h 

0 0 he =i 


where g = 16 exp{2n1i (VI+J/2)D1}, and 
he= ro vexp (271 (Vi=-J37 2) D1). 

Proton decoupling is carried out during two Tl time 
periods, in the middle of the evolution period 
(D142T1+Dl). It should’ be noted that the time period 
Dl = Tl, and it is only the program limitation which 
dictates the use of two different time variables. A 
@arbon 180° pulse is applied in the middle of ‘the 
decoupling period (2T1). The density matrix at the end of 


the decoupling period is given by (Figure 31) 
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sites 0 0 0 
e(DIs2TL) = (1/4) ; [4-24] 
mm 0 0 0 5 

0 0 Nes 0 
where 9) = -iélexp{2wi(=VI-J/2)D1} + ‘exp{2nil-VIed/2)pDI}}. 


Equation [4-24] shows that during the decoupling period 
the carbon magnetization components remain fixed in the 
carbon rotating frame and the proton spin quantum numbers 
of the magnetizations associated with the proton spin 
states to be randomized. Further free precession for a 
period Dl serves to refocus the carbon chemical shifts 
developed during the initial Dl period. The density 
matrix during decoupled acquisition of carbon spin 


Signals sas Given by) (Figure 31) 


0 k 0 0 

k* 0 0 0 
cD 2)) = C178) ’ [4-25] 
at 0 0 0 k 

0 0 ke 0 
where ki= —i6 expionl Vie T2)lexpt2m J Dib + 27+ 


exp{-2ni J Dl}]. The carbon signals are given by (Figure 


30) 
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Which gives rise to the 13231, triplet in the 2p map.” ‘The 
program output given in Figures 32 and 33, show that one 
obtains a quintet with relative intensities 1:2:6:2:1 for 
CHeeondeasscptet swith relat ives@intenottlesmde: 2s jugiec jes Zoek 
for the CH3 group, as obtained in the experimental 2D map 
for 2-butanol in Figure 30. 

When the proton decoupler is gated off during the 
acquisition period, the coupled carbon frequencies present 
during the agquisition period T2, gives rise) tora 
(n+1)x(n+1) parallelogram of peaks for the CH, group in 
the 2D map (Figure 34). This observation can be 
rationalized by considering the density matrix for the CH 


qEoup (Figure 65) during coupled acquisition givens by 


0 n 0 0 

n* 0 0 0 
Coie ne a (la) ’ [4-27] 
a8 0 0 0 m 

0 0 ina 0 


where n = -16 exp{2ni(VItJ/2)T2}[1 + exp{2x1 J Dl}], and 
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OUTPUT FOR SEQUENCE C (COUPLED) ON CH SYSTEM 
Density Matrix During Acquisition ( Decoupler Gated Off _) 
pit ( hy Ga) S sx 0.0 , -1.250)*EXPI(2*PI1*( 17.00 


+ ( OOP alle 250) SEXP 2 = Piel 1.00 


*EXPI(2*P1» | 


Ot Shy CA) Ee a0 Wo) HWA OPRAEINP IU Aerasyeal| 1.00 


*EXPI (2*PI» | 
+1( OOM era 2o Ou MEX P Tite PH 1.00 


Carbon Siqnals During Acquisition ( Decoupler Gated Off ) 


SIGNAL = = (( OoM 4 Sooo er Aetaea( 1.00 
+( OO ese 25 0)) EX Pano Pa | 1.00 
*EXPI(2*P1*( 
+( OO sal ZOO) EXP ING 2 Pals 1.00 
*EXPI (2*P I» ( 
+ ( OR ORs 250) SEXP IN 2% P lise 1.00 


Output data from simulation of 


pulse sequence C (with coupled 
CH system. 
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m = -16 exp{2ni(VI-J/2)T2}[1 + exp{2ni J Dl}]. The carbon 


Signals are then given by 


Signa lwar=1(6/4) [expioni( vied 22) 1 + expt2nL 0 Di )) 


7 exXp{ 201 (Vi=J/2) Tite lexp(=2n1 J Diath | [4-28] 


which gives rise to a 2 x 2 parallelogram of peaks with 
Pieensity. lites Valine them pecontouGemap. ms Tne: density 
matrix Simulations for the coupled spectra for CH» and CH, 
Systems are given in Figures 36 and 37 respectively. For 
CH» and CHz systems, one obtains from the calculation, a 
3x 3 "and 4 x 4 parallelogram of peaks, the positions and 
relative intensities of which, are in accordance with 


experimental results for 2-butanol (Figure 34). 


Vas Summary 


Simulations of several pulse sequences applied to a 
variety of spin systems were included in this chapter, in 
order to demonstrate some of the applications of the 
simulation program. The utility of various features of 
the program, such as rf pulses of any given flip angle, 
broadband decoupler irradiation facility and the phase 


cycling feature, were emphasized in the examples. The 
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QUTPUT FOR SEQUENCE C (COUPLED) ON CH2 SYSTEM 


Carbon Signals During Acquisition ( Decoupler Gated Off ) 


SIGNAL = + ( OF07 SOR 833) EXPT (2*Pi = O08 VI 1.00 J )#*T2 ) 
+( 0.0 , -0.833)*EXP1(2*PI «| OORT 100 Um) * 125) 
*EXPI(2*PI*( Vso) wl Eso); 

A sO 5 ST atsheren caterer (Ars ites | 1.00 V1 Wetele) g) Were 9) 
F EXPT 2 Pil ei 2.00 J )*D1 ) 

+ ( Wel) Walker len ior (Actostea 1.00 VI )*T2 ) 
*EXPI(2*PI x | =—Y,00 Y epi ) 

ta OO |, Cenc eer Asie ie O.OMVAL )*T2 ) 
+ ( Os@ 5 SW afcheete yp (eae sea 1.00 V1 e259) 
*EXPI(2*P1*( TO ORCS) = Dis) 

+ ORO On Sosy EXP IN 2s PIs WoO We —}.00 U Wee }) 
*EXPI(2*PI]*( =o s0.0NU) ) * Die) 

+ ( CoO 5 Wahler dP ster( Veo WN oe a} pe) 
*EXPI(2*PI1*( <1) .(00 a) Vet |] 

+ ( 0.0 , -0.833)*EXPI(2*PI*( VoO@ Wil = Wo) ep ere 


Bigqure 36. Output data Gromesimulation of gated) decoupler 
pulse sequence C (with coupled acquisition) on 
CH» system. 
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OUTPUT FOR SEQUENCE C (COUPLED) ON CH3 SYSTEM 


Carbon Signals During Acquisition ( Decoupler Gated Off ) 
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SIGNAL = 


Figure 37. 
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Output data from simulation of gated decoupler 
pulse sequence C (with coupled acquisition) on 
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625) *EXPI(2*P I] *( 
*EXPI (2*P] x ( 
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*EXPI(2*P1* | 
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experiments simulated in this chapter included DEPT, 
selective 2D heteronuclear shift correlation, INADEQUATE 
and 2D gated J-6 correlation experiments which demonstrate 


certain unexpected effects of spin decoupling. 
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CHAPTER V 


DESCRIPTION OF SPIN RELAXATION MECHANISMS 


Velie Snitroduction 


In this chapter, some of the important spin 
relaxation mechanisms and the application of relaxation 
rate measurements in studying molecular motion in liquids 
are discussed. Various relaxation mechanisms and their 
characteristic correlation times are described. 
Theoretical relaxation rate expressions for various 
mechanisms are presented. A brief review of the diffusion 
models for describing molecular motion in liquids is also 
given. Finally the objectives of the project are 


described. 


V.2 Spin Relaxation and Applications of Relaxation 


Studies 


Magnetic relaxation processes are responsible for 
restoring any non-equilibrium state of the spin system to 
thermal equilibrium. The interactions which cause 


magnetic relaxation give rise to a coupling between the 
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nuclear spin system and the lattice. The lattice acts as 
a heat sink and comprises the collection of translational, 
rotational and other degrees of freedom of the molecule 
(65). The coupling between nuclear spins and lattice acts 
as a pathway for energy transfer from the spins to the 
lattice. The coupling between the spin system and the 
lattice is inherently weak (66), but it is via this 
coupling that nuclear spin relaxation occurs. 

A spin system in a static magnetic field at 
equilibrium, has the various spin states populated 
according” to the Boltzmann probability distribution. | Any 
perturbation from this equilibrium state causes the system 
to respond in such a way that the spin populations return 
to thermal equilibrium values, with the rate of recovery 
characterized by the time constants Ty and T5 (66). The 
spin-lattice relaxation time T, is of more importance in 
the study of molecular mechanics in liquids than 1s the 
Spin-spin relaxation time Tor due to the ease of 
measurement of T,. The spin-lattice relaxation time 
describes the recovery of the longitudinal component of 
the macroscopic magnetization (M,) back to its thermal 
equilibrium value, while the spin-spin relaxation time T5 
determines the recovery of the transverse components (My, 
and My) from a nonequilibrium state (66). The 


longitudinal components and transverse components may have 
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different rates associated with them since the dynamics of 
the interaction processes associated with their recovery 
anemditferent. For the intéraction process ‘to: be 
effective in bringing about relaxation, it must be time 
dependent, must involve the spin of interest and should 
fluctuate on a time scale such that its frequency is 
comparable to or greater than the resonance frequency of 
therm spars ls welaxation is “arkected aby high trequency 
processes, while the T, relaxation is influenced by high 
frequency processes as well as those with low frequencies 
suchivas: chemical "exchange ((3\)-" “Hence "in most ‘cases, T> 
can never be longer than T,, however, for molecules 
tumbling rapidly in liquids, they are often equal. 
Relaxation effects have been introduced phenomenologically 


(67) in the classical Bloch equations: 


aM, a ce ~ Mo) 

dat T, 

fell M 

es = wM - Tes FF [5-1] 
ro i ORY, T, 

aM Se oe My 

ac "0: & Be d 


where Mg is the thermal equilibrium value of M, and wo is 


the resonance frequency of the spin. 
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Magnetic relaxation rate measurements provide a 
powerful and highly versatile tool for the study of 
molecular dynamics and structure (3,65,68-71). In the 
field of molecular dynamics, overall molecular 
reorientation (66,70), "internal motion of molecu lar 
Eragments Pin nonrigid systems (72,73) “and relative 
translational motion between molecules (66,70) can be 
investigated using spin relaxation rate measurements. 
Relaxation data can be used to test models (77-80) of the 
microdynamic behavior of liquids. In the determination of 
Organic molecular structure, nuclear spin relaxation 
measurements play a very useful role, particularly the 13¢ 
Spin relaxation rates. Many of these applications have 
been reviewed by Lyerla et al. (68). The 13¢_ly dipolar 
relaxation rates have been used in spectral assignment 
since they can provide information regarding the number of 
protons and ‘their distances to the carbon nucleus. ‘This 
method is commonly used for distinguishing between 
protonated and non-protonated carbons. Another common 
technique which is used in assignment of complex spectra 
is specific deuteration. 

As mentioned above, nuclear spin relaxation can be 
induced by any mechanism which provides a time dependent 
magnetic field at the site of the nucleus (3,66). Hence 


spin relaxation measurements are an excellent tool COL 
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probing molecular dynamics, since molecular motion renders 


time dependence to many of the spin interactions. More 
specifically, the relaxation rates are functions of 
CeLrtaingequantities known assthes corre lation timess= A 
correlation time is the time constant characterizing the 
decay of a particulars molecular property, dsessit is the 
time in which the system loses memory of its previous 


state. Hence, for short correlation times, the system 


loses phase memory quickly and for long correlation times, 


thes phase +memoryepersists for longer periods; i.e. a 
correlation time can be viewed qualitatively as the time 
constant characterizing the exponential decay of the 
correlation function fon the particular property, 

The physical interactions which give rise to various 
relaxation processes have different correlation 
functions. The relaxation rates arising from these 
interactions are proportional to the spectral density of 
the sunteraction, which is just. the Fourier itranstormmaor 
Giemcorre lativon: function.) (66). The spectraisdensity, 
which is a function of frequency, describes how the power 
available for relaxation is distributed along the 
frequency dimension. The spin relaxation rate is related 
to the spectral density by certain constants 
characterizing the intensity of the coupling between the 


spins and the lattice. The importance of the re laxation 
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rate measurements stems from the fact that if these 
coupling constants are known, the experimental 
determination of relaxation rates enables the measurement 
Of the correlation times. To make further connection of 
the correlation times to the way in which molecules 
reorient in liquids, one has to refer to appropriate 
models for the molecular reorientation. The predictions 
of the particular model, about the relation between 
different correlation times, are then compared with the 
experimentally observed variation. Thus relaxation 
studies are of great importance in testing various 


theoretical models describing molecular rotation (70). 


Vee Spin Relaxation Mechanisms 


The spin-lattice relaxation rate of a nucleus is the 
sum of contributions from various relaxation mechanisms. 
For the general case, the time dependent interaction 
Heer) which gives rise to spin relaxation is of the form 


(66) 


ee) upon (BAN oy [5-2] 


z 
qd 


where Iq are spin operators (spin angular momentum 


operators and bilinear combinations of these operators) 
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and Fg (t) are time dependent lattice functions 
characteristic of the spin-lattice interaction. The time 
dependence of Hint is assumed to be random, and the 


relaxation rate, Ry, in the Vimit of strong motional 


narrowing, is a linear combination of the form 


Ree) hee cee can | 6 (Oni Wee [5-3] 
gig e | 
gq 
where Tg is the correlation time over which the 
correlation function Fg(t) Fg (0) has significant 
amplitude, and ag are simple numerical coefficients. Tg 


is the zero frequency component of the spectral density 
Tha) (a), and will be related to the details of molecular 
motion which modulate Fg(t), and [F4(0) |? will be a 
measure of the strength of the interaction producing spin 
relaxation. 

Some of the important relaxation mechanisms are 
intramolecular and intermolecular dipole-dipole 
interactions, spin-rotation interactions, nuclear 
guadrupolar interactions and chemical shift anisotropy. 
in the case Of spin 1/2 nuclei, the dipolar interactions 
and spin-rotation interactions are the major contributors 
to relaxation, while in nuclei which have a quadrupole 
moment, the quadrupolar relaxation almost entirely 


determines the relaxation rate. In certain cases, the 
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presence of internal motion of groups can influence the 
relaxation rates due to some of the above mentioned 


interactions. 


a. Dipolar Interactions 


Dipolar interactions between two spins on the same 
molecule acquire time dependence when the angle between 
the internuclear vector (which is fixed in the molecular 
frame) and the applied field varies with time as a result 
of molecular reorientation. Hence the intramolecular 
dipolar interaction is characterized by a molecular 
reorientational correlation time tg. The functions which 
describe the orientation of the internuclear vector with 
respect to the applied field are the spherical harmonics 
of order two, and t, then describes the decay time for the 
correlation function of the spherical harmonics. The case 
of intermolecular dipolar interaction is more complicated 
since in this case, the time dependence of the interaction 
originates from two sources: (i) as a result of the 
relative translational motion of neighboring molecules, 
the internuclear vector changes its magnitude and 
G@irection with time, (ii) molecular reorientation imparts 
further time dependence to the dipolar interactions. As a 
result, the translational diffusion coefficient Of the 


molecule enters into the relaxation rate expressions. 
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Since ty, decreases with increasing temperature, the 


contribution of the dipolar mechanism to R,} decreases with 


temperature. 
The expression for the relaxation rate, Rie due to 
dipolar interactions between two like spins I and S, has 


been given by Abragam (66) as 
R,1 = (3/2) y4m?z(141) (5) (wy) fea eitoare | ay [5-4] 


where yz, I and wy represent the magnetogyric ratio, spin 
and Larmor frequency of the nucleus I. The functions 
o£) (uy) represent the spectral densities of the dipolar 
interactions. When I and S are unlike spins, the 


relaxation rate of nucleus I is given by (66) 
Rpt = yz27g2n28 (St1){(1/12)T 09) (wp-wg) + (3/2)T |) (wy) + 
(3/4)5'2) (wrtug)} [5-5] 


where S is the spin of the S nucleus. The above 
expressions for the relaxation rates are general in that 
the spectral densities have not been specified as arising 
from any specific mechanism such as random rotation or 


translation. 
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1. Random Rotation 
Using the Debye model for random rotation of 


molecules, the relaxation rate expression for like spins 


is given by (66) 
Ry = (2/5)r7© y#n71 (141) {t9/ (140272) 

+ 4t6/(14+40972)} , [5-6] 
where ris the internuclear distance. In the limit of 


extreme narrowing, i.e. ipo ses <7 lietanc fon spimel/2 


nuclei, Equation [5-6] reduces to 
Rye = (3/2) rg? yr? n? To is—7 1 


For the case of unlike spins I and S, the relaxation rate 


of the I nucleus at extreme narrowing is given by 
6 
R,t = ayy -yeo-Ni-S(St1)74/(3 nara Re 


The correlation time ty is related to the overall 

Fotational diffusion coefficient Ip, by +1, °= 7 SDE) e 
In the case where internal motion of groups is 

present, the internal motion is superimposed on the 


overall rotation of the molecule and the effective 
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correlation time tare LS ae function: of tg and the internal 
motion correlation time t+< The dynamics of internal 


rotation can be studied if vaecan be isolated. The most 


common method of isolation of 4) 1S by studying the 
relaxation rates of two different nuclei such that Tere 
and tg can be determined independently. The mathematical 
ground work for isolating the internal motion from the 
overall rotation was provided by Woessner (72). The 
applications of this method have been compiled in a review 
by Lambert et al. (81). If the internal motion is 
independent of the overall motion and if the internal 
rotation occurs about the principal axis in the molecular 
system, its effect can be introduced into the dipolar rate 
equation by simply adding an internal diffusion constant 
Dea Dew=ad/.)e to the dirtfusion constant D-vabout the 
principal axis (72). For an isotropic (spherical) system, 


the intramolecular dipolar relaxation rate in the presence 


Ofuinternal motion is (72) 
R,1 = (y2y2n2/r°)CA/(6D,) iB) (oDy bse) 
+ C/(6D,+4D;)], [5-8] 


where the geometrical constants A, B and C are given by 
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Raa S cos?A-1)2/4, B = 3 sin2A cos“A, C = (3/4) sinta, 
where A is angle between the internuclear vector and the 
internal rotation axis. Equation [5-8] is derived 
assuming that the internal diffusion is a stochastic 
process and that all internal rotational orientations have 
equal* probability. 

2. Translational Motion 

Based on a simple diffusion model of translation, 
Abragam (66) has given the internuclear dipolar relaxation 


rate as 


Ry = ayn? /(5 ANI y [5-9] 
where N is the number of spins per em, a is the 
hydrodynamic radius of the molecule and D, is the 
translational diffusion coefficient. This simple model of 
translation presented by Abragam (66) is in fact quite 
erude. Lt fails to consider that, in general, 
intermolecular dipolar interactions depend on both 
translational and reorientational motion and that the 
correlation function for the interaction depends on both 
angles and distances (82). Features should also be 
introduced to account for the fact that nuclei are not 
necessarily at the center of the molecule (82). Another 


drawback of the simple diffusion model is that the 
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microscopic detail of the translational motion is obscured 
Since only the gross features of translation are described 
by the diffusion coefficient DiaGecdic me LOL rey ( 33) i nas 
introduced random walk theory to calculate the 
internuclear dipolar relaxation rate and other models are 
due to Hubbard (84) and Oppenheim and Bloom (85). 
Hubbard's model takes into account the effects of 
molecular reorientation on the intermolecular dipolar 


relaxation rate. 


b Spin-Rotatron Interaction 

The spin-rotation relaxation mechanism is of 
importance in spin 1/2 systems. In this case, the 
rotational anguilar-momentum of the molecule produces 
fluctuating magnetic fields at the site of the nucleus 
Since the electron density in the molecule behaves as a 
charge system with an angular velocity. Hence the 
fluctuations of the magnetic fields so produced are 
determined by the time dependence of the angular momentum 
of the molecule, with the correlation time tas in the 


case of spherical molecules the relaxation rate is given 


by (86,87) 


Pysn KULG 1 5/Nse, [5-10] 
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where I is the moment of inertia and C is the LSOLLOD1IG 
part of the spin-rotation constant. For symmetric top 


molecules the relaxation rate expression is given by (88) 
= 2 
eet (Cyber 2C oun) ay (ances [Ss 


where the symbols || and 1 denote the parallel and 
perpendicular components of the respective quantities. It 
has been shown by Hubbard (84) that when la “2h cand jin 
extreme narrowing, the spin-rotation relaxation rate is 
proportional to t;. Hubbard also showed that Ege bE 
DROpCEtLOnal to De) T- ty increases with increasing 
temperature and hence the spin-rotation mechanism gains 
importance at higher temperatures. 

ty contains information on the molecular collision 
frequency. The spin-rotation interaction is more 
important in small symmetrical molecules with weak 
intermolecular interactions since they tend to possess 
rather large angular velocities. Determination of the 
spin-rotation contribution to the total relaxation rate is 
often difficult since the components of the spin-rotation 
interaction tensor are not usually known. Flygare (89) 
has proposed a method which circumvents the task of 
determining the elements of the spin-rotation tensor 


experimentally. This method takes advantage of the 
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relation between spin-rotation constants and chemical 
shift anisotropy. Spin-rotation constants can be related 
to chemical shifts since both depend on the electron 


densaty distribution in the molecule. 


Ce Nuclear Quadrupole Interaction 


For nuclei with spin I > 1/2 the nuclear quadrupolar 
contribution dominates the total relaxation rate 
(3766 ,70) .  Nucletywith spine! 2.1/2 possess anvelectric 
quadrupole moment (produced by the non-spherical symmetry 
of the nuclear charge distribution) which interacts with 


the electric field gradient produced by the surrounding 


electric charges at the site of the nucleus. The electric 


field gradient becomes time dependent due to molecular 
motion. The quadrupole interaction vanishes when the 
nucleus has an electrical environment which has cubic or 
higher symmetry (70). 

Quadrupolar relaxation is an intramolecular process 
and provides a very convenient method for determining the 
reorientational correlation time ty in many cases. For 


the case of extreme narrowing, the expression for 


quadrupolar relaxation rate has been given by Abragam (66) 


as 
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where I is the spin quantum number, (e2qQ/n) is the 


quadrupole coupling constant, and n is the asymmetry 


parameter 


P a*v/ax- - a*v/a 2 
= oS [ES5=134 
OV jez 

where V(x,y,z) is the electric potential due to the 
electrons and {x,y,Z2} form a molecular coordinate system 
centred at the nucleus with z axis chosen to be the 
principal axis of highest symmetry for the field gradient 
tensor. 

For cases where the asymmetry parameter n and the 
nuclear quadrupole coupling constant can be determined, 
Equation [5-12] provides a more convenient and accurate 
method for determining tg than the dipolar relaxation rate 
expression in Equation [5-7]. This is due’ to the fact 
that for nuclei with a quadrupole moment, quadrupolar 
interaction dominates strongly all other mechanisms and 
there are no intermolecular contributions to the 
relaxation rate. 

If the quadrupolar nucleus is situated in a group 
which executes internal rotation, the effect of the 
internal motion is superimposed on the overall 


reorientational motion. Following Woessner's approach (72) 


usecys asenl inn i routine a aot ‘2.¥ ‘s} bas ns 


> 


im? oot anh oe eases vrsepars a - KeoterlY #taiw | 


voi 6 api ape se wish’ os 9m hota 


iA uF ier sar ne aigeentus saa ws hla arin 


eo 
cyons 


at! tbat 0 aay yr eempyes Gar ceahie eon? 10% 


eee 
fraud Lie Sem igo acl. tS" Snasuntrs grttaves 4 : 1 
vodiaitene lay ntaavete SHER 14, elon nie fo iat 
oie: nfiackeiia deleekts onhs rods a? shaunsegsd 708 a 
— 


cael, 399 45 enh, oh aget NED! netttesonll soteawagee! 
seq piepy, trem “come Tp ana ¢ tae intun 203 foe aa 
Lan gee beteitaeh ewes hte vi ni ao seniino® stata 
aie og aa dag a | 


auane 4 ae niece alae 
antsy acme | 


123 


for internal motion modulation of dipolar interactions, 
the reorientation correlation time in Equation [5-12] can 
be replaced (73,126) by an effective correlation time as 


is done in Equation [5-8]. 


V.4 Brief Review of Rotational Diffusion Models 


mne earlvest Of the models for rotational daifttusion 
in liquids was proposed by Debye (90) to describe 
dielectric relaxation. The Debye model of diffusion was 
applied to spin relaxation by Bloembergen, Purcell and 
Pound (91). In the Debye model, the motion of molecules 
in liquid is pictured as the rotation which macroscopic 
spheres undergo in a continuous viscous medium. The 
molecules are acted upon by torques produced by the 
Brownian motion in the fluid. The model further assumes 
that the duration of each rotational step or the time 
between the action of two consecutive torques, is small 
compared to the time required for the molecule to change 
iis Orientation appreciably. This impltes that the 
molecule experiences many changes in angular momentum 
before it changes its orientation. Since the duration of 
a rotational step is characterized by the angular momentum 
correlation time t,; and the overall rotation by the 


reorientational correlation time ty, the Debye model 
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assumes that Tg >> tz. Hubbard (84) has used the 
classical Langevin approach and shown that, for spherical 
molecules, the reorientational correlation time Tg Can be 


related to the angular momentum correlation time ugg 1D 


V6 ty = Te / OKT ’ [5-14] 


where I.y is the average moment of inertia of the 
molecule, kK is “the Boltzmann constant. and Tas the 
absolute temperature. The assumption that ty >> ty is 
especially objectionable when dealing with the 
reorientation of spherical or nearly spherical molecules 
for which it has been found that To andet7 are Of "thie same 
order of magnitude (79). Thus the Debye model fails to 
describe molecular reorientation in the whole range of Tg 
variation with ty;- 

When ty and tz are of comparable magnitude, the 
molecule reorients through a large angle ina single 
arefusion step, just as in the dilute gas where molecules 
rotate freely between collisions. The concept of free 
rotation and its effect on ty was recognized by many 
workers (77,92-94), but many of the theories developed 
were only approximate and were not valid over the entire 
range Of ty, and t7- 

The first successful theory which did not incorporate 


the assumption of small angle diffusion was due to Gordon 
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(77). This model has been called the extended diffusion 
model and he applied it to the case of linear molecules. 
In the extended diffusion model, no restrictions on the 
Size of reorientational steps were imposed and the model 
follows the succession of diffusive steps in microscopic 
detail. The extended diffusion model has been extended to 
apply for spherical top molecules by McClung (79,95,96), 
Mountain (97), McClung and versmold (80); for symmetric 
top molecules. by McClung (98), Fixman;and Rider (78).,, and 
other workers (99-102); and to asymmetric top molecules by 
Leiceknam et al. (103). 

In the extended diffusion model, the molecular motion 
is examined in detail during each successive diffusion 
step. It is assumed that the molecule rotates freely 
between collisions. The collision events are thought to 
be instantaneous and act as terminators of the diffusive 
Steps}. AS the! result tof a collision event, ~ehe motecu lar 
angular momentum is randomized. This leads to two extreme 
Situations: ~ 1) a case where both the direction and the 
angular momentum are randomized, referred to as the J- 
diffusion model and, 2) a case where only the direction of 
the angular momentum is randomized - the M-diffusion 
model. The durations of the diffusive steps are assumed 
to be random and to follow a Poisson distribution with the 
characteristic time Tye since each collision event 


randomizes the components of the angular momentum. 
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It has been shown (79) that the extended J-diffusion 

model gives the same results as the Debye model when 

Tg >> tz and reduces to the perturbed-free-rotor model 
(104) when tg << ty. In the Debye limit (tp. >> tz), the 
variation of tg with respect to v7 ©Or the V-dvetusion 
model agrees with Equation [5-14] derived by Hubbard 
(S6).08 4 this limit, Tg has an inverse t; dependence. 
Whentun <<s 17,8 1. eat the perturbed-f{ree-roton Namits o7 
dilute gas limit, tg has a direct dependence on ty (104) 


and 


TQ = tz/5 . [5-15] 


The calculations done (79) using the extended diffusion 
model show the transition from the Debye limit to the 
perturbed-free-rotor limit. Calculations using the M- 
diffusion model show that, in the Debye limit, To is 
inversely proportional to ts, but the DHoporteronality 
constant is different from that in the Debye model. This 
Nas? been explained “(/79,105) as*due to’ the (fact thatsDebye 
model assumes randomization of the magnitude of the 
angular momentum at each step termination which isi not 
assumed in the extended M-diffusion model. When the 
anisotropic intermolecular interactions are very weak, the 


M-diffusion model may be a better approximation because 
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randomization of the magnitude of the angular momentum at 
each collisional event is expected to be less probable. 
The J-diffusion and M-diffusion models are equivalent in 
the didute gas limit, and To exhibits a linear 
relationship, with ty 
The application of the extended diffusion models has 
been carried out by comparing the predicted variation of 
T, with tz, With the observed variation determined using 
spin relaxation measurements. Independent determinations 
of Tg and tz can be made experimentally from relaxation 
time measurements of different nuclei in suitable 
molecules. An exhaustive survey of these applications has 
been given by McClung (105). The extended diffusion model 
has been applied to linear, spherical top and symmetric 
top molecules. In most of the applications, t,; was 
determined from the measurements of spin-rotational 
contributions to the relaxation rates of spin 1/2 nuclei 
Such as list 19 and 3lp, and Tg was determined from 
relaxation rates of quadrupolar nuclei. In most cases, 
the J-diffusion limit of the extended diffusion model was 
found to give results in quantitative agreement with the 
experimentally observed ty variation with tj. The M- 
diffusion model did not seem to agree very well with most 


of the experimental observations. The J-diffusion model 


results however, as expected, did not agree well with 
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experimental observations on certain linear molecules 
which have strong intermolecular interactions such as 
hydrogen bonding. The reasonably widespread success of 
thevJ=diftfusion model “is not surprising since it is 
capable of yielding results equivalent to the classical 
Debye diffusion model when ty >> tz and those similar to 
the results given by dilute gas model at the limit 
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Mero Objectives Of the Project 


The primary objective of this work is to investigate 
the relative importance of various relaxation mechanisms 
in the molecule hexamethylphosphoramide, O=PLN(CH3)5]3, 
(HMPA) and its deuterated analog (HMPD). Relaxation 


31p, 14x 2H and ly nuclei have been 


rates for the 
studied over a wide range of temperature. The variation 
ofthe translational” difttusion coerficient with 
temperature for hexamethylphosphoramide has been measured 
to facilitate separation of the contribution to 
relaxation rates from intermolecular dipolar 
interactions. It was intended that the reorientational 
correlation time Tg and the angular momentum correlation 
time Tz be determined from the relaxation data. With the 


knowledge of the coupling constants in the relaxation rate 
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expressions, from other independent sources, the 
correlation times characterizing various relaxation 
processes for different nuclei, can be determined. The 
knowledge of the variation of the reorientational 
correlation time Tg With the angular momentum correlation 
time tz; can then be used to test the applicability of the 
extended diffusion model in describing molecular 


reorientation in hexamethylphosphoramide. 


V2.6 Summary 


A discussion of various spin relaxation mechanisms, 
their characteristic correlation times and the application 
of relaxation data in the study of molecular motion has 
been presented. Quantitative expressions for various 
relaxation rates were given. Two models describing 
molecular rotational motion - the Debye model and the 
extended diffusion model were described. The objectives 
Of fhe project wwere*salso laid out. The next chapter deals 
with the experimental aspects of spin-lattice relaxation 
rate determination and the measurement of translational 
diffusion coefficients. In the following chapter the 


results and data analysis for HMPA and HMPD are presented. 
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CHAPTER VI 


EXPERIMENTAL 


View sine rogcuction 


This chapter discusses the experimental aspects of 
the determination of spin-lattice relaxation times and 
translational diffusion coefficients. The instrumentation 
used for the measurements is described in brief. Two 
methods have been adopted for measuring eeke axa ton 
times: (i) the inversion recovery experiment and (ii) the 
triplet sequence which was used for measuring relatively 
long relaxation times. The last part of the chapter 
describes the measurement of translational diffusion 
coefficients by the stationary (time independent) field 


gradient method. 


Vu.2 sample Preparation 


Hexamethylphosphoramide (Aldrich Chemical Co.) was 
dried by storing over molecular sieves (Fisher Scientific 
Ca. Type 3A). It was then vacuum distilled anto) a, clean 


NMR tube attached to the vacuum assembly. The sample was 
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then degassed using the pump-freeze-thaw technique three 
times. After degassing, the NMR tube was sealed under 
vacuum. An identical procedure was used for preparing the 
sample of the deuterated compound. For the measurement of 
3lp T,; in HMPD, a chemical means of deoxygenation of the 
sample using [cot bipyedilclon). with sodium borohydride 
(106) was employed. In an oxygenated solution these 
compounds give a pale yellow-brown color which changes on 
removal of oxygen to an extremely oxygen sensitive intense 
blue color (106). This mixture was added to the HMPD 
sample and the HMPD was vacuum distilled into an NMR 

tube. The 3!p relaxation times in this sample did not 
differ significantly from those in the sample prepared 


with the pump-freeze-thaw degassing technique. 


Views Instrumentation 


The measurement Of T) values was: carried out with a 
Bruker SXP4-100 high power pulsed NMR spectrometer and a 
14 k Gauss Varian Associates electromagnet and V3506 
magnetic field flux stabilizer. The Bruker SxXP4-100 
pulsed spectrometer was operated under the control of a 
Nicolet 1180 minicomputer along with a Nicolet 293A 
programmable pulser. The data were stored on a hard disk 


using a Diablo Disk Drive Unit under the control of the 
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Nicolet 1180 minicomputer. Analysis of the data was 
Carried out on the minicomputer in the BASIC language and, 
to a lesser extent, in the Nicolet 1180 Assembler 
language. 

The spectrometer has four independent radio frequency 
(rf) channels which can be gated on or Oa by the Nicolet 
293A programmable pulser. The phase of the rf pulse can 
be continuously adjusted in each of the four independent 
Channels. Tuning the spectrometer to the required radio 
frequency is carried out by adjusting two variable 
Capacitors on the high power amplifier and on the probe 
arm. 

The detector part Of the = circuitry consists of a 
preamplifier followed by an amplifier which can be 
operated in both phase sensitive and diode detection 
mode. The amplifier, when operating in the diode 
detection mode, gives a signal which is the square root of 
the sum of the absorption and dispersion signals. The 
signal after diode detection is unaffected by moderate 
Changes in magnetic field or radio frequency. Phase 
sensitive detection has the advantage of linearity of 
response over the full dynamic range of the amplifier and 
better signal-to-noise ratio, particularly at low signal 
intensities. The detector also contains a rf filter with 


band width adjustable from 100 kHz to 100 Hz. 
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The pulse lengths corresponding tc. 890° “and “eget 
angles are determined prior to the T; experiment. Using a 
[wait - pulse - acquire] sequence the diode detected 
Signal observed should be a maximum for a 90° pulse and it 
should be a minimum for a 180° pulse. The spectrometer 
was tuned prior to determination of the pulse lengths. 
Tuning was carried out by adjusting the tuning capacitors 
such that, with an initial pulse length much less than an 
approximate 90° pulse, the signal observed was a 
maximum. The pulse length was then increased so that the 
amplitude of the FID decreased to about three-fourths of 
its maximum value and the capacitors were again adjusted 
to give a minimum signal amplitude. Accurate 
determination of the pulse lengths was carried out uSing a 
computer program based on an optimization algorithm 
(107). This program attempts to minimize the signal 
obtained after a 180° pulse. The 90° pulse length was 
then taken as half of the 180° pulse duration. 

In experiments requiring rf pulses of different 
phases, the relative phases of the pulses on each pulse 
channel were adjusted manually by observing the FID in the 
phase sensitive mode after a 90° pulse. The phases of the 
last three channels are always in reference to the phase 
of the first channel. These relative phases were adjusted 


to the desired values as required in particular 


experiments. 
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In all temperature dependent studies, the temperature 
of the sample was controlled by a Bruker B-ST 100/700 
temperature control unit. A Doric Trendicator 400 digital 
thermometer using a copper-constantan thermocouple was 
used to measure the temperature. The temperature was 
measured before and after the T, determination and the 
variation in temperature was less than + 0.5°C. Sample 
was allowed to equilibrate for 15-20 minutes after 


insertion into the probe. 


VI.4 Measurement of Ty Relaxation Time 


a. Inversion Recovery Experiment 

Relaxation times which were not extraordinarily long 
were meaSured uSing this conventional 180-1t-90 pulse 
sequence (3). The experiment can be represented by, 
Wem co oOem 1290 Watt (Ol tho Oa) Swae 
(SOS ae N90) eee) Where 7. and A are the initial delay 
and the increment to the delay respectively. The wait 
time is approximately 5-6 times the T, value. Typically 
abouteten 7 values swere “used ‘such’ thatyrg + LOA = 93T}- 
The initial value +t, was chosen as 0.7 times the 
approximate Ty, Value. @The PID after sthes902 pulse was 
detected in the diode detection mode, digitized and the 


intensity of the magnetization was calculated by 
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integrating over about 50 data points in the initial part 
of the FID. The total acquisition time for the complete 
FID was typically in the range of 5-25 mS. The signal 
intensity M(t) for each t value was then fit into an 


equation of form 


M(t) =A +B exp (=o) , [6-1] 


using a non-linear least squares analysis algorithm. The 
computer program was designed to make repetitive T, 


measurements and to calculate the variance and the 


standard deviation. The diode detection mode was employed 


for signal acquisition since it was independent of 
moderate variations in the magnetic field over the period 
Grircimerrequiredtorcarny out the repetitive: Ty 


measurements. 


b. Toipley Sequence 


The triplet sequence (108-111) was used for 
determination of T, values which were relatively long. 
This sequence is more difficult to carry out and has more 
stringent requirements concerning the accuracy of the 
pulse lengths. The advantage of using the triplet 
sequence is the saving in time, since the wait time is 


avoided. The triplet sequence can be represented by 
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180 ,-1-[903-6-1802,-5-908-A],,, where 6 must be much 
smaller than Ty and was typically about 1 ms, and 
A = acq.time/n, typically about 100 ms. The part of the 
pulse sequence enclosed by the brackets is repeated n 
times (typically n = 256). Data were acquired after the 
first 90, pulse of the pulse train enclosed in brackets, 
using the phase sensitive detection mode. Thus 256 data 
points were collected during the acquisition time, one 
data point for each pass through the seguence in brackets. 
The triplet sequence is basically similar to the 
inversion recovery sequence, except that after the signal 
acquisition period, instead of allowing the magnetization 
to evolve back to the thermal equilibrium value, the 
dephased magnetization vectors are refocussed by the 1802, 
pulse and rotated back to the -Z axis by the second 907 
pulse. This process is then repeated a number of times, 
thus collecting a data point (M(7),+) during each spin 
echo formation. The data points were fit to Equation [6- 
1] using the non-linear least squares analysis program. 
Typical precision of T, measurements was about 15%; 
however, when the T, values were very long as in 3lp in 


HMPD the reproducibility was between +5% and +102. 
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VI.5 Measurement of Translational Diffusion Coefficients 


The translational diffusion coefficient of hexamethy1l- 
phosphoramide was measured by the stationary field 
gradient technique (3,112-114). The amplitude of the spin 
echo from equivalent spins is modulated by diffusion of 
spins along a magnetic field gradient in the sample. The 
etfective T> relaxation time Tere in tne presence Of-ea 
field gradient, contains an extra term involving the 
diffusion coefficient D,- Hence for the experimental 
determination of Dy, the transverse relaxation time is 
measured using the spin echo technique, in the presence 
and absence of a magnetic field gradient. Both stationary 
and time dependent field gradients can be employed 
(3,115,116); however the experiments carried out here were 
done using the stationary field gradient. 

In the first part of the experiment, the relaxation 
time T> of the sample was determined in the absence of any 
field gradient using the Carr-Purcell-Meiboom-Gill (CPMG) 
pulse sequence (112,117). The CPMG pulse sequence can be 
represented as, 90,-1/2-[180)-tIh, where a train of echoes 
is generated by the 180,, pulses and the t time 
intervals. The 90; pulse creates transverse y 
magnetization which is allowed to dephase for a time 


period t/2. A 180), pulse is then applied which refocusses 
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the chemical shift precession and an echo is formed at 
time +/2 after the 180) pulse. This echo is sampled by 
the A/D convertor of the 1180 computer. A further train 
of echoes is formed by successive 180, pulses: and 7 
delays. The amplitude M(t) of the spin echo at time t in 


the absence of magnetic field gradients is given by (3) 


Mice tab exp (Sty Ta, [6-2] 


where A and B are fitting parameters. In the presence of 
a field gradient of magnitude G, the echo amplitude is 
moditiediias «(3,113 


M(t) = A + B exp[-t(— + ae [6-3] 


2 
where: Di®as* the diffusion coefficient and’ 1'as the ime 
between the 180° pulses. 

The signal was sampled in the phase sensitive 
detection mode at the maximum of the echo amplitude. For 
the case where no field gradient was applied, the Ty value 
was determined by fitting the data to Equation [6-2] using 
a non-linear least squares program. The CPMG sequence 
compensates for pulse length errors by using a 90° phase 
shifted channel for the 180° pulses. Trial runs of the 


experiment were carried out during which the relative 


208 


=) Rees ON 
> 
= 


‘a ; 


2 a tah 
1 Sntqrise at aaa ett .eelog OGL = 
nenut vedshat€ sabe OE wat eva 
~ sila soit Gt senonedoum is heh he 
sake seSetow sedi ods 3a Gale ehanthlqns wT (se 
Lob we nevdyel eis peer 


" eae! / +a 


=o ; jemrrees ae 
¥ 
- 
ts evndeeap ond of .stedenesen onizeid ota a hae A wae 


wy) Stusrigne’ cise sity ,t b seer Lemans 20 chee: bisth 


ten 6) da beltinam Ls 


eh = 


(t-ad et! eet + geri a+ a “ak 


22) orig. al, © SH sine sh et Hg 
cao tads “ORL ala ndewsed ~ 

avidiense eaeiy arty DP ichins iw teeta @ 

vedic eS ae ee 
cutav. of mitt (Redtque wow soe ap oben on 
pazay (2-¥) nereeups ae — A} gies: : 38w 


je 
Recsiscive atic dental mer prea 


=i 3 
ise ee reer 
err i ev 


_ 


ela cemau 
~~ 


phases of pulses were adjusted using the spectrometer rf 
phase control knobs and the pulse durations were optimized 
by adjusting the potentiometer knob settings on the 
Nicolet-1180 computer which control the pulse lengths used 
by the computer program. The criterion for adjusting the 
phases and pulse lengths was the symmetry of the spin 
echoes displayed on an oscilloscope screen. 

The second part of the experiment was carried out by 
measuring the effective relaxation time in the presence of 
a stationary field gradient. The field gradient was 
generated by passing current through a pair of copper 
coils mounted on the aluminum side plates on either side 
of the probe. The coils consisted of several turns of 
wire and were arranged so that the current flowed through 
the two coils in opposite directions. The current to the 
coils was supplied by a Union Carbide 6 volts No. MS31 
cell, and was gated on or off under computer control in 
synchronization with the CPMG pulse sequence, using 
programmable level line 7 of the Nicolet 293A pulser. 

The field gradient was calibrated using samples of 
distilled water and benzene whose diffusion constants were 
Guvem Lm iiterature (1187119) . The field gradient G 
(Gaviss/cm), aS: a. function of the current. | (mA) passing 


through the coils, was obtained using 


209 


inna d ioe Saye on aw 
sit HD fsieeylay 
nwou atzgomd va ud edit Lovseot.ainkee, weeeighos nd bnantco say 7 
ato godgeuihe 16% nckeeessD SAT” megONG sedeqND? ond. 
alae! a4 0 ee ats may wayne satu noe onesie | 
«feeb Y etrnnrd) Pied ihr £65 Bevel calh obddae: 
e@¢ tua AA l TID eABW Sais Neha ewy Va, dwg baowee wT ive 
io ooneeesg sda oF omit mid Anat oe av) snatte ed? pntsuesom 
sev. #nei hess dent atT a énedtene Siod2 ranedecem - 
re eet djiveaits | greemrpee ape bene \el Ootarerep: 
ope qendte wo ete shed et ad? we Gednuem spo 
to eoiid Letovo® So Redetetian alae en? .edo3g edt Yo 
Aque ces bewela Inver sas aan ig Hepieaze ogy bre ore 
uty of! 4AgT (UR BNT <enol Tae it sikavege, ot elioo ows axtd 
LkuM: «det alowed aa Lae!) Hobe © Ya belt ggue zav afsee 
yt (oa ea REM Memes. TmIM te Soong esau new ope yhS OR) | : 
Ortsl  .ebeewpee wee) pied! agus te si 4azinowls oe 
Sel Og WEES IehAote ake, Teall Lovet ‘ebdnam 
ay ne ee 2) 


where m = 0.155 + 0.005 Gauss cm7~!ma7!. The cPMG 
experiment was first carried out in the absence of the 
field gradient to determine TD el oeee was EGen Obtained 
from the experiment done in the presence of field 
Gradients, and fitting the data to Equation [6-2]. “The 
field gradient was then calculated from the relation 


(Equations [6-2] and [6-3]) ~ 


The field gradient G was determined for various values of 
current, through the cells, ‘and a calibration chart or 
field gradient against current was constructed. The 
current passing through the coils was measured with a 
Simpson digital multimeter model 461. 

The effective relaxation time of hexamethyl 
phosphoramide was then determined in the presence of the 
field Gradient. Hrom the value of T57--)the diutfusion 
coefficient D, was calculated by substituting the Known 
values of Ty into Equation [6-5]. The diffusion 
coefficient of hexamethylphosphoramide was measured over a 
temperature range 300-440 K. The precision of this method 


of measurement was approximately +10%. 
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VI.6 Summary 


A description of the experimental aspects of spin 
relaxation time determination was given in this chapter. 
The NMR spectrometer system used was described briefly. 
The spin lattice relaxation time was measured using two 
techniques: (i) inversion recovery experiment and (ii) 
triplet sequence. Both of the above methods were 
described briefly. It was fo llowed bypamdiscussion of the 
measurement of diffusion coefficient employing the 
stationary field gradient technique. Brief accounts of 
the field gradient coil calibration procedure, CPMG 
experiment and relevant equations were given. The results 
of the experiments and discussions concerning the 
separation of relaxation rates into contributions from the 


individual relaxation mechanisms are presented in the next 


chap te i. 
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CHAPTER VII 
RESULTS AND DISCUSSION 


Velo” Introduction 


The results of the nuclear relaxation and diffusion 
experiments on HMPA and HMPD and the analysis of the data 
are =presented in this chapter. ‘Section Vilw2 deals. wien 
the T, results for the 31p Renee in HMPD and the 
determination of the angular momentum correlation time. 
This is followed by the analysis of 14y TD oteast ne HMe A, 
imesection Vil.3. Section VIl.4 deals with the 2H 
relaxation data analysis. The T, data for 31p nuclei in 
HMPA, the translational diffusion coefficient results for 
HMPA, and the delineation of the contributions to the 31p 
relaxation rate in HMPA are presented in Section VII.5. 


The conclusions derived from the study are presented in 


Secte2on Vil «6 


vir.2 3p Tae CEMPD) Results and Analysis 


The variation of the spin lattice relaxation time of 
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315 in the deuterated analog (HMPD) with temperature was 
Studied from 295 K to 460 K. The results are shown in 
Figure 38 and are tabulated in Table 2. The relaxation 
times were relatively long and the measurements were 
carried out using the triplet pulse sequence described in 
Chapter VI. The relaxation time at room temperature was 
approximately 48 sec and it decreased to about 12 sec at 
260K. 9 To atgood approximation, the 3lp relaxation in the 
deuterated compound can be considered to arise solely from 
the spin-rotation interaction, since the dipolar 
interaction between 3!p and 2H would be relatively weak. 
The observed variation of the relaxation rate with 
temperature in Figure 38 shows the typical characteristic 
of spin-rotation relaxation: an increase in relaxation 
rate with increasing temperature. In order to obtain an 
approximate measure of the 31p_-2y dipolar contributlons,.-a 
calculation using the reorientation correlation time from 
the 2p relaxation rate (see Section VI1.4) was carried 
Sum pewii ch ayieldedsayvaluic of 637 0seemtor Sthesdipotar ei, 
at room temperature. The experimentally measured T, was 
approximately 48 sec at the same temperature. The 
contribution of 31p-2y dipolar mechanism to the overall 
relaxation rate at room temperature, where the dipolar 
contribution is expected to be larger than at high 


temperatures, is only slightly larger than the uncertainty 
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Observed relaxation rates of 31p in HMPD. The 


Figure 38. 
ay relaxation times were measured with the 


triplet pulse sequence. 


in the T, Measurements. Similar calculations of the 
dipolar interactions with 1]4y nuclei, yields a relaxation 
time of approximately 830 sec for 3!p nucleus at room 
temperature. Hence dipolar contributions to the 31p 
relaxation in HMPD can be neglected. 

It is possible to extract the angular momentum 
correlation time t, from the observed variation of 3lp ai 
values with temperature in the deuterated compound using 
Equation [5-11], provided that C, and C,, the parallel and 
perpendicular components of the spin-rotation interaction 
tensor, and I and Lis the corresponding moments of 
inertia of HMPD are known. The moments of inertia I, and 
Peo, AMD were determined from the x-ray crystal 
structure data reported by Le Carpentier, Schlupp and 
Weiss (120) for the compounds MoO( 05) 9PO[N(CH3) 9] 3H20 and 
MoO(05) »PO[N(CH3) 5] 3Cc5HoN. Average values of the bond 
lengths and bond angles in the two compounds, the C-H 
ar1stance; (1. L011 x 1078 cm) and the H-C-H tetrahedral angle 
(109°28') were used in the calculation of moments of 
inertia. The average moments of inertia were calculated 
to be: I, = 1.31 + 0.02 x 10737 g cm@ and 
I, = 0.98 + 0.01 x 10737 g cm2. 

The components of the spin-rotation constant tensor 


were determined from 3!p chemical shifts utilizing the 


relationship between spin-rotation constants and the 
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paramagnetic part of the chemical shift (89). Flygare and 
Goodisman (121) have shown that the paramagnetic part of 


the chemical shift and the components of the spin-rotation 


Tensor are related by 


x fe e 
Say = Ga, ete atom) + enc (My le + Myyly + Movie) ; 


Ethan) 


where o,y is the average chemical shift, on ( free atom) is 
the free-atom diamagnetic susceptibility, e is the 

electron charge, m is the electron mass, c is the speed of 
Hedecy ty 


We Le the nuclear (magneton, g.is) tie=nucleaneg— 


PaclOc sand Maa.) and Mj, are the components of the 


yy 
spin-rotation tensor (in ergs). The spin-rotation tensor 
C (in radians/sec) is related to M by M = AC. Equation 


[7-1] can be rearranged to the form 


me 28 d > 
Cereal On On, Eee Ate OM oly [elie dat me Xp Vn Zy 


2 


where y is the nuclear magnetogyric factor and je is the 
Bohr magneton i = Je|n/[2me]). 

The parallel and perpendicular components of the 
spin-rotation constant for HMPD, C, and C,, where C, = Cy, 


= - ined from Equation [7-1] 
and C, = Cy, Cyy can be determine q [ 
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if the anisotropic components of the chemical shift 
tensor, o, and o, are known. The 31p Specerumeon solid 
HMPA at 77 K shown in Figure 39 was measured by Dr. J. 
Ripmeester at the National Research Council, Ottawa. The 
resonance frequency of 31p was 72.87 MHz and the spectrum 
was measured with a sweep width of 125 KHz collected in a 
memory Size of 512 points and zero filled to 4096 points 
before Fourier transformation. The delay between 
Successive scans was 160 sec and four scans were co-added. 
The spectrum shows the anisotropic line shape resulting 
from a system with axial or near axial symmetry. The 
components of the chemical shift tensor o, and o, were 
Obtained by fLitting the spectrum to an anisotropic line 
shape function with axial symmetry and Gaussian line 
broadening. The simulation was analogous to that for 
paramagnetic line shapes of polycrystalline substances 
reported by Ibers and Swalen (122), who used Lorentzian 
rather than Gaussian broadening. The chemical shift 
anisotropy was determined to be o, - o, = 163.8 + 0.1 ppm. 
The average chemical shift of 3lp in HMPA with reference 
to PH; given by Pople et al. (123) and the chemical shift 


of PH3 relative to the free phosphorous atom given by 
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Cay(HMPA) - o5,(PH3) = -261 + 1 ppm , 
[7-3] 
oS Usa Sayil free atom) = -366.4 + 0.1 ppm 
give 
O ay (HMPA) ~ pee oELee atom) = -628 + 1 ppm. [7-4] 


Equation [7-4] and the measured chemical shift anisotropy 


in HMPA given by 


were used to obtain the components of the chemical shift 


tensor relative to free atom value as 


o, (HMPA) - Beat ree atom) -518 + 1 ppm , 


So) 


6, (HMPA) - ogyi(free atom) = -682 + 1 ppm 


Equation [7-5] and the calculated moments of inertia were 


substituted in Equation [7-2] to obtain 
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The angular momentum correlation time Ty was 
determined from the 31p relaxation rate using Equation 
[5-11] and the above values Cir Cy, 1 and I). The values 
of ty, together with the values calculated from the least 


Squares fit to the form 
eX iB yl): by -6e 


are given in Table 2. The non-linear least squares 
program determined the fitting parameters to be: 

Pee eon 005m 100 seceandes 66700 me s0lK me Them ane 
Oa ac) 300 7Ki is 471 10714 sec and it increased to 1.1 


x 10 22 sec at 460 K. 


WIP a ln T, Results and Analysis 


MhesVariablon Ofe iy; of the 14N nuclei in HMPA was 
measured over the temperature range 300 K - 445 K. The Ty 
experiments were carried out using the inversion recovery 
method described in Chapter VI. The experimental results 
are given in Table 3. The measured T, values showed a 
smooth increase with increasing temperature. The value at 
300 K is 0.5 mS and it increases to about 2 mS at the 


highest temperature. Figure 40 shows the variation of 
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Table 3 


Relaxation Data and Results for 14, in HMPA 


Temp. T, obs. ee obs. ties Calc. 
(K) (msec) (x1012 sec) (x1012 sec) 
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Figure 40. 


Temperature (K) 


Observed relaxation rates of 14y in HMPA. The 


1?) -relaxatiron times were measured using the 
inversion-recovery pulse sequence. 


the 14n relaxation rate R,; with temperature. The smooth 
decrease of R, with temperature indicates that the 
quadrupolar relaxation mechanism dominates all other 
processes, as is typical for nuclei with quadrupole 
moments. 

Since the quadrupolar relaxation process entirely 
dominates other relaxation mechanisms, the 14, a data 
provide a very convenient means for determination of the 
reorientational correlation time toe The temperature 
dependence of Ry is) due tosthat or 4g (Equation {[5-12]). 
tg Can be easily determined from R, data if the nuclear 
quadrupolar coupling constant (e2q0/n) and the asymmetry 
parameter (yn) are known. The quadrupole coupling constant 
and the asymmetry parameter were determined from the 14y 
nuclear quadrupole resonance data reported by Krause and 
Whitehead (125). They have reported the temperature 
dependence of 14y quadrupole resonance frequencies of HMPA 
over the temperature range from 77 K to 200 K. The 
quadrupole coupling constant and the asymmetry parameter 
were calculated at the melting point of HMPA using their 
relationships with observed nuclear quadrupole resonance 
frequencies (125). Two sets of 14x quadrupole resonance 
frequencies are reported (125). The separation between 
these frequencies decreases with increasing temperature, 


but the lines do not coalesce below the melting point. It 


224 


finan sae switch Peers iP 
wits ake Ramin anasaneied et (Re 
wilt Fle Sonshine me inasogn walt 


ab eqhete iim Letuaen 169 contaggd ba 


Ylevb tee ceantdg, dGhgeselus Vebequsbaap oat aaet a 
nied © oP ® ante amet an dew colteehe galso sbseninee 
aie “)! ot terior ee 28 eger In6 (peviie yauy 8 sh ive 
cep sso bn, ant qs ertio fei el esten “Lane! oizueq 
(jh2-3) wliddps) , plo ani) of eal) sip 2a sonemnagen 
‘hulaun self al etab/ge ate? Gua leaatet ylieee a NBS gry 
criuamyue only Bibs Goghed “sastenee emgage vatoguahiog 
ineiados usiluuey aeyudbadp Set ~.AaQnd ata 1a) duane . 
wf! gfe evel hetlineedd® Saew 2eeegetng cetomey ye Be 
iw meusie ya beieeqes 2260 6s0Ronees el. tapribens sete 
LvaiAgees s0F HeIneges ever” CwrT ffi) beedes2aw 
he te catassepe? eonepamsa Glbguataup “Pt. te 
Lt OS OF NEY meee Byee Srutesegmss 
“addmet ng ifis amps Ad) Bnew Mae anes bind que. 
bso Qube MW to tateg gotal on oa jo f 
Wienkes: ebugualiaiy sewh aunt ; 
arenowed eheqwttnap ah! he a. 
saew2ed nO Mhingse war) 


io aes ns _ 
a “x 


‘ : ane | ra 
wo 7 


was therefore concluded (125) that the frequency 
separation was due to crystallographically inequivalent 
nitrogen atoms, there being two inequivalent nitrogen 
lattice sités, per unit cell an HMPA. ‘The separation, of 
resonance lines introduces some uncertainty into the 
calculations reported here. The average values of 
quadrupole coupling constant and asymmetry parameter 
evaluated from two sets of frequency data, were 
calculated. The average value of the quadrupole coupling 
constant e2Qq/h was Calculated to be 4.8 + 0.2 MHz and the 
asymmetry, parameter, 7 tO be 0.13 -£ G2. OL. 

The calculated values for the quadrupole coupling 
constant and the asymmetry parameter, along with the 
experimental R, values, were used in Equation [5-12] to 
obtain the reorientational correlation time ue These 


N 


values. are included yin Table 3.) The ue values for various 


temperatures were fit into the equation 
ve = a) exp(B>/T) ' [7-7] 


where Ay and By are fitting parameters. The non-linear 
least squares analysis yielded: A» = Barnicle topes 
sec and Bo = 1360 + 20 K. 

It is now appropriate to attempt a comparison of the 


observed variation of v with tz, with that predicted 
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using theoretical models for molecular reorientation. The 
relative magnitudes of the observed v5 (from 14n data) and 


Ty (from 31 


P data in HMPD) are such that v >> tz and 
hence these values correspond to rotational motion in the 
Debye limit (79,90). In this limit, the extended J — 
diffusion model (77,95-103) and the rotational Fokker- 
Planck-Langevin (FPL) model (78) give results which are 
identical to those given by Debye diffusion model (90). 
The investigated range of cand Uy ecoess Mel extend. 2ito 
the transition region between the Debye limit and the 
dilute gas limit. It is in this transition region where a 
rigorous test of the validity of the extended J-diffusion 


and FPL models could be made. Hence the 4N 


andmr, results 
Obtained in this work donot afford’ an opportunity to test 
models for rotational motion in) liquids. 


The experimentally determined variation of v, 


and ty 
has been compared with the variation predicted by the 
Hubbard relation (Equation [5-14]). The plots of the 
correlation times are shown in Figure 41. It is evident 


from Figure 41 that the observed variation of v 


Wi that 7 
does not agree with the theoretical variation. The 
Hubbard relation predicts substantially larger values for 
Tg than those experimentally observed and the deviation of 


the experimental tg values from the theoretical values 


becomes larger at higher temperatures. These deviations 
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are much larger than the uncertainties in v and ty; due to 
experimental errors in the T,; measurements and the 
uncertainties in the 4n quadrupole constant and the 31p 
spin rotation tensor. The increasing deviation with 
increasing temperature (larger TS values) suggests that 
the internal rotation process about the P-N bond must be 
taken into account in the theoretical model. Although the 
internal rotation about the P-N bond may be expected to 
possess a relatively large activation energy, it could 
affect the efficiency of 144 relaxation in the high 
TCemperatuLe range’. 


The reorientational correlation time +1 


deduced from 
the !4n data (Equation [7-7], Table 3) is therefore an 
effective correlation time (ci ee) which is a function of 
the correlation time for molecular tumbling (tg) and that 
for internal rotation about the P-N bond (TEN )e The 
internal motion about the P-N bond is superimposed on the 
overall molecular tumbling. It is possible to isolate the 
internal motion about the P-N bond if it is assumed that 
the molecular reorientation of HMPA is described by the 
Hubbard relation (Equation [5-14]). The effective 


correlation time alee for the case of isotropic 


reorientation is given by (73,126) 
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TeEE = Atg + Brg ti N/( tatty) 


+ Ctgti/(4tgtte) ; (78) 


where A, B and C are defined in Equation [5-8]. The angle 
A in Equation [5-8] in the case of quadrupolar relaxation 
corresponds to the angle between the internal rotation 
axis and the principal field gradient axis. The principal 
field gradient axis is assumed to lie along the nitrogen 
atom-lone pair direction and hence A is expected to be 
close to the tetrahedral angle. Equation [7-8] assumes 
that the internal motion and overall molecular tumbling 
are independent. The temperature dependence of the 
internal motion correlation time can be described by an 


equation of form 
aa = Ag exp (B3/T) eek 


where A3 and B3 are fitting parameters. The parameters A3 
and Bz were determined by fitting the observed 14y 
correlation times (Table 3) to Equation [7-8], and Tg was 
calculated from Hubbard relation (Equation [5-14]) with Ty 
given by Equation [7-6]. The least squares analysis gave 


the following results: Mee Sis 0.2 x 10713 sec and 
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B3 = 1710 + 80 K. Equation [7-9] shows the characteristic 
Arrhenius type variation and hence the energy barrier to 
internal rotation about the P-N bond can be extracted from 
it. The parameter By corresponds to AE/R where AE is the 
internal rotation barrier and R is the gas constant. 

Hence the barrier to internal rotation has been calculated 
to be 3.4 + 0.2 kcal/mole. No data on activation energies 
for the hindered rotation about P-N bonds are available in 
the literature for HMPA in liquid state, but P-N rotation 
barrier in solid HMPA has been estimated to be 6 kcal/mole 
by Andreeva et al. (127) from studies of 14 nuclear 
quadrupole resonance. The energy barrier in liquid HMPA 
can be expected to be lower than the value estimated in 
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Vir -H T, Results and Analysis 


thei) relaxation times (Of 2H nuclei in the 
deuterated analog (HMPD) were studied as a function of 
temperature. The 2H T; values were measured using the 
inversion recovery sequence over the temperature range 
250-410 K. The relaxation time at 250 K was 0.40 sec and 
it increased to 3.3 sec at 410 K. The Experimental 
results are shown in Table 4. Figure 42 shows the 


variation of the relaxation rate R, with temperature. 
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Figure 42. Observed relaxation rates of 5 in HMPD. The 
T,; relaxation times were measured with the 


inversion-recovery pulse sequence. 


The 2H relaxation in HMPD can be assumed to be 
entirely due to motional modulation of quadrupolar 
interactions. The quadrupolar relaxation rate of 
deuterium in the CD3 group in HMPD is governed by three 
correlation times: (1) that due to molecular tumbling 
Cte) (2) that due to internal rotation about P-N bond 


(+ and :03) that. due (to vinterna l rotation, soouk the CN 


ial 
bond ener The internal rotation about the P-N bond is 
Superimposed on the molecular tumbling motion to yield the 
effective correlation time weg discussed in Section 

Vil woeseeLne internal motion of stne CD, group is 
Superimposed on wee to dive a third) corre lation time tee 
which governs the time dependence of the 24 quadrupolar 


interaction. Since all three motions are independent of 


each other, tee can be written as 


= 1_N N N N 
weg = Attlec + Biche eco /( the ett$%) 


+ Chae cON/( 4c N ects $%) 


' [7-10] 
where tee is given by Equation [7-8] and A‘, B* “and C” 
are analogous to A, B and C respectively, in Equation 
(S-o)).. Gn. the debinitions of Al Bb” and ¢ 7 the angle 
variable is denoted by A‘ and is the angle between the 


internal rotation axis (the C-N bond) and principal axis 
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of the deuterium field gradient tensor (the C-D bond). 
Hence?Ae isutheitetrahedral jangle 1092028 ‘aneThe 
relaxation of the 2H nuclei in HMPD thus provides an 
interesting case where the relaxation interaction is 
modulated by a time dependence characterized by three 
correlation times. 

The quadrupolar relaxation rate in the presence of 
internal motion is obtained from Equation [5-12] by 


replacing tg by pares The deuterium quadrupole coupling 


in HMPD has not been reported, so we have assumed that the 


observed value, 166 + 10 KHz (128), for the CD3 deuterium 
quadrupole coupling constant in methylamine would provide 
an accurate estimate for the HMPD molecule. The value of 
e2Qq/h for D in CD3 group varies from 180 + 5 KHz in CD3I 
(i2°o) (Lom 33h J WKiz ei nechel (130) gendethegveriarioumis 
related to the electronegativity of the group attached to 
the CD3 group (131). Hence CD3NHy seemed the logical 
molecule to use in estimating e“Qqg/h for D in HMPD. The 
asymmetry parameter n for deuterium in CD3 group is very 
close to zero due to the cylindrical symmetry of the 


electric fileld gradient and* cansbe; dropped» torwa good 


approximation (132). Hence the relaxation rate expression 


(Equation [5-12]) reduces to 


poze eee Pee 
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The deuterium T, data (Table 4) were analyzed to 
yield the effective correlation time wee which is also 
given in Table 4. 

ine .errective correlation time wee appearing in 
Equation [7-10] was determined from the !4N data in 
Section VII.3. The temperature dependence of the internal 


rotation correlation time was assumed to have the form 


ota —SAASeXDA Ey T) ay [7-12] 
where Ay and By are fitting parameters. §1he observed 


deuterium correlation times (tees) (Equation [7-11]) were 
fit to Equation [7-10] using Equations [7-8], [7-91], 
[7-12] and the Hubbard model. Fixing A, and Bs at the 
values obtained from least squares analysis (Section 
VEIl@o)s the best. fit values torn Ajwand BZ were 
determined. The non-linear least squares analysis gave 
the following results 25) Anes Sil tea x 10714 sec and 
Bye 0082 40s eine DarlLerytomtne= yOorationvol the Cb. 
group in HMPD was determined from the parameter By to be 
3.0 + 0.1 kcal/mole. This compares favourably with the 
set of methyl rotation barrier values for various organic 
molecules compiled by Lambert et al. (81), for example, 


2.9 kcal/mole in 1,1,1-trichloroethane [CH3CCl3] and 3.5 
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kcal/mole in tertiarybutylchloride C(CH3)3cCc1] 
(133,134). Andreeva et al. (127) have reported the CH3 
rotational barrier in solid HMPA to be 2 kcal/mole, but 


have stressed that their work gives a lower bound. 


vir.5 1p 7, (HMPA) Results and Analysis 


The longitudinal relaxation times of 3lp in HMPA were 
measured over the temperature range 299 K to 460 K. Ty 
measurements were carried out using the inversion recovery 
sequence described in Chapter VI. The experimental 
results are shown in Pagure 43 and Table 5. The T; value 
at 296 K was 18 sec which increased to a maximum of 20 sec 
at 311 K and then decreased steadily to 10 sec at 461 K. 
The occurence of a maximum in the temperature variation of 
T, indicates the presence of two competing relaxation 
mechanisms: the spin rotation and dipolar contributions, 
which have complimentary temperature dependence. The 


total relaxation rate can be written as 


Ry = Riser + RiintrapDpD + RiinterDD ' ENS 


where Ryspr- RiintraDD and Ryinterpp 2° given by Equations 
[5-11], [5-8] and [5-9] respectively. The spin rotation 


relaxation rate increases with increasing temperature 
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Figure 43. 
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d 3'p relaxation 

The asterisk points show the observed T, 
rates with +10% error bounds. fThe solid lines 
calculated relaxation rates: curve A 

the spin-rotation relaxation rate; curve B 

the sum of spin-rotation and intramolecular 


Gipolar relaxation rates; and curve C represents the sum 


of spin-rotation, 


intramolecular and intermolecular 


dipolar relaxation rates. 
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Table 5 
Relaxation Data for ?!p in HMPA 
Temp. (K) eGcec) 
295.0 18 
295.6 18 
Be 20 
Vuilvs 20 
340.2 19 
3:5 Obeez 18 
SOM Ss) 18 
384.5 16 
AOOR 2 L5 
409.9 14 
420.4 14 
430.4 13 
441.2 12 
450.5 11 


460.9 10 
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while the intramolecular and intermolecular dipolar 
relaxation rates decrease with increasing temperature, but 
at different rates. An attempt was made to delineate the 
contributions due to spin=rotation, inter= and 
intramolecular dipolar interactions over the temperature 
range ‘studied: “To attain this objective, the individual 
relaxation rates were computed using parameters determined 
from the analysis of 14n, 24 and 3!p in HMpp data in 
previous sections. The correlation times weg and +N 
obtained from the !4n and 2H data were used to compute the 
Riintra term. The angular momentum correlation time oF 
determined from 3!p in HMpD was used to calculate RisR° 
The translational diffusion coefficients which were 
measured over the entire temperature range were used in 


the calculation of Riinter: 


a. Calculation of Spin-Rotation Relaxation Rates 

Equation [5-11] was employed to calculate the spin- 
rotational contribution to the total relaxation) rate. “The 
moments of inertia were calculated as described in Section 
Vit.4 using .crystallographicudata (120). ¥eThe=tellowing 
moments of inertia were obtained: I, Sek Guat COO. Tone! 
g cm? and Lh Orel 280. 01ex 10722 g cm2. The components 
of the spin-rotation tensor were evaluated using Equation 


{7=2], with®the results: (OF ja ie =) 999C+ 0 302 SKHzeand 
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Cy = -1.74 + 0.02 KHz. The angular momentum correlation 
time t; in HMPA was assumed to be equal to that in HMPD. 
This should prove to be a reasonable assumption since 
there moments of inertiawaresquite simijar.. ihe 
contribution Riga to the total relaxation rate was then 
calculated using Equation [5-11]. The results are given 


ithetabp ler 6c 


Db. Calculation of Intermolecular Dipolar Relaxation Rates 
The translational diffusion coefficient D, of HMPA 
was measured using the stationary field gradient technique 
described in Chapter VI. The experimental results are 
given, in Table 7. The diffusion coefficient at 304 Kk 1s 
3.8 x 107© cm4sec~! and it increases to 2.3 x 107° 
cm2sec7! at 450 K. The experimental data were fit to the 


equation 
De =nAc exp(B./T). , [7-14] 


where the fitting parameters A, and Bs were determined to 
be: As = 1.2 + 0.2 x 107% em#sec™! and Be = -1740 + 50 K. 
Equation [5-9] was used to compute Rj interpp: Thee tacvor 
a in Equation [5-9] represents the hydrodynamic radius of 
HMPA, i.e. the effective radius of the HMPA molecule if 


the molecule is assumed to be spherical in shape and to be 
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immersed in a continuous viscous medium. The hydrodynamic 
radius was taken to be 4 x 1078 cm. The calculated 


Riinterpp Values are given in Table 6. 


ec.) Calculation of the Intramolecular Dipole Relaxation 


Rates 

The intramolecular dipolar relaxation rates were 
calculated taking into account the effects of internal 
motion of methyl groups superimposed on the overall 
molecular reorientation. Temas been shown (135) that 
when internal rotation modulates the internuclear 
separation, the intramolecular dipolar relaxation rates 
can be affected in certain molecular geometries. Hence in 
such cases, Equation [5-8] which does not take into 
account the modulation of internuclear distance, does not 
adequately determine the relaxation rate. In the work 
(135) mentioned above, the internuclear distance term 
appearing in the correlation functions, was expressed as a 
function of the internal rotation angle and the 
correlation functions were evaluated (135) assuming that a 
rotational diffusion model describes both internal and 
overall molecular rotation. The spectral density 
functions resulting from these correlation functions were 
substituted in Equation [5-5] to give the relaxation rate 


expression [Equation [11] of Reference (135)] for dipolar 
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relaxation with modulation of the internuclear distance by 
internal rotation. 

In HMPA, the 34p-ly internuclear separation is 
modulated by rotation of the methyl groups. Hence the 
approach described above was employed for the calculation 
Of Riintrapp: The effective reorientational correlation 
time GN ee determined from the 14n data in Section VII.3 
and the internal rotation correlation time GN obtained 
from the analysis of the 24 data in Section VII.4, were 
substituted in Equation [11] of Reference (135) to 
evaluate the Rint ra Domes Dee Lesuits sare igivon in 


Table 6. 


ahs Comparison of Calculated and Experimental Relaxation 


Rates 

The calculated values for spin-rotation, inter- and 
intramolecular dipolar relaxation rates along with the 
experimental relaxation rates of 3lp in HMPA are given in 
Table 6. Graphical illustration of the data in Table 6 is 
shown in Figure 43. At 311 K, the spin-rotation mechanism 
contributes approximately 48% of the total relaxation rate 
while the inter- and intramolecular dipolar rates are 
approximately 41% and 19% respectively. Although the 
calculated relaxation rate is only grossly consistent with 


the experimental value at 311 K, it exhibits a reassuring 
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feature: the contribution of dipolar relaxation 
mechanisms and spin-rotation mechanisms are approximately 
equal at the minimum in the plot of relaxation rate with 
temperature (136). The calculations (Table 6) show that 
Une <spln=rotation contribution iss48¢ "ar sll Kewhich is 
approximately the temperature at which the experimental 
plot shows a minimum (Figure 43). The spin-rotation 
contribution increases to about 80% at 400 K while the 
inter- and intramolecular dipolar contributions decrease 
to about 9% and 5% respectively. At 461 K the spin- 
rotation interaction contributes approximately 81% of the 
total relaxation rate and the dipolar mechanisms account 
for about 6% of the total rate. Although the calculated 
relaxation rates do not agree with experimental rates in a 
rigorous quantitative manner, the delineation of the 
contributions of various relaxation mechanisms remains 
valid considering the limits of experimental errors and 
the various approximations used. 

Another interesting feature which emerges from Table 
6 is that the intermolecular dipolar relaxation rates are 
larger than the intramolecular dipolar rates. Although 
this result seems surprising, it can be rationalized when 
the implications of the internal rotation of methyl groups 
are recognized. It has been shown (137) that the effect 


of internal motion is most pronounced when the internal 
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rotation correlation time Ten and the overall 
reorientational correlation time Tg are comparable. In 
HMPA, co and ty correlation times are comparable 
(Equations [7-7] and [7-12]) and hence the internal motion 
can be expected to alter the intramolecular dipolar 
contribution when compared with situations where internal 
motion is absent or when met iewdrt ferenterrom 4 by 


several orders of magnitude. 


Vir ov uconclustitons 


The work described in Chapters V, VI and VII of this 
thesis demonstrates the use of relaxation data for the 
various nuclei to obtain an overall picture of the 
molecular motion in HMPA. 

The study of 3lp relaxation in the deuterated 
compound (HMPD) has yielded the angular momentum 
correlation times 77. The) study of the deuterated form of 
the molecule has been useful in isolating the spin- 
POtatloneinteract ton, Ore tile 31p nucleus. 

From the 14n data analysis, it has been found that 
the effective correlation time for reorientation is a 
function of the molecular tumbling correlation time and 
that for internal motion about the P-N bond. It was found 


that the inclusion of the internal motion about the P-N 
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bond was essential for agreement between the experimental 


Varlatwron vof a 


with ty, and the corresponding variation 
predicted by the Hubbard relation. The internal rotation 
about the P-N bond was found to be quite hindered, with an 
acti Vationivenergy .of (34 = 042 kKeal/moles, The internal 
rotation correlation time PN was about 5 times larger 
than the correlation time for molecular tumbling tg (given 
by the Hubbard relation) at room temperature. At 450 K 
this ratio decreased to about 1.4. Thus at high 
temperatures the internal motion time scale was comparable 
to that for molecular tumbling and hence altered the 
effective reorientation correlation time Nee 
significantly. 

The relative magnitude of Tg and tz was such that 
Tg >> tz, which is the Debye limiting situations = Insthe 
Debye limit, the results from various models for molecular 
reorientation, such as the J-diffusion model (77,95-103) 
and the Fokker-Planck-Langevin model (78), essentially 
reduce to those given by the Hubbard relation. The 
applicability of models for molecular reorientation can 
only ‘be tested using the results at the transition stage 
between the Debye limit and the dilute gas limit. Hence 
it was not possible to test the validity of the 


J-diffusion or FPL models from the results on HMPA 


obtained in this work. 
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the-srelaxatiton study ot 2H in HMPD provided an 
interesting situation where the effective correlation time 
was affected by three types of motion: (1) overall 
molecular tumbling, (2) internal motion about the N-P bond 
and (3) internal motion of the CD3 groups. It was found 
that the rotation of CD3 groups was also relatively 
hindered ‘with an activation energy or 3.0 + 0.1 
keal/mol. This result seems to be in the correct range 
when compared with the results in the literature for 
activation energies for rotation of methyl groups. A 
comparison of the relative magnitudes of the correlation 
times, Give the ratio of Cotas to be 40-1531 at room 
temperature. At 450 K this ratio becomes approximately 
30:20:1. The relative ratio of each internal correlation 
Lime \tOSua 1s) such that, the effect of weach internals motion 
process on the overall motion is appreciable. 

It has been possible to delineate the individual 
COntri butions to 31p relaxation in HMPA at various 
temperatures. The separation of the contributions was 
Carried out by calculating the relaxation rates due to 
various mechanisms, using results derived from the 
analysis of 315 in HMPD, 14, and “H data. It was found 
that at about 311 K the spin-rotation and the total 
dipolar contributions were approximately equal at higher 


temperatures, the spin-rotation contribution became much 
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more dominant and accounted for most of the relaxation 
rate at 460 K. Another interesting observation was that 
the contribution of intermolecular dipolar interactions 
was larger than the intramolecular dipolar component. 

This can be rationalized on the basis that the internal 
motion of methyl groups, with correlation time comparable 
to that for molecular tumbling, reduces the effective 
Spectral density and hence makes the intramolecular 
dipolar relaxation rate smaller than one would find in the 


absence of internal rotation. 


Vib teed Summary 


The results of the relaxation rate studies in HMPA 
and the data analysis were presented in this chapter. The 
use of 31p Ty values in HMPD to obtain angular momentum 
correlation time was described. The analysis of 144, data 
to yield the correlation times for reorientation and 
internal motion about the N-P bond was presented. It was 
followed by 2H data analysis yielding the correlation time 
for internal rotation of methyl groups. The delineation 
of contributions to the relaxation rate of 3lp in HMPA was 
discussed. Finally, the conclusions derived from this 


study were presented. 
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